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SUMMARY 


This report contains a detailed description of the Langley computer 
program PLOTWD which plots and fairs experimental wind-tunnel data. The 
program was written for use primarily on the Langley CDC computer and CALCOMP 
plotters. The fundamental operating features of the program are that the 
input data are read and written to a random-access file for use during program 
execution, that the data for a selected run can be sorted and edited to delete 
duplicate points, and that the data can be plotted and faired using tension 
splines, least-squares polynomial, or least-squares cubic-spline curves. The 
most noteworthy feature of the program is the simplicity of the user-supplied 
input requirements. Several subroutines are also included that can be used to 
draw grid lines, zero lines, axis scale values and labels, and legends. A 
detailed description of the program operational features and each subprogram 
are presented. The general application of the program is also discussed 
together with the input and output for two typical plot types. A listing of 
the program code, user-guide, and output description are presented in appen- 
dices. The program has been in use at Langley for several years and has 
proven to be both easy to use and versatile. 
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INTRODUCTION 


Although the primary function of the existing NASA Langley wind-tunnel 
complex is the same as that of its earlier predecessor NACA, the capabilities 
of the existing facilities are far superior to those of its predecessor with 
regard to the relative speeds obtainable, Reynolds number range, and the types 
and quality of data that can be taken. The advent of modern strain-gage 
balances, pressure and acoustic transducers, flow visualization systems, 
engine simulators, attitude measuring devices, and tunnel control and data 
acquisition systems have resulted in a tremendous increase in the amount of 
data generated during a typical test program. The recent "explosion" in mini- 
computer technology has further enhanced the ability to acquire even more 
data. Without the availability of both on- and off-line high-speed digital 
and analog computers, the test engineers' ability to acquire, reduce, and 
analyze this large volume of data would be totally inadequate. 

Following a typical wind-tunnel test program, the test engineer is faced 
with a large amount of data that must be plotted and faired for analysis and 
report presentation. During the earlier NACA period, these data were plotted 
and faired manually which required the services of a rather large supporting 
staff. This manual operation was slow, very tedious and monotonous, and 
often very inaccurate. However, since the introduction of high-speed com- 
puters with peripheral mechanical and optical plotting devices, manual 
plotting is no longer necessary or desirable. The Langley computing complex 
currently consists of some of the advanced CDC computers and CALCOMP plotting 
devices available. The speed, quality, and accuracy of the plots generated by 
the CALCOMP devices are generally far superior to those generated manually. 

The available CALCOMP computer software package is extremely versatile and 
relatively easy to use. 


2 



The purpose of this report is to describe the features and application of 
a computer program written specifically to plot and fair wind-tunnel generated 
data. This program has the internal designation "PLOTWD" and was written 
strictly for use on the Langley CDC computer and CALCOMP plotting systems. 

The program can, however, be adapted to other computer and plotting systems 
provided the computer has random-access file capability and the plotters have 
comparable system software routines. This computer program has been used 
extensively by personnel in the Langley 4— by 7— Meter and 7— by 10-Foot High- 
Speed wind tunnels for several years and has proven to be very reliable, 
flexible, and easy to use. In fact, one of the most noteworthy features of 
this program is the simplicity of the input requirements. 

This report contains a discussion of the overall program capabilities and 
operational structure. A detailed description is included of each subprogram 
contained in the program. A listing of the computer code is presented in 
appendix A. Only the subprograms that were written specifically for this pro- 
gram are listed in the appendix. The external subroutines called that are 
part of the Langley CALCOMP software package are not listed. A description of 
the user-supplied input requirements and a description of a sample output are 
presented in appendixes B and C, respectively. 
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SYMBOLS 


a f a 2 , . •/ a^ polynomial coefficients 

c d , e , f . coefficients of spline equation 

i' i i i 


increment in spline independent variable such as 
(x - x i ), (y - y i )/ or (t - t^ 
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Subscripts; 

o 

P 

sf 

wt 


order of polynomial 

number of points 

polynomial equation 

sum-of- squares of differences 

cubic-spline smoothing parameter 

weighting factor 

independent variable 

x— axis dependent variable 

d 2 x/dt 2 

y— axis dependent variable 
d 2 y/dt 2 

spline tension factor 
increment 

origin 

plotting sheet 
scale factor 
wind-tunnel data 
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PROGRAM OPERATIONAL FEATURES 


The computer program PLOTWD was originally formulated to plot and fair 
sequentially obtained wind-tunnel data on the large 33-inch CALCOMP plotter 
paper with either 10 or 20 grid lines per inch. The plots generated could be 
used either as working figures or, after proper border masking and labeling, 
as final report figures. The later change in the Langley editorial figure 
presentation standards due to the enhancement of the CALCOMP software package 
allows for the acceptance of final report figures that are almost totally 
generated with computer graphics. As a result, additional subroutines were 
added to the original program version that can be used to draw grid lines, 
axis scale values and labels, and a legend containing the run numbers and 
corresponding plot symbols. These plots can be drawn on the plain white 
CALCOMP paper and, after additional legend information has been added, are 
suitable report quality figures. 

The data on a figure for a particular run are identified by one of 22 
available symbols as illustrated in Table I. Each symbol can be drawn with 
three basic sizes: (1) small (0.100 in.), (2) medium (0.132 in.), and (3) 

large (0.168 in.). The data are faired using the standard curve fairing 
technique which utilizes tens ion- splines to compute a maximum of 100 enhance- 
ment points between each pair of input data points. All enhancement points 
are plotted as a continuous solid line which intersects but does not pass 
through the symbol surrounding each input data point. During a typical wind- 
tunnel test, the engineer will often request that repeat data points be taken 

to either establish anomalies in performance trends or to provide a routine 

„ Loofe [ 

check value. '("Fitting a spline curve through the input data with the repeat 
points included will almost always produce curves with unrealistic oscilla- 
tions, especially between original and repeat data points^ To overcome this 
problem, a subroutine has been included in the plot program which first sorts 
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the data so that the data points are in monotonic increasing order and then 
deletes all but one data point in each cluster of duplicate data points. The 
deleted points are not used as input to the spline curve fit subroutines, but 
they are plotted as individual data points with the appropriate symbol. 

Additional subroutines are available in the plot program that can compute 
and plot either least-squares polynomial or least-squares cubic-spline curve 
fits to the input data. These least- squares curve fits are plotted with a 
continuous dashed line corresponding to the particular symbol as illustrated 
in Table I. If the user desires, a group of subroutines are available to draw 
grid lines, axis scale values and labels, and legends. The axis labels that 
can be drawn consist of a set of nine standard aerodynamic coefficients as 
illustrated in Table II. A convenient place has been designated in the sub- 
routine that draws labels where the user can add the coding needed to draw 
other axis labels. The character sets available in the CALCOMP software pack- 
age are very versatile and the user should be able to draw almost any imagin- 
able label. 

The user-supplied input requirements were formulated to be simple and at 
the same time provide as much flexibility as possible. This flexibility 
consists of the ability to plot multiple figures on a single sheet, to plot 
almost any combination of data variables, and to plot almost any combination 
of runs on a sheet. The input requirements for the basic type of plot without 
drawn grid lines, axis scale values and labels, and legends are divided into 
two parts. The first part, designated as the plotting setup information, 
identifies the plotting sheet boundaries, the symbol size and starting number, 
curve fairing and editing options, the data array location of the independent 
test variable, and the combinations of test variables to be plotted with cor- 
responding axis origin locations and scale factors. The second part, desig- 
nated as the plotting sheet information, is a set of namelists, each 


6 



containing a sheet identification number and a list of the runs to be plotted 
on the sheets. 

The wind-tunnel data to be plotted must be written on a file designated 
as TAPE1 and attached as part of the job execution procedure. The data must 
be written on TAPE1 with the following unformatted write statements 
WRITE ( 1 ) I RUN, ITEST, IPT, (DATA(I), 1=1, IMAX) 
where IRUN is the run number, ITEST is the test number, IPT is the data point 
number, DATA is an array containing the wind-tunnel data, and IMAX is the max- 
imum number of elements in the DATA array. This input data requirment should 
be able to accommodate almost any type of wind-tunnel data that is taken 
sequentially. After reading the setup information during program execution, a 
subroutine is called that reads the data on TAPE1 and reloads it on the 
random-access file TAPE2 with an equivalent run number indexing. Thereafter, 
anytime a particular run is to be plotted, the corresponding data for that run 
is read from the random-access file and loaded into a plotting array. 
Retrieving the data in this manner simply means that the file read of the data 
for a particular run begins immediately at the first data point rather than 
each time having to read the data file from the beginning to locate the first 
point. This random-access feature greatly reduces the total computing costs 
by reducing the execution time and the number of I/O operations. 

DESCRIPTION OF COMPUTER PROGRAM 

The wind-tunnel data plotting program PLOTWD consists of a main program, 
18 subroutines, and 2 function subprograms. A computer listing of the program 
coding is presented in Appendix A. A description of the input requirements for 
the program is presented in Appendix B and a description of the output for a 
sample input case listed in Table III is presented in Appendix C. The primary 
function and execution sequence for the main program and each subprogram are 
discussed in this section. A supplemental list and corresponding description 
of the input variables and internal parameters for each routine is included. 



Program PLOT 


The function of program PLOT is to control the overall execution of the 
plotting process. Listed at the beginning of the program are a group of com- 
ment statements that describe the program input requirements. A group of 
dimension and common statements appear next consisting of the arrays contain- 
ing the run numbers and plotting variables, axis origin locations, and axis 
scale and label information. Several global program constants are then 
defined and calls made to subroutines PSEUDO and LEROY to initialize the plot 
vector file SAVPLT for subsequent post-process plotting on a variety of 
plotters at Langley. A call is then made to subroutine SETUP which reads the 
input plotting setup data. This is followed by a call to subroutine DATADK 
which loads the input wind-tunnel data that is sequentially written on the 
input file TAPE1 into the random-access working file TAPE2 • This subroutine 
is called only once during the program execution? therefore, if the plotting 
setup information is changed during the program execution, the test number may 
not change. 

The data array locations of the independent and dependent variables to be 
plotted are then loaded into the work array NLOC. A check is then made to 
determine if any parts of the plotting setup have common abscissas or ordi- 
nates and, if so, the corresponding internal XPLT and YPLT array elements are 
set equal to 1.0. These internal arrays are checked during the axis labeling 
and zero-line drawing execution step to prevent duplication. The next execu- 
tion step is the read of the input namelist SHEET which contains the sheet 
number and corresponding run numbers to be plotted. The test, sheet, and run 
numbers are then drawn at the bottom of the plotting sheet and the sheet frame 
drawn. Arrowheads are then drawn along the bottom and left-hand edge of the 
sheet frame to indicate the origins. The corresponding axis scale and label 
are then drawn adjacent to each arrowhead. Zero lines which run from the 
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arrowhead to the opposite frame border are then drawn if the parameter LZERO 
is equal to zero. 

Execution then proceeds to plot the desired wind-tunnel data array 
variables for each run* The first step is an identification of the plotting 
symbol to use. The 22 symbols available are presented in Table I in the 
standard order. The number of the symbol to use is the same as the input 
order of the run number in the SHEET namelist. If the user desires to skip a 
particular symbol, simply input a run number of zero at the appropriate place 
in the run number list in the SHEET namelist. If the user desires to set up a 
symbol order different from the standard order, this can easily be accom- 
plished by redefining the sequential values in the internal array LSYM. 

The next step is to read the particular wind-tunnel data corresponding to 
the variable to be plotted from the random-access file TAPE2 and to load these 
data into the internal array VAR. The next step is to convert the wind-tunnel 
data- from cofficient form to the equivalent Xp and yp locations on the 
plotting sheet using equations 


x = x . /x - x 
p wt / sf o 


^p ^wt^sf ^o ' 


( 1 ) 


( 2 ) 


and to load the converted values into arrays X and Y, respectively. At the 
same time, the values for the independent data variable are loaded into array 
T and a check is made of the Xp and y p data to determine if any points are 
outside the sheet frame. If any points fall outside the sheet frame, a mes- 
sage will be printed stating the run and plot number and the number of points 
outside the frame. The next step is to sort and edit the Xp and y p data 
if the user specified a nonzero value for the input parameter IEDIT. The 
final step is to call subroutine CURPLT which plots and fairs the x p and 
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Yp data using tension splines. After plotting all of the specified variable 
combinations for a particular run, the above sequence of steps is repeated for 
the next run. 

After plotting all the specified runs for a particular sheet, the plotter 
is advanced to the start location for the next sheet, the next SHEET namelist 
read, and the entire sequence of scaling and plotting steps are repeated. The 
next execution step in the program is to print a summary of the sort and edit 
information. The final execution step is to draw a statement indicating that 
all plotting has been completed and to call subroutine CALPLT to close the 
plot vector file SAVPLT. The following is a list and description of the 
parameters used in this program: 

LSYM array containing integer values corresponding to the symbol 

order 

RUN array containing run numbers to be plotted 

NPT - array containing the directory or index information on the 

random-access file TAPE2 (dimensional size equivalent to 
number of data points on input wind-tunnel data file TAPEl) 
XTAPE, YTAPE arrays containing wind-tunnel data array indices of x wt 

and y wt variables 

XOFFSET, YOFFSET arrays containing origin x Q and y Q locations on plotting 
sheet, in. 

XSCALE, YSCALE arrays containing axis scale factors x s f and y s f# change 
in coefficient per inch 

XLABEL, YLABEL arrays containing x- and y-axis labels 

XPLT, YPLT arrays containing a value of 1.0 if the particular x- or 

y-axis is repeated 

X,Y arrays containing scale Xp and yp values, in. 

T array containing value of independent variable 
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VAR 

NLOC 

NO 

NEWCASE 

JREAD 

NRNMAX 

NPLMAX 

NDMAX 

NPMAX 

IPRT 

LZERO 

NORG 

IDISK 

ITEST 

SHEETW, 

SPACE 

ISYM 

ISIZE 

IOP 

TENSION 


two-dimensional array containing y wt data in one level and 
x w t data in the next level 

array containing consecutive DATA array indices of y^ and 
x wt variables 
sheet number 

if = 1, new setup deck follows the current SHEET namelist 
if = 0, new SHEET namelist follows 

tape number of file containing input setup and namelist data 

maximum number of allowable RUNS per sheet 

maximum number of allowable plots per sheet 

maximum number of allowable data points per run 

maximum number of total data points on TAPE1 

if = 0, print sort - and- edit summary data 

if = 1, do not print sort — and— edit summary data 

if = 0, draw zero lines 

if = 1, do not draw zero lines 

if = 0, draw arrowheads at axis origins 

if = 1, do not draw arrowhead at axis origins 

if = 0, random-access file TAPE2 has not been loaded 

if = 1 , TAPE2 has been loaded 

test number 

SHEETH plotting sheet width and height, in. 

space between plotting sheets, in. 

starting symbol number 

symbol size 

data fairing option 

spline tension factor a 
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IEDIT data sort-and-edit option 

TOLR editing tolerance of independent variable 

IW DATA array index of independent variable 

NPLOT number of desired plots per sheet 

IERR if = 1 after call to subroutine SETUP, indicates that last 

case has been read or error occurred 

if = 1 after call to subroutine DATADK indicates that error 
occurred during either a read or write of TAPE2 
NVAR number of variables read from TAPE2 for each data point 

(2*NPL0T+1) 

NP number of data points to be plotted for a specified run 

Subroutine SETUP 

The function of subroutine SETUP is to read and print the plotting setup 
information as described in Appendix B. A sample of the print formats are 
presented as page 1 output in Table IV. As various parameters are read from 
the input file, checks are performed to insure that the parameter is within 
the program limitations. If the particular parameter cannot be safely 
redefined within the program limitations, an error message is printed and pro- 
gram execution terminated. All the parameters in the subroutine argument list 
are defined in the description of program PLOT. 

Subroutine DATADK 

The primary function of subroutine DATADK is to read data from or write 
data to the random-access file TAPE2 • If the ICODE parameter equals 0, the 
wind-tunnel data on the input file TAPE1 is read and then written on the 
random-access file using the mass storage write statement WRITMS. Each data 
point is read from TAPE1 using the following read statement: 

READ ( 1 ) IRUN, ITEST, IPT, ( DATA ( I ) , 1=1, I MAX) 
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The information in the data array is then shifted forward three array loca- 
tions and the values of IRUN, ITEST, and IPOINT loaded into the first three 
locations. The information in the resultant DATA array is then written on the 
random-access file TAPE2 and the access location loaded into the array NPT. 

This access location will be needed during all subsequent reads of this file; 
therefore, the user should not add any additional coding that would destroy 
the contents of this array. If the data being written on TAPE2 corresponds to 
that for the first point of each run, the run number and the array index in 
the NPT array are written on a local file TAPE12 for use during the later 
reads of TAPE2. The default value for the parameter IMAX is 300. If the 
input data tape contains either more or less variables for each data point, 
the value of IMAX must be redefined and the dimension size of the DATA array 
redefined to a value equal to IMAX + 3. After transferring the data on TAPE1 
to TAPE2, a summary of the run numbers available and the total number of data 
points transferred is printed as illustrated in Table IV as page 2 output. 

If the ICODE parameter equals 1, all the wind-tunnel data points for a 
specified run number IRUN are read from the random-access file TAPE 2 and the 
specified variables defined by the user-supplied input values for XTAPE , 

YTAPE , and IW are then loaded into the two-dimensional array VAR. The random- 
access index location for the first point of the specified run is determined 
by reading the run and NPT array index information written previously on 
TAPE12. During the read of each data point on TAPE2, the data are loaded into 
the array DATA. The subsequent transfer of the information in the DATA array 
to the VAR array is shifted back three array locations to prevent loading the 
values for IRUN, ITEST, and IPT into the VAR array. The test number read from 
TAPE2 is, however, checked against the test number specified by the user in 
the setup information to ensure that the correct test data were input on 
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Several types of errors can occur which may cause the program execution 
to either terminate or proceed to the next specified run number* The types of 
errors that may occur include (1) there are no data on TAPE1, (2) the total 
number of input data points exceeds the allowable maximum, (3) the input and 
specified test numbers are different, (4) the specified run number is not on 
TAPE2, (5) number of data points for a specified run number exceeds the allow- 
able maximum, (6) the number of data points on TAPE2 for a specified run can- 
not be determined from the information written on TAPE12, and (7) the speci- 
fied XTAPE or YTAPE value exceeds the allowable maximum IMAX. If any of these 
errors occur, an appropriate error message will be printed. Most of the 
important parameters used in this subroutine are defined in the description of 
program PLOT. The following additional parameters are used in this sub- 
routine: 

DATA array containing and y wt data read from either TAPE1 or 

TAPE2 (must be dimensioned by IMAX+3) 

IMAX maximum number of elements in data array 

IRUN specified run number 

IPT input data point number 

ICODE if = 0, read data on TAPE1 and write to random-access file TAPE2 

if = 1, read data on TAPE2 for a specified run and load into the 
VAR array 

OPENMS system subroutine to open and initialize random-access file 

READMS system subroutine to read information from random-access file 

WRITMS system subroutine to write information on random-access file 

Subroutine EDIT 

The function of subroutine EDIT is to sort and edit a set of x p and 
y p data to be plotted for a specified run. If the user desires to fair 
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the wind-tunnel data with the tension spline method, the set of data to be 
faired must be arranged so that the corresponding values for the independent 
variable, t, are monotonically increasing, single-valued, and reasonably 
spaced to prevent large oscillations in the faired curve. The first time this 
subroutine is called with a specified run number, the results of the sort and 
edit procedure are written on the output file TAPE? and the run number on out- 
put file TAPE11 for future printing during the call to subroutine PRTEDT • 
Therefore, the first execution step in this subroutine is a check of the 
information on TAPE1 1 to determine whether or not summary data have previously 
been written on TAPE?. If so, the parameter NSTORE is set equal to a value 

of 1. 

The next execution step is to sort the input values of the independent 
variable, t, so that the final set of values are monotonically increasing. 

The sorting is accomplished in an iterative manner with the following steps: 

(1) Load remaining unsorted values of the independent variable and 
their corresponding indices into temporary arrays TEDT and IEDT , 
respectively. 

(2) Determine the minimum TEDT value and its corresponding index and 
load index value into array ISAVE. 

(3) Delete the minimum TEDT value from the TEDT array. 

(4) Repeat steps 1, 2, and 3 until all input values of the independent 

variable have been processed. 

(5) Load sorted independent variable values and corresponding y p values 
into temporary arrays using index information in the ISAVE array. 

If the user specified on edit tolerance TOLR that is greater than zero, the 
sorted data will be edited next to delete all but one data point in each 
cluster of repeat points. Repeat points are defined as any two adjacent 
points whose difference in the absolute values of the independent variable 
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is less than a specified tolerance. The edit procedure consists simply of 
connecting the two data points adjacent to a particular cluster of repeat 
points with a straight line and then determining the point in the cluster 
closest to the line. The closest point will be saved and the remaining points 
in the cluster will be deleted. These deleted points are, however, plotted 
immediately on the plotting sheet with the appropriate symbol. 

Although this editing procedure is simple, four problems may be encoun- 
tered which will prevent the sorted data from being properly edited. If any 

of these particular problems occur, the error flag IERR will be set equal to 

the problem number, all data points input for the particular run plotted with 

the appropriate symbol, and control returned to the main program PLOT. Pro- 
gram execution will then proceed to plot and fair the next specified run 
number. The first error which may occur is that all the values for the sorted 
independent variable are within the specified tolerance TOLR. This error 
usually occurs because the user specified a tolerance which is too large. The 
remaining errors occur due to the presence of several adjacent clusters of 
repeat points which prevents formation of the straight line between adjacent 
data points. The second flagged error indicates that all data points except 
the last are clustered; the third flagged error, that all data points prior to 
the particular point being edited are clustered; and the fourth flagged error, 
that all data points except the first are clustered. If any of the latter 
three flagged errors are printed in the output, the user should carefully 
examine the data and decide which data points to delete. 

Following the editing procedure, a summary of the sort and edit informa- 
tion is written on TAPE7 provided the value of the NSTORE parameter is zero. 
The last execution step is to reload the independent variable t, x p , and 
Yp input arrays with the sorted and edited values. The following additional 
parameters are used in this subroutine: 
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KEDT 

TIN 

TEDT 

TSAVE, YSAVE 
IEDT, ITEMP, 
I SAVE 
TMIN 
NSTORE 

NSYM 


number of run numbers written on TAPE 11 
array containing input values of independent variable t 
array containing sorted and edited values of independent 
variable 

temporary arrays containing values of t and yp 

temporary arrays containing indices of sorted or edited t 

values 

minimum value of t 

if = 0, summary of sort and edit information for particular 

run is to be written on TAPE7 

plot symbol number for deleted points 


Subroutine PRTEDT 

The function of subroutine PRTEDT is to read the sort and edit 
information generated during the calls to subroutine EDIT# This information 
is read from TAPE7 and consists of the input, sorted, edited, and deleted 
values of the independent variable t for all of the runs requested during the 
job execution. A sample of the output is illustrated as page 4 output in 
Table IV. As previously noted, the sort information is the same regardless of 
the number of times the particular run is plotted. However, the edit informa- 
tion is only applicable to the first plot of the run. All of the important 
parameters used in this subroutine are defined in the description of subrou- 
tine EDIT. 

Subroutine ZEROLN 

The function of subroutine ZEROLN is to draw heavy lines on the plot 
perpendicular to the axis at each XOFFSET and YOFFSET location to indicate 
origins or zero lines for the variables to be plotted. Each heavy line is 
generated by simply drawing a straight line perpendicular to the sheet border 
at the specified origin location and then drawing two additional lines 
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parallel to the original line and offset to either side of it. The following 
subroutine arguments and parameters are used: 

Xo, YO starting Xp and yp plot locations of the zero line, in. 

DIST length of zero line, in. 

line if = 1, zero line parallel to x-axis 

if = 2, zero line parallel to y-axis 

D offset distance of additional lines drawn parallel to zero 

line (width of zero line = 2 x D + plotting pen width) , in. 

Subroutine CURPLT 

The function of subroutine CURPLT is to plot and fair a tension spline 
curve through an input set of Xp and yp values. The subroutine contains a 
plot option to plot the symbols only (IOP = 0) and a plot option to plot the 
symbols and draw a faired curve between them ( IOP = 1 ) • For both plot 
options, the symbols are plotted by calls to subroutine PNTPLT which has been 
attached as part of the CALCOMP software package. If the curve fairing option 
has been chosen, the input data points x p and y p are faired with a tension 
spline that has a tension factor equal to the value of the input parameter 
TENSION. The coordinates of the faired curve are determined by fitting a ten- 
sion spline through both the x p and y p as a function of the independent 
variable t and then interpolating the two curves to obtain the resultant y p 
as a function of x p faired curve which is then drawn on the plot. If the 
Xp is a scaler function of t (XTAPE value equals IW) , the double interpola- 
tion process is not necessary. The input value of the parameter IVAR must be 
equal to 0 if Xp and t are equal and equal to 1 if they are not equal. The 
first execution step is a check to ensure that the independent variable values 
are monotonically increasing. The next step is a call to subroutine CUBSPL 
which determines the value of the second derivative x p " and y p " of the 
tension spline at the input points. The final execution step is to 
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interpolate a specified number of spline coordinates between each pair of 
input points and to draw the curve so that it does not cross the borders of 
the symbol surrounding each data point. This capability represents a unique 

feature of this program. 

The following discussion outlines the procedure used to determine the 
intersection points of the spline curve with the border of each end point 
symbol and the number of coordinates to interpolate between the pair of 
points. The straight line distance between the pair of points is computed 
first. If the length of the line is less than the nominal width of the 
symbol, the spline curve cannot be drawn between the pair of symbols and 
execution advances to the next pair of points. If the length of the line is 
greater than the nominal width of the symbol, the interpolation increment of 
the independent variable, t, between the pair of points is determined using 


the formula 




(-S C$. 
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(3) 


(NH-1) I(t) 


where NH - 101 and I(t) is the integer portion of (t i+1 - . Spline 

coordinates are then computed at consecutive 6 increments between the pair of 

points. 

The next step is to determine the intersection points between the spline 
curve and the borders of the symbols drawn at each point. This step is accom- 
plished utilizing an iterative technique that is centered around the capabili- 
ties of subroutine SYMBOL which determines the coordinates of the border of a 
specified symbol as a function of the angular position of a radial line drawn 
through the symbol center. The iteration cycle begins at a point on the 
spline curve from the symbol center corresponding to a value of t equal to an 
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incremental value, H, based on the symbol size. The coordinates of the spline 
curve are interpolated using the function subprogram FUNC and then the 
straight line length and angular location computed from the center of the 
symbol. The angular value is then input to subroutine SYMBOL which determines 
the corresonding coordinates on the symbol border. The straight line length 
from the symbol center to the border is then computed. If the spline straight 
line length is less than that from the center to the border, the value of t is 
incremented by H and the above steps are repeated. If the spline straight 
line length is greater than that from the center to the border, the increment 
H is halved and added to the previous t value and the above steps repeated. 
This procedure is repeated until the difference between the spline and symbol- 
border straight line lengths is within a tolerance equal to the value of the 
parameter EPS. 

After determining the intersection coordinates of the spline curve and 
the symbol borders for a pair of points, the final step is to plot the spline 
curve between the intersection coordinates. The LIME subroutine, which is 
part of the CALCOMP software package, is used to plot this curve. This sub- 
routine draws a series of connecting straight lines between an input set of 
coordinates. By plotting a relatively large number of closely spaced coordi- 
nates, the smoothness of the curve between each pair of input points is 
retained. The following additional subroutine arguments and parameters are 
used in this subroutine: 

IVAR if = 0, independent variable t is plotted on x-axis 

if = 1 , a dependent variable is plotted on x-axis 
IOP if =» 0, plot symbols only 

if = 1, plot symbols and spline curve 
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MX, MY 


DS1 , DS2 


DUMX, DUMY 


NH 


EPS 


H 


arrays containing second derivative values x" and y^ 
of tension spline curve fit 

arrays containing interpolated Xp and yp spline curve 
coordinates between a pair of input points 
arrays containing interpolated Xp and yp spline 
coordinates between the borders of the symbols surrounding 
a pair of input points 

number of interpolated spline coordinates between a pair of 
input points 

tolerance used to determine the intersection point between 
each symbol border and the spline curve 
initial increment in t value 


Function ATANF 

The purpose of function subprogram ATANF is to compute the angular posi- 
tion of a specified coordinate point in degrees with the quadrants defined 
from 0 to 360 degrees. This function is used with subroutine CURPLT to deter 
mine the radial angle of a straight line connecting a specified point on the 
spline curve and the center of the nearest symbol. 


Subroutine SYMBOL 

The function of subroutine SYMBOL is to compute the coordinates of the 
intersection between the border of a specified symbol and a radial line drawn 
from the center of the symbol. This requires that the border of each symbol 
presented in Table I be defined as a function of the angular equivalent of 
border coordinates with the origin of the axis system at the symbol center. 

The shapes of 22 symbols are defined in this subroutine. The first 10 symbols 
are the NASA standard open symbols: (1) circle, (2) square, (3) diamond, 

(4) triangle, (5) right triangle, (6) quadrant, (7) dog house, (8) fan, (9) 
long diamond, and (10) house. The next 10 symbols are the standard centered 
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(center of symbol indicated by +) symbols corresponding to the standard open 
symbols* The 21st symbol is a period and the 22nd symbol is a plus sign. The 
symbol border for the last two symbols is an imaginary circle with a radius 10 
percent larger than the symbol height. The shapes of the 10 basic open and 
centered symbols are made up of various combinations of circles, squares, and 
triangles. Most of the symbols are symmetric about the local y-axis except 
the right triangle and the quadrant. The shape of each symbol is determined 
by computing the straight line or circular arc coordinates of the various 
segments that make up the symbol. The equations and angular positions of the 
segments that make up each of the 10 basic symbols are relatively easy to 
derive and, therefore, are not presented in this report. The following argu- 
ments and parameters are used in this subroutine: 

NO symbol number 

IS symbol size (1 for small, 2 for medium, and 3 for large) 

X, Y x and y coordinates of intersection of symbol border and 

j radial line drawn from the center of the symbol 
T angular position of radial line, deg. 

SCALE relative height of three basic symbol sizes, in. 


Subroutine CUBSPL 


The function of subroutine CUBSPL is to fit a cubic or tension spline 
through an input set of x and y values. The cubic spline represents a special 
case of the tension spline where the tension is zero. The input x and y data 
are used to compute a matrix of simultaneous equations in the form of 


+ d. y. + f.y.. i 

l l li+l 


(4) 


for i = 2, 3, • • ., N— 1 where for the cubic spline 




and f^ 



6 


(5) 
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and for the tension spline 


c . = 
1 


d. = 
l 


] 2 [ h i_i ’ sinh tah i-i’ ] 


acosh ( a‘n i _ 1 ^ 
sinh (ah i _ 1 ) h i-i 


and f . = 
l 


_l r_^ . 

2 h 

a L l 


sinh ( ah 


7-] 


For both cubic and tension spline, 

•. ■ (^) • ftr 1 ) 


acosh (ah^) 
sinh (ahj 


-d- 


(6) 


h i = X i+1 “ X i ' 


(7) 


and h._ 1 - x. - x i _ 1 . 

For the tension spline, the tension parameter a is the average tension per 
length of the input x values and determined using the equation 


a (N-1) (8) 

where a is the user-supplied value of the input parameter TENSION and N is the 
number of input values. A complete derivation and description of cubic and 

tension splines are presented in reference 1. 

The second derivatives are the unknowns in the system of equations 
generated by equation (4). The matrix of terms is tridiagonal with two less 
unknowns than equations; therefore, the second derivative at the end points 
must be specified. In this subroutine, the second derivatives at the end 
points are computed by fitting a second-order polynomial of the form 
y = a^ 2 + a 2 x + a 3 

to each end point and its two adjacent input points and then differentiating 
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to determine the second derivative which is 

y" = 2a t . (10) 

The resultant tridiagonal matrix can be solved using the Crout reduction 
method which is described in detail in reference 2» Using this method, the 
solution becomes a simple back substitution 
y N = C N for 1 = N 

(ID 

and y" = c . - f, y" for i = N-1, N-2, ... 1 

J i i i i+ 1 

where d . = d , - e . f , „ 

l i l i- 1 , 

f i - f ± /d i , (12) 

and 


c . — 
l 


c, - e.c. 
i i 1-1 


The following is a description of the parameters in the argument list for this 
subroutine: 


X,Y 


N 

YPP 

TENSION 

TENS 

A 


input arrays containing table of x- and y-values 
number of input coordinates 
output array containing y" values 
input tension factor a 

output dimensional ized tension factor o 

internal work array used during back-substitution process 


Function FUNC 

The function of subprogram FUNC is to compute the y-value of a spline 
curve at a given x-value in a specified interval x^ < x > x^ + -j. 
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The interpolation equation for the cubic spline is 


y(x) 



< X i+1 


- x) 
6 



+ 


i+1 


px - xj 3 (x - x ± ) h i “| 

4, 

r*i (x i+i • x) + y i+i <x ■ *i v 

[eh. 6 J 

T 

L J 


( 13 ) 


and for the tension spline is 



0 


“sinh [o(x i+1 - x)J 

1 

“x 

i 

+ 

-H 

X 

w* 

sinh (oh^) 

h. J 

1 


y'.' 

y i+1 

sinh 

a(x - x i >] 

i — 

•H 

X 

1 

>< 

4. 

’ y i (x i+i • x) + y i+i (x ‘ x i’l 

2 

0 

sinh 

(ahj 

h. 

1 J 

T 

L ^ J 


where h^ = (x - x^). 


Subroutine DASHLN 

The function of subroutine DASHLN is to draw a dashed line spline curve 
through an input set of x p and yp values. The spline curve fit procedure 
for this subroutine is very similar to that described for subroutine CURPLT 
with the major difference being that symbols are not drawn around each symbol 
with this subroutine. Therefore, the complex procedure to determine the 
intersection between the spline curve and the symbol border is not included in 
this subroutine. The dashed lines that can be drawn consist of various 
combinations of long and short dashes as illustrated in table I. 

The subroutine CUBSPL is called initially to determine the y” values of a 
tension or cubic spline fit of the input data. A specified number of equally 
spaced points are then interpolated between the input end points and the total 
length of the curve computed using a simple triangular summation method. The 
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basic length of the long and short dash combination and the space between 
combinations are then scaled so that the final plotted dashed line will start 
and end with a complete long and short dash combination. The scaled dashed 
line is then drawn using linear interpolation of the spline curve to determine 
the end points Xp and yp values for each consecutive long or short dash. 

The following additional subroutine arguments and parameters are used in this 
subroutine: 


NSYM 

IOP 

NL 

NS 

XI, YI 
S 

NP 

SL 

SS 

SP 

NLT 

DTN/DT 


symbol or dashed line number 

if = 0, plot symbols only 

if = 1, plot dashed line only 

number of long dashes in each combination 

number of short dashes in each combination 

arrays containing interpolated Xp and yp values of 

spline curve 

array containing interpolated length of spline curve 

number of interpolated spline curve values 

length of long dash, in. 

length of short dash, in. 

length of space between dashes, in. 

total number of long and short dash combinations to 

be plotted 

scaling factor for each dash combination 


Subroutine LINEAR 

The function of subroutine LINEAR is to determine the x and y values at a 
specified t value using linear (straight line) interpolation. The first 
execution step is to determine the interval containing the specified t value. 
The final step is to compute the corresponding x and y values using the 
following linear interpolation equations: 
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X 


y = 

where 6 = 


x. . + 6 

(X. - x._,) 


r-1 


y. - + 6 
2 i-l 

(y i - y i-1> 

( 15 ) 


(t - - V,). 


Subroutine LSQPLT 

The function of subroutine LSQPLT is to draw a least-squares polynomial 
or least-squares cubic spline curve through an input set of x and y values. 

For each set of input data, the input points are plotted with the specified 
symbol and the least-squares curve plotted with the corresponding dashed 
line. The curve-fitting procedure used in this subroutine is similar to that 
used in subroutine CURPLT. A separate least-squares curve is determined for 
each set of input x and y values as a function of the independent variable t 
and a y(x) curve then obtained by interpolation. If the user specifies that 
a least-squares polynomial curve be drawn through the input data, subroutine 
LSQ is called which determines the coefficients of the specified-order 
polynomial. If the user specifies that a least-squares cubic-spline curve be 
drawn, subroutine CSDS is called which determines the coefficients of the 
third— order piecewise polynomials that constitute the spline curve. 

Following the call to the selected least-squares subroutine, the x and y 
values at each input t value are computed and sum-of-the-squares of the 
differences between the input and least— squares values are computed. If the 
user specifies a nonzero value for the parameter IPRINT, a summary of the 
least— squares curve fit process will be printed as illustrated in table V and 
VI. The final execution step is a call to subroutine DASHLN which draws the 
least— squares compute x and y values with the specified dashed line. The 
following additional arguments and parameters are used in this subroutine: 
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ILSQ 

NPOL 

DF 

COEF 

CX,CY 

XNEW,YNEW 

WT 

IPRINT 

ERRX,ERRY 


if = 0, fit input data with least-squares polynomial 
if - 1 ( fit input data with least-squares cubic spline 
order of polynomial if ILSQ = 0 

standard deviation for least-squares cubic spline if 
ILSQ = 1 

two-dimensional array containing coefficients of piecewise 
third-order polynomials of cubic spline 

arrays containing coefficients of least-squares polynomial 
arrays containing new x and y values after curve fit 
weighting factor used in least-squares polynomial curve fit 
if - 0, do not print summary data 
if = 1, print summary data 

sum-of-squares of the differences between input and 
least-squares x and y data 


Subroutine LSQ 

The function of subroutine LSQ is to determine the coefficients of a 
polynomial that best fits an input set of weighted x and y data* The 
polynomial is in the form 


2 M 

p(x ) = a + ax + a„x. + . • • + ax. 
i o 1 i 2 i Mi 


(16) 


where M is the order of the polynomial. The sum-of-the-squares differences 
between the weighted input and polynomial fit values are 


S 


2 (P(x ) - w y ) 2 
i=l 


(17) 


i 

i 
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The least or minimum error is 


as 

9a. 


N r -i 

. = 2 .2 «*!> - w.y.l 


= 0 


(18) 


where j = 0, 1, 2 • • . , M. This reduces to a set of simultaneous equations 
in the form 


5) x. j + a 2 x. 

" l 1 l 


j+1 


+ a V x. ^ +M = Tw.y.x. 
M " 1 " 11 1 


(19) 


which can be solved for the polynomial coefficients aj using a simplified 
solution technique that takes advantage of the symmetric properties of the 
matrix of terms generated by the left-hand side of equation (19). The 
following arguments and parameters are used in this subroutine: 

X# Y arrays containing input x^ and y^ data 

W array containing input weighting factor w^ 

NP number of input points N 

N order of polynomial M 

C array containing values of the coefficients of the 

least-squares polynomial 
A two-dimensional work array 

Subroutine CSDS 

The function of subroutine CSDS is to fit a least-squares cubic-spline 
through an input set of x and y data. The method used in this subroutine 
defines a continuous cubic-spline function in the form of the third-order 
polynomial 


2 3 

P(x) = a + a„h. + a^h. + a_h. 

o 1 l 2 l 3 l 


( 20 ) 


for i = 1, 2, 3, . . ., N - 1 where h^ = (x - x^) and N is the number of 
input points. The coefficients a Q , a^ , a 2 , and a^ are computed in a least- 
squares manner so that 
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where the smoothing parameter SP is in the interval N - /2N^ SP _< N + and 
6y^ is specified allowable standard deviation in the error in y^. 

A detailed discussion of the least-squares cubic-spline method is pre- 
sented in reference 3 and, therefore, will not be included in this report. 

This subroutine is also a part of the standard math-library subprogram package 
available on the Langley CDC computer system and is identified by the same 
name and argument list. A complete description of the input and output 
parameters in the argument list are presented at the beginning of the listing 
of the subroutine in appendix A. 


Subroutine AXISLB 

The primary function of subroutine AXISLB is to draw scale values on the 
x or y axis. The scale value is drawn in one-inch increments starting at the 
origin and continuing for a specified length. Each scale value drawn is 
centered adjacent to the inch mark for the y-axis or centered below the inch 
mark for the x-axis. The following is a description of the subroutine argu- 
ments; 


IAXIS 


XORG, YORG 

ORG 

SCALE 


if = 1, draw x-axis scale values 
if = 2, draw y-axis scale values 
Xp and yp values of axis origin, in. 
scale value at origin 
increment in scale value per inch 
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HT height of scale values, in. 

NDIG number of significant figures to draw on the right side of 

the decimal point (NDIG = -1 will drop decimal point) 

Subroutine COEFSY 

The function of subroutine COEFSY is to draw a specified standard 
aerodynamic label as illustrated in table II. The first seven labels are the 
standard designations for the lift, drag, pitching-moment, rolling-moment, 
yawing-moment, and side-force coefficients and lift-drag ratio. The last two 
labels are the standard designations for the angle of attack and angle of 
sideslip. The width-height ratio for each label is also listed in table II 
and is useful to the user when determining the Xp and yp values needed to 
position the label. If the user desires to add coding to draw other labels, 
the additional coding can be inserted following the statement 3 CONTINUE. The 
following is a description of the subroutine arguments: 

XO x D location of left-hand edge of main character in label, 

XT 

in. 

YO yp location of centerline of main character in label, in. 

HT height of main character, in. 

ISYM identifying number of label 

Subroutine LAM 

The function of subroutine LAM is to draw the script "jf" character used 
as the subscript for the rolling-moment coefficient label generated by 
subroutine COEFSY. This special character is not available as part of the 
standard character series contained in the CALCOMP plotting subprogram 
package. 


Subroutine GRIDLN 

The function of subroutine GRIDLN is to draw a series of horizontal and 


vertical grid lines within a specified frame size. This process is 
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accomplished using simple straight line plotting pen movements. The following 

is a description of the subroutine arguments: 

x p and y p values for the lower left-hand corner of the 

frame, in. 
width of frame, in. 
height of frame, in. 
number of grid lines per inch 

Subroutine RUNKEY 

The function of subroutine RUNKEY is to draw a legend containing a list 
of run numbers for the data plotted and a corresponding list of symbols. The 
run numbers are spaced so that a corresponding list of descriptors can be 
typed on a sheet of gum— bacic paper using the IBM Executive model typewriter 
and then cut and attached to the legend adjacent to the list of run numbers. 
This particular typewriter is widely used at Langley to type figure titles and 
legends. After drawing the symbols and run numbers the words "Symbol" and 
"Run" are drawn and underlined above the corresponding lists. The user often 
desires to scale a series of plots to either a smaller or larger size using 
the magnification factors available as part of the postprocessor PLOT control 
card. Using these magnification factors would also cause a corresponding 
change in the legend which would prohibit the use of the IBM Executive Model 
typewriter to type the list of descriptors. To overcome this problem, a 
magnification factor, XM, has been included in this subroutine that adjusts 
the letter and number heights and spacings to compensate for the postprocessor 
scaling. The following subroutine arguments and parameters are used in this 
subroutine : 

XO, YO The Xp and y p values of the lower left-hand corner of the 

legend, in. 


XO, YO 

XL 

YH 

NDIV 
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RUN 


LSYM 

ISIZE 

HT 

NRNMAX 

XM 


array containing input run numbers 
array containing input symbol order 
symbol size 

letter and number character height, in. 
maximum allowable number of runs per sheet 
magnification factor for character and spacings 


APPLICATION OF COMPUTER PROGRAM 

The basic program variables have been dimensioned for a maximum of 10 
runs per sheet, 10 plots per sheet, 50 data points per run, 5000 data points 
per test, and 300 wind-tunnel test parameters per data point. The following 
procedure outlines the changes necessary to either increase or decrease these 
basic program variables: 

(1) To change the maximum number of allowable runs per sheet, change the 
dimension of the variable RUN and the value of the parameter NRMAX in the pro- 
gram PLOT. 

(2) To change the maximum number of plots per sheet, change the dimension 
of the variables XPLT, YPLT, YTAPE, YOFFSET, YSCALE , YLABEL, XT APE, XOFFSET , 
XSCALE, XLABEL, VAR, and NLOC and the value of the parameter NPLMAX in program 
PLOT. 

(3) To change the maximum number of points per run, change the dimension 
of the variables VAR, ARRAY 1 , ARRAY 2 and the value of the parameter NDMAX in 
program PLOT? change the dimension of the variables TIN, TEDT, TSAVE, YSAVE, 
IEDT, ITEMP, and ISAVE in subroutine EDIT? change the dimension of the vari- 
ables TIN, ISAVE, TSAVE, IEDT, TEDT, IUN, and TUN in subroutine PRTEDT? change 
the dimension of the variables MX, MY, and A in subroutines CURPLT and DASHLN? 
and change the dimension of the variables XNEW, YNEW, COEF, and SDV and the 
value of the parameter NMAX in subroutine LSQPLT. 
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(4) To change the maximum number of allowable data points per test, 
change the dimension of the variable NPT and the value of the parameter NPMAX 
in program PLOT. 

(5) To change the maximum number of allowable wind-tunnel test parameters 
per data point, change the dimension of the variable DATA and the value of the 
parameter IMAX in subroutine DATADK. 

The best way to demonstrate the application of the program is to present 
and discuss the output for several sample cases. Two sample cases will be 
presented (1) a single plot with three variables per sheet and (2) four plots 
with a single variable per sheet. The input data for these two cases are 
presented in table III. Both cases contain one sheet namelist with identical 
sheet and run numbers. The tabulated output for the first case is presented 
in table IV and the plotted output, in figure 1. If the user desires to plot 
dashed lines instead of symbols with connecting solid lines, the call to sub- 
routine CURPLT in program PLOT can be replaced with a call to subroutine 
DASHLN which will produce the plot presented in figure 2. In a like manner, 
if the user desires to least-squares curve fit the input data, the call to 
CURPLT can be replaced with a call to subroutine LSQPLT which will produce the 
plot presented in figure 3(a) for the least-squares polynomial option and in 
figure 3(b) for the least-squares cubic-spline option. The programming 
changes required to obtain these optional plots are very simple and require 
very little programming knowledge. 

Some programming knowledge is required if the user desires to generate a 
report quality plot on blank (no grid) paper by drawing the appropriate grid 
lines, axis scale values and labels, and the legends. The input requirements 
for the subroutines that generate the various parts of the plot are relatively 
simple and, after a short period of time, the average user should become very 
proficient in their use. The computer code required to generate report 
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quality plots for the two sample cases presented in table III are presented in 
table VII and VIII, respectively, and the corresponding plots in figures 4 and 
5. Each set of code was inserted in the PLOT program following card number 
233 (see appendix A) and the original code from cards number 237 to 269 
deleted. Careful examination of the code for both cases will reveal that the 
same general procedure was used to generate a particular plot* This procedure 
is outlined as follows: 

( 1) position plotting pen and draw grid lines (call subroutine GRIDLN) 

(2) draw all zero lines (call subroutine ZEROLN) 

(3) draw scale values on x-axis (call subroutine AXISLB) 

(4) position and draw label for x-axis (call subroutine COEFSY) 

(5) repeat steps 3 and 4 for y-axis 

(6) position and draw legend (call subroutine RUNKEY) 

The two types of plots presented represent only a small example of the 
wide variety of the types of plots that can be generated using the PLOT pro- 
gram as the baseline. Most engineers that test in the same facility generally 
prefer to present the experimental data taken during a typical test program in 
a well-established standard format. The computer code to generate these stan- 
dard report quality plots can be written by a programmer with more experience 
and simply inserted in the FLOT program code by the user with either the EDIT, 
XEDIT, UPDATE, or MODIFY editing postprocessor commands. 

CONCLUDING REMARKS 

The wind-tunnel data plotting program PLOTWD described in this report has 
been used successfully at Langley for several years with very few reported 
execution problems. The required format for the input data is somewhat 
restrictive, but the required user input is very simple and the variety of 
plots generated numerous. With relatively minor code changes, the original 


version of the program can be changed to handle either larger or smaller 



amounts of test data. The program code presented in this report was written 
especially for use on the Langley CDC computer system. The program can be 
adapted to other computer systems provided the system selected has a 
random-access file capability and a similar CALCOMP software package. A copy 
of the basic source code for this program can be obtained using the following 
control statements : 

GET, PLOTWD/UN = 621288N 

UPDATE, F,P = PLOTWD, S = SOURCE, I = 0, L = 0. 
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APPENDIX A 


COMPUTER LISTING OF WIND-TUNNEL DATA PLOTTING PROGRAM PLOTWD 
This appendix contains a computer listing of the wind-tunnel data plott- 
ing program PLOTWD which consists of a main program, eighteen subroutines, and 
two function subprograms. 
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LISTING OF DECK J PLOT 


PAGE 


1 


CARD NO. 

1 
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1TA PEI 1-1001, TAP FI 2-1001) 

PROGRAM TO PLOT AND FAIR WIND TUNNEL DATA 

CODED BY — HARRY L. MORGAN NASA/L ARC/TAO/AAB 


********************************** **** *********** ******************** 
♦ 

* 

* 

* CARD NO. 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

♦ 

♦ 

* 

♦ 

* 

' * 
* 
* 
* 
* 
♦ 
* 
* 
* 
* 
* 


— DEFINITION OF SETUP DECK CARDS — 
DESCRIPTION 
TEST - TEST NUMBER 

SHEETV - WIDTH OF PLOTTING SHEET (INCHES) 
SHEETH - HEIGHT OF PLOTTING SHEET (INCHES) 
SPACE - SPACE BETWEEN PLOTTING SHEETS (INCHES) 


ISYM - STARTING SYMBOL NUMBER 
ISIZE - SYMBOL SIZE (1 ■ SMALL 


2 ■ MEDIUM 


IOP - PLOTTING OPTION I0P-0 PLOT DATA ONLY 

IOP-l PLOT AND FAIR DATA 

TENSION - SPLINE TENSION FACTOR (FOR IOP-1 ONLY) 

I EDIT - EDIT OPTION I EDTT»0 NO EDITING 

IEDIT.l €OIT WITH GIVEN TOL 
TOLR - EDIT TOLERANCE FOR INDEPENDENT VARIABLE 


NPLOT - NUMBER OF PLOTS PEP SHEET 

YTAPE r ARRAY LOCATION OF Y VARIABLE (2(3F 

YOFFSET - ORIGIN OFFSET OF Y VARIABLE (INCHES) 

YSCALF - SCALE VALUE PER INCH FOR Y VARIABLE 

YLA9EL - LABEL FOR Y VARIABLE 

XT APE - ARRAY LOCATION OF X VARIABLE 

XQFFSET - ORIGIN OFFSET OF X VARIABLE (INCHES) 


■1001, 


PT 

1 



PT 

2 



PT 

3 



PT 

4 



PT 

5 

1983 


PT 

6 



PT 

7 

******** 

** 

PT 
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♦ 

PT 
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* 

PT 

10 


* 

PT 

11 

format 

* 

PT 

12 


* 

PT 

13 

(F10.0) 

. * 

PT 

14 


' * 

PT 

15 

(3*10.0) 

* 

PT 

16 


* 

PT 

17 


* 

PT 

18 


* 

PT 

19 

(2F10.0) 

* 

PT 

20 

LARGE) 

* 

PT 

21 


* 

PT 

22 

(2F10.0) 

* 

PT 

23 


* 

PT 

24 


* 

PT 

25 


♦ 

PT 

26 

(2F10.0) 

* 

PT 

27 

ERANCE 

* 

PT 

28 


* 

PT 

29 


# 

PT 

30 

(*10. 0) 

* 

PT 

31 


* 

PT 

32 

(F10.0) 

* 

PT 

33 


* 

PT 

34 

. 0* A10 ) ) 

* 

PT 

35 


* 

PT 

36 


* 

PT 

37 


* 

PT 

38 


* 

PT 

39 


* 

PT 

40 
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CARO NO • 


41 

C 

* 

XSCALE - SCALE VALUE PER INCH FOR X VARIARLF 

* pt 

41 


c 

* 

XLABEL - LABEL ROR X VARIARLE 

* PT 

42 


c 

* 

NOTE - CARO 8 IS RFPEATED NPLOT TIMES 

* PT 

43 


c 

* 


* PT 

44 

45 

c 

* 

— DEFINITION OF SHEET NAMELIST PARAMETERS — 

* PT 

45 


c 

* 


★ PT 

46 


c 

♦ 

iSHEET - NAMELIST LABEL 

* PT 

47 


c 

* 

NO - SHEET NUMBER 

* PT 

48 


c 

* 

RUN - PUN NUMBERS TO BF PLOTTED ON SHEET (MAXIMUM 

OF 10) ’ 

* PT 

49 

50 

c 

* 

NEWCASE - IF NFV SETUP DECK FOLLOWS SET NEWCASF-1 

* PT 

50 


c 

♦ 

$END - FND OF NAMELIST 

* PT 

51 


c 

* 


3 

► PT 

52 


c 

*********************** ********************************************** py 

53 


c 




PT 

54 

55 



INTEGER RUN 


PT 

55 


c 




PT 

56 


c 


DIMENSION ARRAY LSYM WHICH CONTAINS ORDER OF PLOTTING 

SYMBOLS 

PT 

57 


c 




PT 

58 




DIMENSION LSYM ( 22 ) 


PT 

59 

60 

c 




PT 

60 


c 


DIMENSION RUNCNRNMAX) WHERE NRNMAX IS THE MAXIMUM NUMBER OF 

PT 

61 


c 


ALLOWABLE RUNS PER SHEET 


PT 

62 


c 




PT 

63 




DIMENSION RUN ( 10 ) 


PT 

64 

65 

c 




PT 

65 


c 


DIMENSION NPT(NPMAX) WHERE NPMAX IS THE MAXIMUM NUMBER 

OF DATA 

PT 

66 


c 


POINTS FOR TOTAL TF ST 


PT 

67 


c 


NOTE — 00 NOT DESTROY THIS ARRAY DURING PROGRAM EXECUTION 

PT 

68 


c 




PT 

69 

70 



DIMENSION NPT ( 5000) 


PT 

70 


c 




PT 

71 


c 


DIMENSION EACH OF THE FOLLOWING ARRAYS BY NPLMAX WHICH 

is the 

PT 

72 


c 


MAXIMUM NUMBER OF ALLOWABLE PLOTS PER SHEET 


PT 

73 


c 




PT 

74 

75 



DIMENSION XPLT(IO), YPLT(lo), YTAPEC10)# YOFFSETUO), YSCALE(IO), 

PT 

75 




lYLABfL(lO), XTAPE(IO), XOFFSET { 10) * XSCALE(IO), XLABEL ( 10) 

PT 

76 


c 




PT 

77 


c 


DIMENSION EACH OF THE FOLLOWING ARRAYS BY NDMAX WHICH 

IS THE 

PT 

7B 


c 


MAXIMUM NUMBER OF ALLOWABLE POINTS FOR A GIVEN RUN 


PT 

79 

80 

c 




PT 

80 
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CARD NO, 
81 

85 

90 

95 

100 

105 

110 

115 

120 . 


C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DIMENSION X ( 50) , Y(50), T ( 50 > 

DIMENSION VAR(NDMAX,NVAR) AND NLOC(NVAR) WHERE N VAR-2+NPLM AX+l 

DIMENSION VAR { 50# 21)# NL0C(21) 

DIMENSION ARRAY1(7*NDMAX,GE,500) AND ARRA Y2 ( 3*NDMAX ) 

COMMON /WORK/ ARRAYK500) 

COMMON /PLT / ARRAY? (150) 

SHEET NAMELIST PARAMETERS 

NAMELIST /SHEET/ NO, RIJN, NEWC ASE 

DEFINE INPUT TAPE NUMBER 

DATA J READ/5/ 

DEFINE OROER OF PLOTTING SYMBOLS 

DATA LSY M/1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,17,18,19,20,21,22 
1 / 

DEFINE MAXIMUM NUMBER OF ALLOWABLE RUNS PER SHEET 
DATA NRNMAX/10/ 

DEFINE MAXIMUM NUMBER OF ALLOWABLE PLOTS PER SHEET 
DATA NPLMAX/10/ 

DEFINE MAXIMUM NUMBER OF ALLOWABLF DATA POINTS EOR A GIVFN RUN 
DATA NDMAX/50/ 

DEFINE MAXIMUM NUMBER OF ALLOWABLF TOTAL DATA POINTS 
NP MA X *5000 


PT 

81 

PT 

82 

PT 

83 

PT 

84 

PT 

85 

PT 

86 

PT 

87 

PT 

88 

PT 

89 

PT 

90 

PT 

91 

PT 

92 

PT 

93 

PT 

94 

PT 

95 

PT 

96 

PT 

97 

PT 

9e 

PT 

99 

PT 

100 

PT 

101 

PT 

102 

PT 

103 

PT 

104 

PT 

105 

PT 

106 

PT 

107 

PT 

108 

PT 

109 

PT 

110 

PT 

111 

PT 

112 

PT 

112 

PT 

114 

PT 

115 

PT 

116 

PT 

117 

PT 

118 

PT 

119 

PT 

120 
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4 


CARD NO. 


121 C DEFINE PROGRAM OPTIONS PT 121 

C PT 122 

C IF IPRT-O* PRINT EDITED DATA PT 123 

C IF IPRT-1* DO NOT PRINT EDITED DATA PT 124 

125 IPRT-0 PT 125 

C LZERO-O* DRAW ZERO LINES PT 126 

C LZERO-1* DO NOT DRAW ZERO LINES PT 127 

LZERO-1 PT 128 

C NORG-O* DRAW ARROWHEAD AT ORIGIN OF PLOTTING VARIABLES PT 129 

130. C NORG-1* DO NOT DRAW ARROWHEAD AT ORIGIN OF PLOTTING VARIABLES PT 130 

NORG-O PT 131 

C PT 132 

C . SET PROGRAM CONSTANTS PT 133 

c PT 134 

135. IDISK-0 PT 135 

0 PT 136 

C INITIALIZE PLOTTING DEVICE PT 137 

C PT 138 

CALL PSEUDO PT 139 

140. CALL LEROY PT 140 

C PT 141 

C READ SETUP DECK PT 142 

c PT 143 

1 CALL SETUP (ITEST*SHEETW*SHEFTH, SPACE* ISYM*ISI ZE*IOP*TENSION, IEDIT PT 144 

145. 1* TCLR*IW*NPLOT* YTAPE* YOFFSET* YSCALE* YLABEL* XTAPE* XOFFSET* XSCALE* XL PT 145 

2ABEL*JREAD*NPLMAX*IERR) PT 146 

IF (IERR.NE.O) GO TO 23 PT 147 

c PT 148 

REWIND 7 PT 149 

150 REWIND 11 PT 150 

KEDT-0 PT 151 

C PT 152 

C LOAD DATA DISK PT 153 

C PT 154 

155 IF ( I DISK ) 3*2*3 PT 155 

2 CALL DATADK (0* ITEST* IRUN, VAR* NVAR* NLOC* NP* NPT* NPMAX* NOMAX* IERR) PT 156 

IF (IFRR.NE.O) GO TO 23 PT 157 

C PT 158 

3 IDISK-1 PT 159 

160 PRINT 24 PT 160 
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161 

165 

170 

175 

180 

185 

190 

195 

200 


C 

C 

C 


A 


C 

C 

C 


5 

6 

7 

8 

9 
C 
C 
C 

10 
11 

C 


12 

C 

C 

C 


LOAD PLOTTING VARIABLES 

DO 6 1*1# N PLOT 

J2-2+I 

J1-J2-1 

NLOC( Jl)-IFIX(YTAt»c (T ) ♦0.0001) 

NLTC( J2)-IFIX( YTAPECD+O.OOOl) 

NVAR-2*NPL0T+1 
NLOC (NVAR )-lW 

CHECK FOR DUPLICATE PLOTTING VARIABLES 

DO 9 I«l> NPLOT 
YPLT( I )-XPLT(I )-0. 

IF (I.EG.l) GO TO 9 

Jl-I-1 

DO 8 J-1,J1 

IF (XTAPE CJ )-XTAPE(I ) ) 6.5#6 

IF (XOFFSET(J).EO.XOFFSET(I ) . AND . XSC A LE ( J J.EQ.XSCALE(I)) 

1.0 

IF {YTAPE(J J-YTAPE(I) ) 8,7, 8 

IF (YOFFSETt J) .EO.YOFFSETm .AND.YSCALEt J).EO.YSCALF(I) ) 

1.0 

CONTINUE 

CONTINUE 

READ SHEET NUMBER AND PUNS TO BE PLOTTED 

DO 11 I«l# NRNMAX 
RUN ( I )»0 
NEWCASE-NO-O 

READ ( JREAD* SHEFT ) 

IF ( EOF ( JRE AD) ) 23,12 
CONTINUE 

PRINT 25, NO, PUN 
PERFORM GRID CHFCK 


PT 161 
PT 16? 
PT 163 
PT 166 
PT 165 
PT 166 
PT 167 
PT 168 
PT 169 
PT 170 
PT 171 
PT 17? 
PT 173 
PT 176 
PT 175 
PT 176 
PT 177 
PT 178 
PT 179 
XPLT(I)-1 PT 180 
PT lftl 
PT 182 
YPL T( I ) *1 PT 183 
PT 186 
PT 185 
PT 186 
PT 187 
PT 188 
PT 189 
PT 190 
PT 191 
PT 192 
PT 193 
PT 196 
PT 195 
PT 196 
PT 197 
PT 198 
PT 199 
PT 200 
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CAPO NG * 
201 

205 

210 

215 

220 

225’ 

230. 

235' 

240 


c PT 201 

CALL CALPLT (0.»0.#-3> PT 202 

CALL GPIDCK p T 203 

c PT 204 

C POSITION PLOTTING D E N AT LOWER LEFT HAND CORNER OF PLOTTING SHEFT PT 205 

c PT 206 

CALL CALPLT (2*#2.#-3) PT 207 

c PT 208 

C LABEL TEST# SHEFT# AND RUN NUMBERS p T 209 

c PT 210 

SZ-0.15 PT 211 

D— 1.9 PT 212 

C LABEL TEST NUMBFR PT 213 

CALL NOTATE (D#-2 .0, SZ#5HTEST #0.#5) PT 214 

D-30.+SZ/7.+D PT 215 

TP«FLOAT ( ITEST) PT 216 

CALL NUMBER <!>#-? .0# SZ #TP#0. #-l ) PT 217 

0-36,*SZ/7.+D PT 218 

C LABEL SHEET NUMBER PT 219 

CALL NOTATE (n,-2.0#SZ#7H,SHEET #0.#7) PT 220 

D«42 • *SZ /7#+D p T 221 

TP-FLOAT(NO) PT 222 

CALL NUMBER (D#-2.0# SZ#TP#0.#-1> PT 223 

D-36.*SZ/7.+D PT 224 

C LABEL RUN NUMBERS PT 225 

CALL NOTATE ( 0#-? .0# SZ # 6H, RUNS #0.#6) ■ PT 226 

D-42.*SZ/7.+D p T 227 

DO 13 I-1#NRNMAX PT 228 

T P«FL CAT ( RUN(I) ) PT 229 

CALL NUMBER ( D#-2 *0# SZ#TP,0. #-l ) p T 230 

13 D-42 . + SZ/7. + D PT 231 

IX*IF IX ( D) PT 232 

DMAX-FL0AT(IX)+SPACE+2.0 PT 233 

c PT 234 

C DRAW PLOTTING SHEET FRAME PT 235 

c PT 236 

CALL CALPLT (0.#0.#3) PT 237 

CALL CALPLT < SHEETV# 0. # 2) PT 238 

CALL CALPLT (SHEETW# SHEETH#? ) PT 239 

CALL CALPLT (0.#SHEETH#2) PT 240 
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CARD NO. 


241 

245 

250 

255 ! 

260 

265. 

270 

275 

280 


CALL CALPLT (0.;0.»2) 

C 

C LABEL ORIGIN* VARIABLE DESCRIPTION* AND SCALE 
C 

DO 15 I»1*NPL0T 

IF (YPLTCI).EQ.l*) GO TO 14 

D«YOPFSET<n 

IF (NORG.EQ.O) CALL PARROV (-0. 3* D* 0 . * D* -2 . 0* . 3 ) 
CALL NOTATE (-1. 8* D* 0. 15* YL AB EL (I )* 0.* 10) 
D*Y0FF$ET(I)-0.2 

CALL NUMBER (-1. 8* D* 0 . 15* YSC ALE ( T ), 0 . * 4 ) 

14 IF {X PLT (I)« 60.1 .) GO TO 15 
D«X0FFSET{ I) 

IF (NORG.EQ.O) CALL PARROV (D#~0.3*D*0.*~2.0*.3) 
D*X0FFSET(I)-0.65 

CALL NOTATE ( D* -0 .6* 0 . 1 5* XL ABEL ( I ) * 0. *10 ) 

CALL NUMBER (0*-0.8*0.15* XSCALE ( T )* 0. * 4) 

15 CONTINUE 
C 

C DRAW ZFRO LINES 

C 

IF (LZERO.NE.O) GO TO 1? 

DO 17 I-l*NPLnT 

IF ( YPLT (I).EO.l.) GO TO 16 

CALL ZEROLN (0,*YOFFSET(I). SHEETV*1 ) 

16 IF (XPLT(I).EQ.l. ) GO TO 17 

CALL ZEROLN ( XOFFSET (I )*0.*SHFETH* 2) 

17 CONTINUE 

18 CONTINUE 
C 

C READ DATA FROM DATA DISK AND PLOT 

C 

KS YM* ISYM-1 
C 

DO 21 IR-1*NRNMAX 
C DEFINE RUN NUMRFR 

IRUN-PUN(IR) 

C DEFINE PLOTTING SYMBOL 

KSYM-KSYM+1 

IF (KSYM.GT .22) KSYM-ISYM 


PT 241 
PT 242 
PT 243 
PT 244 
PT 245 
PT 246 
PT 247 
PT 248 
PT 249 
PT 250 
PT 251 
PT 252 
PT 253 
PT 254 
PT 255 
PT 256 
PT 257 
PT 258 
PT 259 
PT 260 
PT 261 
PT 262 
PT 263 
PT 264 
PT 265 
PT 266 
PT 267 
PT 268 
PT 269 
PT 270 
PT 271 
PT 272 
PT 273 
PT 274 
PT 275 
PT 276 
PT 277 
PT 278 
PT 279 
PT 280 
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CARO NO, 


281 


IF (KSYM.LT.U KSYM-1 


PT 

281 



NSYN-LSYM(KSYH) 


PT 

282 


C 



PT 

283 



IF (IRUN.EQ.O) GO TO 21 


PT 

284 

285 

c 

READ DATA FROM DATA DISK 


PT 

2 85 



CALL DATADK ( 1 # I TEST# I RUN# VAR # NVAR# NLOC# NP#N PT# 

NPMAX,NDMAX#TERR) 

PT 

286 

, 


IF (IEPR.NE.O) GO TO 21 


PT 

287 


c 

LOAD AND PLOT DATA 


PT 

288 



DO 20 I«1,NPL0T 


PT 

289 

290. 


J2-2+I 


PT 

2 90 



J1-J2-1 


PT 

291 



IX«NL0C(J2) 


PT 

292 


c 

LOAD PLOTTING ARRAYS 


PT 

293 



N*0 


PT 

294 

295 


DO 19 J«1#NP 


PT 

295 



XP«VAR(J,J2)/XSCALE< IJ+XOFESETCI) 


PT 

296 



YP«VAR(J#Jl)/Y$CALF(I)+YOFFSETm 


PT 

297 



TP-VAR ( J#NVAR> 


PT 

298 


c 

CHECK TO SFF IF DATA IS WITHIN PLOTTING 

FRAME 

PT 

299 

300 


IF (XP.LT.O.O.OR.XP.GT.SHFETV) GO TO 19 


PT 

300 



IF (YP. LT.O.O.nR.YP.GT.SHEEYH) GO TO 19 


PT 

301 



N-N+l 


PT 

302 



X(N)*XP 


PT 

303 



Y(N) ■ YP 


PT 

304 

305 


T(N)*TP 


PT 

305 


19 

CONTINUE 


PT 

306 


C 

PRINT MESSAGE IF SOME POINTS NOT WITHIN 

PLOTTING FRAME 

PT 

307 



JS-NP-N 


PT 

308 



IF (N.NE.NP) PRINT 26# IRUN#JS#I 


PT 

309 

310. 

C 

EDIT DATA 


PT 

310 



I ERR "0 


PT 

311 



IF (I EDI T. N E. 0 ) CALL EDIT (T, X# Y# N» NS YM# ISI ZE# IRUN, TQLP# KEDT, I ERR ) 

PT 

312 



IF (IERR.NE.O) GO TO 20 


PT 

313 


C 

PLOT VARIABLES 


PT 

314 

315 


JS-0 


PT 

315 



IF (IX.NE.IW) JS-1 


PT 

316 



CALL CURPLT <T# X# Y#N # J S# NSYM# ISIZE# TOP# IRUN# TENS ION ) 

PT 

317 


c 

DO NEXT VARIABLE 


PT 

318 


20 

CONTINUE 


PT 

319 

320 

C 

DO NEXT RUN 


PT 

320 
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CARD NO. 


321 

325 

330 

335 

340 

345 

350 

355 


21 

C 


C 

C 

C 


C 

C 

C 

22 


23 

C 

C 

C 


C 

c 

c 


c 

24 

25 

26 

27 

28 


CONTINUE 

ADVANCE TO ORIGIN OF NEXT FRAME 
D-SHEETW+SPACE+2.0 
IF (D.LT.DMaX) D-DMAX 
CALL NFRAME (D#0.) 

READ NEXT SHEET NAMELIST UNLESS NEW 

IF (NEWCASE.NF.O) GO TO 2? 

GO TO 10 

PRINT EDITED DATA 


PT 321 
PT 322 
PT 323 
PT 324 
PT 325 
PT 326 

SETUP DECK FOLLOWS PT 327 

PT 328 
PT 329 
PT 330 
PT 331 
PT 332 
PT 333 
PT 334 
PT 335 
PT 336 
PT 337 
PT 338 
PT 339 
PT 340 
PT 341 
PT 342 
PT 343 
PT 344 
PT 345 


IF CIEDIT.NE.O.AND.IORT.EO.O) CALL PRTEDT 
PRINT 27 
GO TO 1 

IF (I EDIT. NE.O. AND. IPRT.EQ.O) CALL PRTEDT 
FINALIZE PLOTTING DEVICE 

CALL NOTATE (2,# 2*#0.5#15HEND OF PLOTTING# 90 •# 15 ) 
CALL NFRAME (4.#0.) 

CALL CALPLT (0.#0.#99Q) 

PRINT CLOSING MESSAGE 


PRINT 

STOP 


28 


FORMAT (1H1#3X#12HDATA PLOTTED) 

FORMAT { /4X# 10HSHEFT NO ■ # 15# 5X# 6HRUN $ . 
1 # ) # 17 ) ) 

FORMAT ( /4X#9HF0R RUN «#I7»3H # #I4#55H 
1TTING SHEET FRAMF FOR PLOT «#I3) 

FORMAT ( /4X# 22HNE W CASE SHOULD FOLLOW) 
FORMAT ( /4X# 22HALL PLOTTING' COMPLETED ) 
END 


#9(I7#1H#)#I7/(30X#9CI7#1H 
POINTS WERE NOT WITHIN PLO 


PT 

PT 

PT 

PT 

PT 

PT 

PT 

PT 

PT 

PT 

PT 

PT 


346 

347 

348 

349 

350 

351 

352 

353 

354 
3 55 
356 
357- 
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1 


5 


10 


15 


20 


25 


30 


35 


40 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


SUBROUTINE SETUP (TTEST* SHEETV, SHEETH, SPACE,ISYM, TSIZE* IPP, TENSION ST 1 
I ED IT* TOLR, IW*NPLOT*YTAPE*YOFFSFT*YSCALE*YLAB£L*XTAPE*XOFF$ET,XSC ST 2 


2ALE* X LABEL# J READ* NPLMAX*IERR> ST 3 

ST 4 

ROUTINE TO READ PLOTTING SETUP DFCK ST 5 

ST 6 

CODED BY — HARRY L. MORGAN NASAVL ARC/TAD/AAB 1983 ST 7 

ST 8 

PARAMETER DEFINITION ST 9 

ITEST - TEST NUMBER ST 10 

SHFETW - WIDTH OF PLOTTING SHEET (IN) ST 11 

SHEETH - HEIGHT OF PLOTTING SHEET (IN) ST 12 

SPACE - SPACE 8FTWFEN PLOTTING SHEETS (IN) ST 13 

ISYH - STARTING SYMBOL NUMBER ST 14 

IOP - PLOTTING OPTION - IDP-O* PLOT SYMBOLS ONLY ST 15 

IOP-1* PLOT AND FAIR DATA ST 16 

TENSION - SPLINE TENSION FACTOR ST 17 

IEDIT - EDIT OPTION - IEOTT-0* 00 NOT SORT AND EDIT DATA ST 18 

IEOIT-1* SORT AND EDIT DATA ST 19 

TOLR - EDIT TOLERANCE OF INDEPENDENT VARIABLE ST 20 

IW - DATA ARRAY LOCATION OF INDEPENDENT VARIABLE ST 21 

NPLOT - NUMBER OF PLOTS PPR SHEFT ST 22 

YT APE - DATA ARRAY LOCATION OF Y VARIABLE ST 23 

YOFFSET - ORIGIN OF Y VARIABLE (IN) ST 24 

YSCALE - VALUE PER INCH FOR Y VARIABLE ST 25 

YL ABEL - LABEL FOR Y VARIABLE ST 26 

XT APE - OATA ARRAY LOCATION OF X VARIABLE ST 27 

XOFFSET - ORIGIN OF X VARIABLE (IN) ST 28 

XSCALE - VALUE PER INCH FOR X VARIABLE ST 29 

XL ABEL - LABEL FOR X VARIABLE ST 30 

JREAD - TAPE NUMB FR OF INPUT FILE ST 31 

IE RR - ERROR INDICATOR - IERR-O* NO ERRORS OCCURRED ST 32 

IERR-1* ERRORS OCCURRED ST 33 

ST 34 

XT APE* XOFFSET* XSC A L E* XL ABEL* YTAPE* YOFFSET* YSCALE* AND YLABEL ST 35 

MUST BE PROPERLY DIMENSIONED IN CALLING PROGRAM ST 36 

DIMENSION XTAPE(l)* XOFFSET(l)* XSCALE(l)* XLABEL(l)* YTAPE(l), YO ST 37 

IFFSETd)* YSCALE ( 1 )* YLABEL(l) ST 38 

DIMENSION AND DEFINE SYMBOL SIZE ARRAY AND NMAX ST 39 

DIMENSION IS (3 ) ST 40 
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CARD NO. 


41 


DATA IS/6HSMALL #6HME0IMM#6HLARGE / 

ST 

41 


C 


ST 

42 


C 

INITIALIZE ERROR INDICATOR 

ST 

43 


C 


ST 

44 

45 


IERP-0 

ST 

45 


C 


ST 

46 


C 

READ AND PRINT TEST NUMBER 

ST 

47 


C 


ST 

48 



READ ( JREAD#6) D1 

ST 

49 

50 


IF (EOF(JREAD)) 4# 1 

ST 

50 


1 

ITEST-IFIX(D1+0.0001) 

ST 

51 



PRINT 7 # ITEST 

ST 

52 


C 


ST 

53 


C 

READ AND PRINT SHEET V* SHFETH# AND SPACE 

ST 

54 

55 

c 


ST 

55 



READ ( JR£ AD# 6) SHEETV# SHEETH# SPACE 

ST 

56 



IF (SHEFTW.LT. 11. 0) $HEFTW*11 .0 

ST 

57 



PRINT 8# $HEETW#SHEFTH, SPACE 

ST 

58 


c 


ST 

59 

60 

c 

READ AND PRINT TSYM AND ISIZE 

ST 

60 


c 


ST 

61 



READ ( JRE AD# 6 ) D1#D2 

ST 

62 



ISYM-IFIX(D1+C.0001) 

ST 

63 



ISIZE-IFIX(D2*0.0001) 

ST 

64 

65 


IF (ISYM.LT.1.0R.ISYM.GT.22) ISYM-1 

ST 

65 



IF (ISIZE. LT.l.rjR. ISIZE. GT. 3) ISIZE-2 

ST 

66 



PRINT 9# I$YM#IS( ISIZE) 

ST 

67 


c 


ST 

68 


c 

READ AND PRINT 10* AND TENSION 

ST 

69 

70 

c 


ST 

70 



READ ( JRE AD# 6) D1#TENSI0N 

ST 

71 



IOP-I FIXCD1+0.0001) 

ST 

72 



IF (IOP.LT.O.OR.IOP.GT.l) I0P-0 

ST 

73 



IF (IOP.EQ.O) PRINT 10 

ST 

74 

75 


IF (IOP.EQ.l) PRINT 11# TENSION 

ST 

75 


c 


ST 

76 


c 

READ AND PRINT IEDIT AND TOLR 

ST 

77 


c 


ST 

78 



READ (JRE AD# 6) D1#T0LR 

ST 

79 

80 


IEDIT-IFIX(D1+0.0001) 

ST 

80 
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i 


LISTING OF DECKS SFTU 0 PAGE 3 

CARD NO. 


01 


IF (IEDIT.LT.O.OPilFDIT.GT.l) IFniT.O 

ST 

01 



IF (IEDIT.EQ.1) PRINT 12, TflL R 

ST 

8? 


C 


ST 

83 


c 

READ AND PRINT IW 

ST 

84 

85 

c 


ST 

85 



READ { JRFAD, 6) 01 

ST 

86 



IW*IFIXtDl+-0*0001 ) 

ST 

87 



PRINT 13, IW 

ST 

88 



IF (IV.LE.O) GO TO 4 

ST 

89 

90 

c 


ST 

90 


c 

READ AND PRINT NPLOT 

ST 

91 


c 


ST 

92 



RF AD ( JRFAD, 6 ) D1 

ST 

93 



NPL0T-IFIX(D1+0.0001) 

ST 

94 

95 


PRINT 14, NPLHT 

ST 

95 



IF (NPLOT. LE.O. OR. NPLOT. GT.NPLMAX) GO TO 3 

ST 

96 


c 


ST 

97 


c 

READ AND PRINT YT APE , YOFF S ET, YSC AL E * YL ABE L, XTAP E, XT F F $ ET» XS CAL E, 

ST 

98 


c 

AND XLAREL 

ST 

99 

100 

c 


ST 

100 



DO 2 I «1, NP LOT 

ST 

101 



RFAD ( JR EAD, 1 5 ) YTAP F (I) , YQ FFSET (I) , YSC AL E ( I ) , YL ABFL ( I ) , YTA PF ( I ) , X 

ST 

102 



10FFSET(I)#XSCALF(I), XLA8EL(I) 

ST 

103 



IX-IF IX ( XT APE ( I) ♦0.0001) 

ST 

104 

105 


I Y«IF IX <YTAPE( D + 0.0001 ) 

ST 

105 



PRINT 16, I # I Y# YOFFS ET ( I ) , V SC ALF ( I ) , YL ABE L ( I ) , I X, XOFF SF T ( I ) , XSC A LF 

ST 

106 



1 ( I ) , X LABEL< I ) 

ST 

107 


2 

CONTINUE 

ST 

108 



GO TO 5 

ST 

109 

11C 

c 


ST 

110 


c 

PRINT ERROR MESSAGES THAT MAY HAVE OCCURRED 

ST 

111 


c 


ST 

112 


3 

PRINT 1 7, NPLOT 

ST 

113 



IERP-1 

ST 

114 

115 


GO TO 5 

ST 

115 


A 

I ERR * 1 

ST 

116 


C 


ST 

117 


C 

RETURN TO CALLING PROGRAM 

ST 

118 


c 


ST 

119 

120 

5 

RETURN 

ST 

120 
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LISTING 
CARO NO 
121 

125 

130 

135 


P AGP 


OF DECK: SETUP 


10 

11 

12 

13 

14 


15 

16 
17 


1 FnRhAT E {4x!24HSTARTING SYMBOL NUMBER ■> 13, 5X# 15HSYMnm SIZ* IS > A6 

‘is:;;; :s «*..« ...... ...» .. 

^ FORMA t* ( ax> ^AHOATA WILL BE FITTED WTTH AN EDIT TOLERANCE -,EB.3> 

S«;S^S^!i 0 Sr5^TrpS D iSl^-:xX/"“SHSLb^ 1 5S..,* i9 HVT 

1A»E,5X>THY0FF<ET,6*»6HY'CAL6»6X»6HYIABEL»').X,5H*TAPF,‘ ; X,7HX0FF?E . 

2X,6HX$CALE,6X,6HXLAPFL> 

FORMAT (3F10.4,A10,3F10.4>A10) ... 

FORMAT (8X,I3»1Y,2(8'X,I4,5X,F7,2,2V,f1E.4,2X>A10)) 

FORMAT </4X,40HFRR0R OCCUPREO — IN SETUP DECK NPLHT ,i*j 

end 


ST 121 
ST 122 
ST 123 
ST 124 
ST 125 
ST 126 
ST 127 
ST 128 
ST 1?Q 
ST 130 
ST 131 
ST 132 
ST 133 
ST 134 
ST 135 
ST 136 
ST 137 
ST 138 
ST 13R- 
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LISTING OF DECKt DATADK PAGE 1 

CARO NO, 


1 


SUBROUTINE DATADK (I CODE* T TEST, IPUN# VAR# N VAR* NLQC#N P#N PT# NPMAX# NDM 

DK 

1 



1AX # I E RR ) 

DK 

2 


C 


DK 

3 


C 

ROUTINE TO READ FROM OR WRITE TO A DATA DISK 

DK 

4 

5 

c 


DK 

5 


c 

CODED BY — HARRY L. MORGAN NAS A/LARC/ TAD/ AAB 1983 

DK 

6 


c 


DK 

7 


c 

PARAMETER DEFINITION 

DK 

8 


c 

► 

DK 

9 

10 

c 

IC ODE - READ OR WRITE OPTION 

DK 

10 


c 

ITEST - TEST NUMBER 

DK 

11 


c 

ICOOE-O# WRITE IN DATA DISK 

DK 

12 


c 

ICODF-1# READ DATA FROM DATA DISK 

DK 

13 


c 

IRUN - RUN NUMRER TO BE READ FROM DATA DISK 

DK 

14 

15 

c 

VAR(I#J) - TWO DIMENSIONAL ARRAY CONTAINING DATA REAO FROM 

DK 

15 


c 

DATA DISK, ( 1*1# , • • • # NP AND J«l# .. ,,#NVAR ) 

OK 

16 


c 

NVAR - NUMBER OF VARIABLES READ FROM DISK FOR FACH POINT OF 

DK 

17 


c 

DATA 

OK 

18 


c 

NLOC(J) - ARRAY CONTAINING DATA ARRAY LOCATIONS OF DESIRED 

DK 

19 

20 

c 

VARIABLES. ( J*l# • . . NVAR ) 

DK 

20 


c 

NP - NUMBER OF DATA POINTS READ FROM DISK FOR DESIRED RUN 

DK 

21 


c 

NPT - ARRAY CONTAINING COUNT NUMBER OF DATA POINTS STORED ON 

DK 

22 


c 

THE DATA DISK 

DK 

23 


c 

NPMAX - MAXIMUM NUMB FR OF DATA POINTS 

DK 

24 

25 

c 

NDMAX - MAXIMUM NUMBER OF ALLOWABLE DATA POINTS FOR ANY GIVEN 

DK 

25 


c 

RUN' 

DK 

26 


c 

IERR - ERROR INDICATOR 

DK 

27 


c 

IERR*0# NO ERRORS OCCURRED 

DK 

28 


c 

IERR-1# ERROR OCCURRED DUPING READ OR WRITE 

DK 

29 

30 

c 

IERP-2# IRUN NOT ON DATA DISK 

DK 

30 


c . 


DK 

31 


c 

NOTE — VAR# NLOC# AND NPT MUST BE PROPERLY DIMENSIONED IN THE 

DK 

32 


c 

CALLING PPOGPAM 

DK 

33 


c 


DK 

34 

35 


DIMENSION VAR ( NDMAX# 1 ) # NLOC(l)# NPT(l) 

DK 

35 


c 


DK 

36 


c 

WORK ARRAY 

DK 

37 


c 


DK 

38 



COMMON /PLT/ NPRT(IO) 

DK 

39 

40 

c 


DK 

40 


51 


LISTING 
CARD NO 

41 

45 

50 

55 

60 

65 

70 

75, 

80 ! 


Kt 

DATADK 

PAGE 

2 

C 

DIMENSION DATA ( IM AX + 3 ) WHERE TNAX IS THE MAXIMUM NUMBER OF 

DATA DK 

41 

C 

VARIABLES STORED ON THE INPUT DATA TAPE 

DK 

42 

c 


DK 

43 


COMMON /WORK/ OATA (303 ) 

DK 

44 

c 


OK 

45 

c 

DEFINE ITAPE (INPUT DATA TAPE NUMBER) AND NDISK (DATA DISK 

DK 

46 

c 

NUMBFR) 

DK 

47 

c 


DK 

48 


DATA ITAPE/1/>NDISK/?/ 

DK 

49 

c 

DEFINE IMAX 

DK 

50 


DATA IMAX/300/ 

DK 

51 

c 


DK 

52 

c 

COMPUTE ROUTINF CONSTANTS 

DK 

53 

c 


DK 

54 


IERP-0 

DK 

55 


13 *IM AX+3 

DK 

56 


REWIND 12 

DK 

57 

c 


DK 

58 

c 

SELECT READ OR WRITE PROCEDURE 

DK 

59 

c 


DK 

60 


IF (ICODE) 1# 1# 16 

DK 

61 

c 


DK 

62 

c 

WRITE DATA ON DATA DISK 

DK 

63 

c 


DK 

64 

c 

INITIALIZE NPT 

DK 

65 

1 

DO 2 1*1, NPHAX 

DK 

66 

2 

NPT( I )*0 

DK 

67 

c 

OPEN OATA DISK 

DK 

68 


CALL OPENMS (NDISK,NPT,NPMAX,0) 

DK 

69 

c 

COMPUTE NEEDED CONSTANTS 

DK 

70 


11 *1 MAX+1 

DK 

71 


K*0 

DK 

72 


IRUNP-0 

DK 

73 


IC-0 

DK 

74 

c 

REWIND INPUT TAPE 

DK 

75 


REWIND ITAPE 

DK 

76 

c 

READ DATA FROM INPUT TAPE 

DK 

77 

3 

READ (ITAPE) IRUN, ITESTC, IPT, (DATA ( I ), I-l, IMAX ) 

DK 

78 


IF (EOF(ITAPE)) 10,4 

DK 

79 

4 

CONTINUE r 

DK 

80 
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LISTING OF DECKS DATADK 


PAGE 


3 


CARD NO, 
81 

85 

90' 

95 

100 

105 

110 

115 

120 


IF (IRUN.GT.99999) GO TO 10 
IF ( ITESTC. NEATEST) GO TO 15 
C 

c NOTE — POINTS OR RUNS CAN BE VOIDED HERE 

C 

C LOAD DATA ARRAY 

DO 5 1-1* IMAX 
J-Il-I 
J3-J+3 

DATA ( J3)» DATA ( J ) 

5 CONTINUE 

DATA (1 )*FLOAT (I RUN) 

DATA (2)«FL0ATCITE$TC) 

DATA (3)»FL0AT(IPT> 

K«K+1 

C CHECK TO SEE IP POINT COUNT EXCEFDS ALLOWABLE MAXIMUM 

IF (K-NPMAX) 6*14*14 

C WRITE ON TAPE 12 RUN NUMBER AND COUNT NUMBER OF FIRST DATA 

C POINT OF RUN 

6 IF ( I RUN P ) 7*8*7 

7 IF C IRUN-IRUNP ) 8*9*8 

8 IC-IC+l 

WRITE (12) IRUN*K 
IRUNP-IRUN 

C WRITE DATA ON DATA DISK 

9 CALL WRITMS (NDI SK* D ATA* 13* K* 0*0) 

C READ NEXT DATA POINT ON INPUT TAPE 

GO TO 3 

C WRITE ON TAPE 12 RUN NUMBER AND COUNT NUMBER OF LAST DATA 

C POINT ♦ 1 OF RUN 

10 K*K+1 

WRITE (12) IRUNP*K 
END FILE 12 
C 

C PRINT SUMMARY OF RUNS STORED ON DATA DISK 

C 

PRINT 29* ITESTC 
IF (IC.EQ*0) GO TO 13 
REWIND 12 
K*0 


DK 81 
DK 82 
OK 83 
DK 84 
DK 85 
DK 86 
DK 87 
DK 08 
DK 89 
DK 90 
DK 91 
DK 92 
DK 93 
DK 94 
DK 95 
OK 96 
DK 97 
DK 98 
DK 99 
DK 100 
DK 101 
DK 102 
DK 103 
DK 104 
DK 105 
DK 106 
DK 107 
DK 108 
DK 109 
DK 110 
DK 111 
DK 112 
DK 113 
DK 114 
DK 115 
DK 116 
DK 117 
DK 118 
DK 119 
DK 120 
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LISTING OP DECK* DATADK 


PAGE 4 


CARD NO. 










121 


DO 12 I-1*IC 






DK 

121 



READ (12) IRUN*JS 






OK 

122 



K-K+l 






DK 

123 



NPRT ( K)*IRUN 






DK 

124 

125 


IF (K-10) 12*11*11 






DK 

125 


11 

PRINT 30* (NPRT(J)*J»1*10) 






DK 

126 



K-0 






OK 

127 


12 

CONTINUE 


* 




DK 

128 



IF (K.GT.0.AND.K.LT.10) PRINT 

30* (NPRT( J>* J«1*K) 



DK 

129 

130 


READ (12) IRUN*K 






DK 

130 



K-K-l 






DK 

131 



PRINT 31* K 






DK 

132 



REWIND 12 






OK 

133 



GO TO 28 



1 



DK 

134 

135 

C 







DK 

135 


C 

PRINT ERROR THAT MAY HAVE OCCURRED 1 

DURING WRITE TO 

DISK 


DK 

136 


C 







DK 

137 


13 

PRINT 32 






DK 

138 



IERR-1 






DK 

139 

140 


GO TO 20 






DK 

140 


14 

PRINT 33* NPMAX 






DK 

141 



IERR*1 






DK 

142 



GO TO 28 






DK 

143 


15 

PRINT 34* ITEST* ITESTC 






DK 

144 

145 


IERR-1 






DK 

145 



GO TO 28 






DK 

146 


C 







DK 

147 


C 

READ DATA FROM DATA DISK 






DK 

148 


C 







DK 

149 

150 

C 

READ TAPE 12 TO FIND POINT 

COUNT 

OF 

FIRST DATA 

POINT 

FOR 

DK 

150 


C 

DESIRED RUN 






DK 

151 


16 

READ (12) J*K 






DK 

152 



IF ( EOF ( 12 ) ) 23*17 






DK 

153 


17 

IF (J-IRUN) 16*18*16 






DK 

154 

155 

C 

READ TAPE 12 TO FIND POINT 

COUNT 

OF 

FIRST DATA 

POINT 

FOR 

DK 

155 


C 

NEXT RUN ON DATA DISK 



* 



DK 

156 


18 

READ (12) J*JS 






DK 

157 



IF (E0F(12>) 24*19 






DK 

158 


C 

COMPUTE NUMRER OF DATA POINTS ON 

DATA DISK FOR 

DESIRED RUN 

DK 

159 

160 

19 

NP-JS-K 






DK 

160 
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LISTING OF DECK* DAT AD K PAGE 5 

CARD NO. 


161 


IF (NP.GT.NDMAX) GO TO 26 

DK 

161 


C 

READ DATA FROM DATA DISK 

DK 

162 



DO 22 1*1# NP 

DK 

163 



CALL READMS (NDISK#DATA# I3#K ) 

DK 

164 

165. 

c , 

CHECK TEST NUMBER 

DK 

165 



ITESTC«I FIX (DATA (2)+0. 0001) 

DK 

166 



IF (ITESTC.NE.ITEST) GO TO 25 

DK 

167 


c 

STORE DATA NEEDED FOR PLOTTING 

DK 

168 



DO 21 J *1# NVAR 

DK 

169 

170 


J 3*NL0C t J ) +3 

DK 

170 



IF ( J 3-13 ) 20# 20# 27 

DK 

171 


20 

VAR(I#J)*DATA(J3) 

DK 

172 


21 

CONTINUE 

DK 

173 


22 

K-K+l 

DK 

174 

175 


GO TO 28 

DK 

175 


C 


DK 

176 


C 

PRINT ERROR THAT MAY HAVE OCCURRED DURING READ FROM nATA DISK 

DK 

177 


c 


DK 

178 


23 

PRINT 35# IRUN 

DK 

179 

180 


IERR-2 

DK 

180 



GO TO 28 

DK 

181 


24 

PRINT 36# IRUN 

DK 

182 



I ERR*1 

DK 

183 



GO TO 28 

DK 

184 

185 ! 

25 

PRINT 37# ITE$T#ITESTC 

DK 

185 



IERR-1 

DK 

186 



GO TO 28 

DK 

187 


26 

PRINT 38# IRUN#NDMAX 

DK 

188 



IERR*1 

DK 

189 

190 


GO TO 28 

DK 

190 


27 

PRINT 39# JS# IRUN# IM AX 

DK 

191 



IERR-1 

DK 

192 


C 


OK 

193 


C 

RETURN TO CALLING PROGRAM 

DK 

194 

195 

c 


DK 

195 


28 

RETURN 

DK 

196 


C 


DK 

197 


29 

FORMAT <1H1#5X#71HTHE FOLLOWING LIST OF RUNS ARE AVAILABLE ON THE 

DK 

198 



1DATA DISK FOR TEST NO. #15/) 

DK 

199 

200 

30 

FORMAT (5X#10T10) 

DK 

200 
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LISTING CF DECK* DATADK 


PAGE 


6 


CARO NO. 


201 

205' 

210 

215 

220 


31 

32 

33 


34 

35 

36 

37 

38 

39 


FORMAT (/5X,29HT0TAL NUMBER OF DATA POINTS 17 ) 

FORMAT ( /5X, 67HER ROR OCCURRED DURING WRITE TO DATA DISK NO OA 

IT A ON INPUT TAPE) 

FORMAT ( /5X, 92HERR0R OCCURRED DURING WRITE TO DATA DISK POINT 

1 COUNT HAS EXCFFDED ALLOWABLE MAXIMUM OF ,17) 

FORMAT ( / 5X, 71HE RROR OCCURRED OURING WRITE TO DATA DISK REOUF 

1STFD TEST NUMRFR IS >I5,34H AND TEST NUMBER ON INPUT TAPE IS ,15) 

FORMAT ( /5X,5PHERR0R OCCURRED DURING READ FROM DATA DISK RUN 

INUMBER ,17, 15H IS NOT ON DISK) 

FORMAT C /5X, 91HERR0R OCCURRED DURING READ FROM DATA DISK SEAR 

1CH FOR STARTING POINT FOR RUN FOLLOWING ,I7,14H WAS NOT FOUND) 

FORMAT { / 5X, 72HERR0R OCCURRED DURING READ FROM DATA DISK REOU 

1ESTED TEST NUMBER IS , I5,?8H AND TEST NUMBER ON DTSK IS ,15) 

FORMAT ( /5X,47HERR0R OCCURRED DURING READ FROM DATA DISK /5V, 

125HNUMBER OF POINTS FOR RUN ,T7,30H EXCEEDS ALLOWABLE MAXIMUM OF , 
214) 

FORMAT ( / 5X, 47HERR0R OCCURRED DURING READ FROM DATA DISK /5X, 

118HVARIABLE LOCATION ,I4,23H REQUESTED FOR RUN NO. ,I7,30H EXCEEDS 

2 ALLOWABLE MAXIMUM OF ,14) 

END 


DK 201 
DK 202 
DK 203 
DK 204 
DK 205 
DK 206 
DK 207 
DK 208 
DK 209 
DK 210 
DK 211 
DK 212 
DK 213 
DK 214 
DK 215 
DK 216 
DK 217 
DK 21B 
DK 219 
OK 220- 
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LISTING OF DECK* EDIT PAGE 1 

CARD NO. 


1 


SUBROUTINE EDIT ( T, X , Y, N, NS YM, I S IZE# I RUN, TOLR, KEDT, I FRR ) 

ED 

1 


C 


ED 

2 


c 


ED 

3 


c 

ROUTINE TO SORT AND FDIT DATA 

ED 

A 

5 

c 


ED 

5 


c 

CODED BY — HARRY L. MORGAN NASA/LARC/TAD/AAB 1983 

ED 

6 


c 


ED 

7 


c 

PARAMETER DEFINITION 

ED 

8 


c 

T - ARRAY CONTAINING INDEPENDENT VARIABLE 

ED 

9 

10 

c . 

X - ARRAY CONTAINING X VARIABLE 

ED 

10 


c 

Y - ARRAY CONTAINING Y VARIABLE 

ED 

11 


c 

N - NUMBER OF POINTS TO BE SORTED AND EDITED 

ED 

12 


c 

NS YM - SYMBOL NUMBER TO BE USED TO PLOT UNEDITFD POINTS 

ED 

13 


c 

ISIZE - SYMBOL SIZE 

ED 

1A 

15 

c 

IRUN - RUN NUMBER 

ED 

15 


c 

TOLR - EDIT TOLERANCE OF INDEPENDENT VARIABLE 

ED 

16 


c 

KEDT - NUMBER OF RUNS ALREADY EDITED AND STORED FOP FUTURE 

ED 

17 


c 

PRINT OUT 

ED 

18 


c 

IERR - EDIT ERROR INDICATOR 

ED 

19 

20 

c 

IERR«0. NO ERRORS OCCURRED 

ED 

20 


c 

IERR-1, DURING EDIT OF FIRST SORT POINT ALL GREATER 

ED 

21 


c 

POINTS ARE WITHIN EDIT TOLERANCE 

ED 

22 


c 

IERR-2, DURING EDIT OF LAST SORT POINT ALL LESSER 

ED 

23 


c 

POINTS ARE WITHIN EOIT TOLERANCE 

ED 

2 A 

25 

c 

I ERR*3 , DURING EDIT OF INTERMEDIATE SORT POINT ALL 

ED 

25 


c 

LESSER POINTS ARE WITHIN EDIT TOLERANCE 

ED 

26 


c 

I ERR* A, DURING EDIT OF INTERMEDIATE SORT POINT ALL 

ED 

27 


c 

GRFATFR POINTS ARE WITHlfo EDIT TOLERANCE 

ED 

28 


c 


ED 

29 

30 

c 

T 9 X# AND Y MUST BE PROPERLY DIMENSIONED IN CALLING PROGRAM 

ED 

30 



DIMENSION T ( 1) * X<1>, Y(l) 

ED 

31 


c 

WORK ARRAYS 

ED 

32 



COMMON /WORK/ TIN ( 50 ), TEDT( 50 ) , TSAVE < 50) * YSAVE ( 50 )> IEDT ( 50 ), ITEHP C 

ED 

33 



150 )* I SAVE ( 50 ) 

ED 

3A 

35. 

c 


ED 

35 


c 

CHECK TO SEE IF THIS RUN HAS ALREADY BEEN EDITED AND STORED FOR 

ED 

36 


c 

FUTURE PRINT OUT 

ED 

37 


c 


ED 

38 



IERR-0 

ED 

39 

A0 ; 


IF (N.LT.l) RETURN 

ED 

AO 
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LISTING CF DECK* EDIT 


PAGE 2 


CARO NO. 


41 


NSTORE-1 

FD 

41 



REWIND 11 

ED 

42 



IF (KEPT) 3,3,1 

ED 

43 


1 

READ (11) NRUN 

ED 

44 

45 


IF ( EOF ( 11 ) ) 3,2 

FO 

45 


2 

IF (IRUN-NRUN) 1,4,1 

ED 

46 


3 

KEDT-KEDT+1 

ED 

47 



NSTORE-O 

ED 

48 



IF CKEHT.GT.l) BACKSPACE 11 

ED 

49 

50 


WRITE (11) IRON 

ED 

50 



END FILE 11 

ED 

51 


4 

IN»N 

ED 

52 


C 


ED 

53 


C 

SET INITIAL VALUES OF TIN, TEDT, I EDT , AND ITEMP 

ED 

54 

55 

c 


ED 

55 



DO 5 I«1,N 

ED 

56 



TIN(I )-T(I) 

ED 

57 



TEDT ( I )«0.0 

ED 

58 



IEDT ( I )»0 

ED 

59 

60’ 

5 

ITEMP ( I ) *1 

ED 

60 


C . 


ED 

61 


c 

SORT OATA 

ED 

62 


c 


ED 

63 



K1»N 

ED 

64 

65’ 


DO 12 I-1,N 

ED 

65 



IF (I.EQ.l) GO TO 7 

ED 

66 



K*0 

ED 

67 



DO 6 JJ-1,K1 

ED 

68 



IF ( JJ.EQ.JSAVE) GO TO 6 

ED 

69 

70 


K-K+l 

ED 

70 



TEDT(K)«T(JJ) 

ED 

71 



IEDT(K)-ITEMP( JJ) 

ED 

72 


6 

CONTINUE 

ED 

73 



TMIN-TEDT(l) 

ED 

74 

75 


Kl-Kl-1 

ED 

75 



GO TO 8 

EO 

76 


7 

TMIN-T(l) 

ED 

77 


8 . 

JSAVE-1 

ED 

78 



DO 11 J«1»K1 

ED 

79 

80' 


IF (I.EQ.l) TEDT(J)-T(J) 

ED 

80 
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LISTING OF DECKt EDIT 


° AGE 


3 


CARD NO. 
81 

85 

90 

95 

100 

105 

110 

115 

120 



IF (T.EO.l) IEDT(J)-ITEMP( J) 

ED 

81 


IF ( TEDTt J ) • LE .TM IN ) GO TO 9 

ED 

82 


GO TO 10 

ED 

83 

9 

J SAVE *J 

ED 

84 


TMIN«TEDT(J ) 

ED 

85 

10 

T< J)-TEDT(J ) 

ED 

86 


ITEMP (J)-IEDT(J) 

ED 

87 

11 

CONTINUE 

ED 

88 


I$AVE(I)-IEDT( JSAVE) 

ED 

89 

12 

TSAVE (I ) »TMIN 

ED 

90 


DO 13 I*l> IN 

ED 

91 


Kl-ISAVE(I) 

ED 

92 

13 

YSAVE (I )»Y(K1) 

ED 

93 

C 


EO 

94 

C 

EDIT CATA 

ED 

95 

C 


ED 

96 


Nl-N-1 

ED 

97 


IF (Nl.LT.l.OP.TOLP.LE.O.O) GO TO 31 

ED 

98 


K1 *0 

ED 

99 


00 15 I-1#N1 

ED 

100 


ITENP ( I ) *0 

ED 

101 


OIFF-ABS(TSAV c (I+l)-TSAVEfI) ) 

ED 

102 


IF (DIFF.LE. TOLR) GO TO 14 

ED 

103 


GO TO 15 

ED 

104 

14 

ITEMP( I)-I 

ED 

105 


K1 *1 

ED 

106 

15 

CONTINUE 

ED 

107 


IF ^(Kl.EO.O ) GO TO 31 

ED 

108 


N-l 

ED 

109 


I»1 

ED 

110 

16 

Kl-I 

ED 

111 


J*0 

ED 

112 


IF (I.GT.N1) go TO 27 

ED 

113 


IF (ITEHP m.NF.01 GO TO 17 

ED 

114 


IE DT (N ) « IS AVE (I ) 

ED 

115 


N-N41 

ED 

116 


I-I+l 

ED 

117 


I EDT ( N ) -ISA VE( I ) 

ED 

118 


GO TO 16 

ED 

119 

17 

I-I+l 

ED 

120 
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CARO NO. 


121 

125 

130 

135 

140 

145 

150 

155 

160 


J*J + 1 

IF (I .EO.IN) GO TO 21 
IF (ITEHP(I).EO.O) GO TO 18 
GO TO 17 

18 K2-K14J 

. . IF (Kl.EQ.l) GO TO 19 

11- Kl-1 

12- K241 
GO TO 23 

19 II -K2 

20 11-1141 
12-1141 

IF (I2.GT.IN) IFRR-1 
IF (I2.GT.IN) GO td 29 
OIFF-ABS <TSAVF(I2)-TSAVF< ID ) 

IF CDIFF. LE.TOLP) GO TO 20 
GO TO 23 

21 I2-K1 

22 12-12*1 
II -1 2-1 

IF (Il.LT.l) IERR-2 
IF (Il.LT.l) GO TO 29 
0IFF-ABS(TSAVE(I2)-TSAVE( ID ) 

IF (DIFF.LE.TOLP) GO TO 2? 

k?«IM + I 

23 IF (Il.LT.l) TERR-3 

IF (Il.LT.l) GO TO 29 

IF (I2.GT.IN) IERR-4 
IF (I2.GT.IN) GO TO 29 
I EOT ( N)-ISAVE(Kl) 

T(1>-YSAVE(I1)*T$AVE(I2)-YSAVE(I2)*TSAVE<I1> 
T(2)-YSAVE( 1 2)-YSAVF (ID 
T(3)-TSAVE(I2)-TSAVE (ID 

TMIN-ABS(YSAVE(Kl)-(T(l)4T(?)*TSAVE(Kl))/t(3) ) 

K1-K141 

DO 25 K-K1>K2 

DIFF-ABS ( YS AVE (K)-(T (1 )4T(2)*TSAVE(K) ) /T ( 3 ) ) 

IF (DIFF.LE.THlN) GO TO 24 
GO TO 25 

24 THIN-DIFF 


ED 121 
ED 122 
ED 123 
ED 124 
ED 125 
ED 126 
ED 127 
ED 12B 
ED 129 
ED 130 
ED 131 
ED 132 
ED 133 
ED 134 
ED 135 
ED 136 
ED 137 
ED 138 
ED 139 
ED 140 
ED 141 
ED 142 
ED 143 
EO 144 
ED 145 
ED 146 
ED 147 
ED 148 
ED 149 
FO 150 
ED 151 
ED 152 
ED 153 
ED 154 
ED 155 
ED 156 
ED 157 
ED 158 
ED 159 
ED 160 
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161 


IEDT(N)«ISAVE(K) 

ED 

161 


25 

CONTINUE 

ED 

162 



IF (K2.EQ.IN) GO TO 27 

ED 

163 



IF (I.EO.Nl) GO TO 26 

ED 

164 

165 


N-N4-1 

ED 

165 



I-I+l 

ED 

166 



GO TO 16 

ED 

167 


26 

N-N + l 

ED 

168 



IEDT<N)«ISAVE(IN) 

ED 

169 

170 

27 

DO 28 K*1>N 

ED 

170 



Kl-IEDT(K) 

ED 

171 


28 

TEDT(K)-TIN(K1) 

ED 

172 



GO TO 33 

EO 

173 


29 

IF (NSTORE* EQ#0) PRINT 40* I RUN*TOLR> IERR 

ED 

174 

175 


DO 30 I-1*IN 

ED 

175 



Kl-ISAVE(I) 

ED 

176 



CALL PNTPLT {X ( K1 ) > Y (K 1 ) , NSYH, ISTZE > 

ED 

177 



IF (NSTORE.EO.O) PRINT 41* I* KljTSA VE( I ) * X(K1 )* Y <K1 ) 

ED 

178 


30 

CONTINUE 

EO 

179 

180. 


RETURN 

ED 

180 


31 

DO 32 I»1*I N 

ED 

181 



Kl-ISAVE(I) 

ED 

182 



TEDT (I) •TINCK1 ) 

ED 

183 


32 

I EOT { I )*K1 

ED 

184 

185 

C 


ED 

185 


C 

WRITE EDITED DATA ON TAPE 7 FOR USE IN SUBROUTINE PRTEDT 

ED 

186 


C 


ED 

187 


33 

IF CNSTORE.EQ.l) GO TO 34 * 

ED 

188 



WRITE (7) IRUN* IN*N* TOLR 

ED 

189 

190 


WRITE (7) (TIN(I)*I-1*IN) 

ED 

190 



WRITE (7) ( IS AVE (I)*TSAVE(I)#I*1#TN) 

ED 

191 



WRITE (7) (IEDT{I»,TEDT(I)*I-1*N) 

EO 

192 


C 


ED 

193 


c 

REORDER T* X* AND Y FOR PLOTTING 

ED 

194 

195 

c 


EO 

195 


34 

DO 35 1*1# IN 

ED 

196 


35 

T(I)«TEDT(I) 

ED 

197 



DO 36 I«l» IN 

EO 

198 



TSAVECI)-X(I) 

EO 

199 

200^ 

36 

TEDT ( I )*Y ( I ) 

ED 

200 
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201 

205 

210 

215 

220 

225 


OF DECKi EDIT 


37 

C 

c 

c 


38 


39 
C 

C 

C 

C 

C 

40 


41 


DO 37 I«1,IN 
K-IEDT(I) 

X(I)-TSAVE(K) 

Y( I)-TEDT(K) 

PLOT DATA POINTS NOT EDITED 

IF (N.EQ.IN) RETURN 

DO 39 1*1* IN 

Kl-0 

DO 38 J»l,N 

IF (IEOT(J).EO.I) Kl-1 

CONTINUE 

IF (KliEOaO) CALL PNTPLT (TSAVE(T)#TECT(I) #NSYM# IST7E) 
IF ((Kl.EQ.O). AND. (N STOP E.FO*0)) WPITF (7) I#TIN { T ) 
CONTINUE 


RETURN TO CALLING PROGRAM 
RETURN 


C //2X#19HATTFHPT TO 
■#F8*3#10X#9HERR0R f 


FORMAT 
1LERANCE 
21HY/) 

FORMAT ( 215# 3^10 . 4 ) 
END 


EDTT RUN# 17# IX# 29HFAILE0 F PR AN EDIT 


PAGE 

6 

ED 

201 

ED 

202 

ED 

203 

ED 

204 

ED 

205 

ED 

206 

ED 

207 

FD 

208 

ED 

209 

FD 

210 

ED 

211 

ED 

212 

FD 

213 

FD 

214 

ED 

215 

ED 

216 

ED 

217 

FD 

218 

ED 

219 

FD 

220 

ED 

221 

FD 

222 

TO ED 

223 

>X# ED 

224 

FO 

225 

ED 

226 

ED 

227- 
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CARD NO * 


1 SUBROUTINE PRTEDT PR 1 

C PR 2 

C ROUTINE TO PRINT EDITED DATA STORFD ON TAPE 7 PR 3 

C PR A 

5 C CODED BY -- HARRY l • MORGAN NASA/L ARC/TAO/AAB 1983 PR 5 

C PR 6 

C WORK ARRAYS PR 7 

COMMON /WORK/ TI N ( 50 ) , I S A VE < 50 ) , TS A VE ( 50 ) , IE DT ( 50 ) , TEDT ( 50) , TUN ( 50 PR 6 

1) , TUN ( 50 ) PR 9 

10 C p R 10 

END FILE 7 PR 11 

REWIND 7 PR 12 

PRINT 7 PR 13 

1 READ (7) IRUN, TN,N,TOLR PR 14 

15 IF ( E OF (7 ) ) 6,2 PR 15 

2 NR* IN-N PR 16 

READ (7) (TIN(I),I*1,IN) PR 17 

READ (7) (ISAVF(I >,TSAVE(I)f I«1,IN) PR 18 

READ (7) ( IEDT( I) >TFDT(I >, I-l.N) PR 19 

20 IF (NR.LE.O) GO TO A PR 20 

DO 3 1*1, NR PR 21 

3 READ (7) IUN(T),TUN(T) PR 22 

A PRINT 0, IRUN, TOLR PR 23 

DO 5 1*1, IN PR 2 A 

25 IF < ( I.LE.N) .AND. (T.LF.NR) ) PRINT 9, I, TIN { I ), IS AVE ( I ) , TS AVF ( I ) , IE PR 25 

IDT { I ) #TE DT ( I ) j. IUN<I) ,TUN(I) PR 26 

IF ( (I.LE.N). AND. (I. GT.NP) ) PRINT 10, I, TI N ( I ) , I SAVF ( T ) , TS AVF ( I ) , I PR 27 

1FDT ( I ), T FDT ( I ) PR 28 

IF ( ( I.GT.N).AND. (I. LF.NR ) ) PRINT 11, I , TIN ( I) t ISAV F ( I ) , TS AVF ( I ) , I PR 29 

30 1UN (I ) , TUN ( I ) PR 30 

IF ( (I.GT.N).AND. (I.GT.NR )) PRINT 12, I, TI N ( I ) , ISAVF ( I ) , TS A VE < T ) PR 31 

5 CONTINUE PR 32 

GO TO 1 PR 33 

6 RETURN PR 3A 

35 C p R 35 

7 FORMAT ( 1H1, 10X,35HTHE FOLLOWING RUNS HAVE BEEN EDITED) PR 36 

8 FORMAT ( //IX, 9 HR UN NO. *, 17, 10X, 1 INTOLERANCE *,F8.3//9X,5H INPUT, 19 PR 37 

IX, 6HS0RTED, 19X,6HEDI TFD , 20X, 7HDEL ETFD // A ( IX, 9HP0INT NO. , 3X, 7HT-V AL PR 38 

2UE,5X )/) PR 39 

AO 9 FORMAT (A(5X,T3,2X,F10.A,5X) ) PR AO 
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41 


10 

11 

12 


FORMAT (3(5X#I3#2*#F10.4#5X>) 

FORMAT (2(5X#I3*2XfF10»4j5X)>30X#I3>2X*F10.4) 
FORMAT (2(5X,I3,2X,F10.4, 5X) ) 

END 


PR 41 
PR 42 
PR 43 
PR 44 
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LISTING OF PFCKt ZFROLN 
CARD NO* 


PAGE 


1 


1 

5 


10 


15 


20 


25 


30 


35 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1 


SUBROUTINE ZEROLN ( X 0, YO* DI ST» LI NE > 

ROUTINE TO DRAW ZERO LINFS 

CODED BY — HARRY L* MORGAN NASA/LARC/TA^/ AAR 1983 

PARAMETER DEFINITION 

XO t YO - STARTING LOCATION OF 7ER0 LINE 
DIST - LENGTH OF 7FP0 LINE 
LINE - AXIS OPTION 

L INF- 1, ZERO LINE PARALLEL TO X-AXIS 
LINE-2, 7FR0 LINE PARALLEL TO Y-AXIS 

NOTE - THICKNESS OF 7FR0 LINE - 2.*D 
0-0*014 

CALL CALPLT (X0,Y0,3) 

IF (LINE.EQ.2) GO TO 1 
X-XO+DIST 

CALL CALPLT <X,Y0>2) 

Y-YO+D 

CALL CALPLT (X,Y,2) 

CALL CALPLT (X0#Y,2) 

Y-YO-O 

CALL CALPLT <XO,Y,R) 

CALL CALPLT (X,Y,2) 

CALL CALPLT (X,Y0,2) 

GO TO 2 
Y- YO+DIST 

CALL CALPLT (X0,Y,2) 

X-XO+D 

CALL CALPLT (X,Y,2) 

CALL CALPLT (X,Y0,2> 

X-XO-D 

CALL CALPLT <X,Y0,2) 

CALL CALPLT <X,Y,2) 

CALL CALPLT (X0,Y,2> 

CALL CALPLT <XO,YO,3) 

RETURN 

END 


ZN 1 
ZN 2 
ZN 3 
ZN 4 
7N 5 
ZN 6 
ZN 7 
ZN 8 
7N 9 
7N 10 
ZN 11 
ZN 12 
ZN 13 
ZN 14 
ZN 15 
ZN 16 
ZN 17 
ZN IB 
ZN 19 
ZN 20 
ZN 21 
7N 22 
ZN 23 
ZN 24 
ZN 25 
ZN 26 
ZN 27 
ZN 28 
ZN 29 
ZN 30 
ZN 31 
ZN 32 
ZN 33 
ZN 34 
ZN 35 
ZN 36 
ZN 37 
ZN 38 
ZN 39- 
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PAGE 


1 


1 


5 


10 


15 


20 . 


25 


30 


35 


40. 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

1 

2 


SUBROUTINE CURPLT <T,X, Y,n, IVAR, NS YM, ISI ZE, IOP, IRON. TENS ION ) 
ROUTINE TO PLOT AND FAIR DATA 

CODED BY — HARRY L. MORGAN NASA/L ARC/TAD/AAB 1983 

PARAMETER DEFINITION 

T - ARRAY CONTAINING INDEPENDENT VARIABLE 
X - ARRAY CONTAINING X VARIABLE 
Y - ARRAY CONTAINING Y VARIABLE 
N - NUMBER OF PntNTS TO BF PLOTTED 

IVAR - VARIABLE CODE IVAR-0 IF X AND T ARE THE SAME 

IVAR-1 IF X AND T ARE DIFFERENT 

NS YM - SYMBOL NUMBER 
IS IZE - SYMROL SIZE 

IOP - PLOTTING OPTION I0P*0 PLOT SYMBOLS ONLY 

I0P»1 PLOT CURVE AND SYMBOLS 

IRUN - RUN NUMBER 
TENSION - SPLINE TENSION 

T * X* AND Y MUST BE PROPERLY DIMENSIONED IN CALLING ROUTINE 
DIMENSION T( 1)# X<1), Y<1), HC3) 

WORK ARRAYS 

COMMON /WORK/ DS1 (105 ) ,D$2< 105) ,DUMX ( 105 ), DUMY < 105) 

COMMON /PL T / MX(50), MY (50), AC50) 

REAL MX, MY 

ROUTINE CONSTANTS 

DATA NH/101/,XNT/100./,EP$/0.0001/,NTA/100/,H/.02,.03,.04/ 

SELECT PLOTTING OPTION 

IF (I0P.EQ.0.0R.N.LT.3) GO TO 1 
GO TO 3 

PLOT SYMBOLS ONLY 
DO 2 I-1,N 

CALL PNTPLT (X(I), Y(I ) ,NSYH, I SIZE ) 

RETURN 


CU 

1 

cu 

2 

CU 

3 

cu 

A 

cu 

5 

cu 

6 

cu 

7 

cu 

8 

cu 

9 

cu 

10 

cu 

11 

cu 

12 

cu 

13 

cu 

14 

cu 

15 

cu 

16 

cu 

17 

cu 

18 

cu 

19 

cu 

20 

cu 

21 

cu 

22 

cu 

23 

cu 

24 

cu 

25 

cu 

26 

cu 

27 

cu 

28 

cu 

29 

cu 

30 

cu 

31 

cu 

32 

cu 

33 

cu 

34 

cu 

35 

cu 

36 

cu 

37 

cu 

38 

cu 

39 

cu 

40 
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CARD NO, 


41 

C 


CU 

41 


C 

PLOT CURVE AND SYMBOL 

cu 

42 


C 


cu 

43 


C 

CHECK TO SEE IF X OR T IS STRICTLY INCREASING 

cu 

44 

45 

3 

DO 5 I«2#N 

CU 

45 



IF (IVAR.EO.l) GO TO 4 

cu 

46 



IF (Xm.lT.XII-m GO TO 6 

cu 

47 



GO TO 5 

cu 

4P 


4 

IP (Till .LT.T(I-l)) GO TO 6 

cu 

49 

50 

5 

CONTINUE 

cu 

50 



GO TO 7 

cu 

51 


C 

PRINT ERROR MESSAGE IF X OR T NOT STRICTLY INCREASING 

cu 

52 


6 

IF (IVAR.EO.O) PRINT 26* IRUN 

cu 

53 



IF (IVAR.EO.O) PRINT 27, ( X ( I ) # Y( I) # I »1# N> 

cu 

54 

55 


IF (IVAR.EO.l) PRINT 28# IRUN 

cu 

55 



IF (IVAR.EO.l) PRINT 29# (T(I)#X(I)#Y(I)#I*1#N) 

cu 

56 


C 

PLOT SYMBOLS ONLY IF X OR T NOT STRICTLY INCREASING 

cu 

57 



GO TO 1 

cu 

58 


C 

FIT SPLINE CURVE THROUGH DATA POINTS 

cu 

59 

60 

7 

IF (IVAR.EO.O) CALL CUBSPL ( X, Y# N, MX, TENS I ON# TENS# A ) 

cu 

60 



IF (IVAR.EO.l) CALL CUBSPL ( T# X#N# MX# TENSION# TENS# A ) 

cu 

61 



IF (IVAR.EO.l) CALL CUBSPL (T# Y#N# MY# TENS ION # TENS# A ) 

cu 

62 


C 

PLOT FIRST POINT 

cu 

63 



CALL PNTPLT (X(l)# Y( 1)#NSYM#ISIZE) 

cu 

64 

65 

C 

PLOT AND FAIR REMAINING POINTS 

cu 

65 



NN«N-1 

cu 

66 



DO 25 I* 1# NN 

cu 

67 


C 

COMPUTE STRAIGHT LINE DISTANCE BETWEEN TWO POINTS. IF DISTANCE 

cu 

68 


c 

LESS THAN SYMBOL DIAMETER# PLnT POINTS ONLY. 

cu 

69 

70 


Xl-X( I+1)-X(I) 

cu 

70 



Yi-Y(in)-Yd) 

cu 

71 



DS»SGRT(X1*X1*Y1*Y1) 

cu 

72 



T1*AT ANF ( Yl# XI ) 

cu 

73 



CALL SYMBS (NSYM#ISIZE#XS1#YS1#T1) 

cu 

74 

75 


DSS1«SQRT(XS1*XS1+YS1*YS1) 

cu 

75 



XI— XI 

cu 

76 



Yl— Yl 

cu 

77 



T1-ATANF(Y1#X1) 

cu 

78 



CALL SYMBS (NSYM# ISI ZF# XS2* YS2#T1) 

cu 

79 

80 


DS$2«SQRT(XS2*XS2+Y$2*YS2) 

cu 

80 
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CARD NO. 
01 

05 

90 

95 

100 

105 

no' 

115. 

120 


IF ((DSS1+DSS21.GE.DS) GO TO 24 

C COMPUTE DISTANCE ALONG CURVE AS A FUNCTION OF X OR T BETWEEN 

C ‘ POINT I AND 1+1 

IF (IVAR.EO.O) NT-IFIX(XNT*tX(I+l)-X(I)) )+l 
IF (IVAR.EQ.l) NT-IFIX(XNT*(T(I+1)-T(I)))+1 
IF (NT.LT.3) NT-3 
IF CNT.GT.NH) NT-NH 

IF (IVAR.EQ.O) DELTA-(XCI41)-XCI))/FL0AT(NT-1) 

IF (IVAR.EQ.l) DELTA-(T(I+1)**T(I) l/FLOAT (Nt-1 ) 

DUMX(1)-X(I) 

DUMY( 1)-Y( I ) 

DS1 (1 1-0.0 

IF (IVAR.EQ.O) D A-X ( I ) 

IF (IVAR.EQ.l) DA-T(I) 

C 

DO 9 J-2# NT 
DA-DA+DELTA 

IF (IVAR.EQ.l) GO TO B . 

DUMX ( J )-DA 

0UMY(J)-FUNC(DA,X(I),X(I+l)*Y(I)>Y(U-l)*MX(I)*MX(l4l)>TENS) 

GO TO 9 

0 DUf*Y( J l-FUNC (DA* T ( I ) »T ( I+l )# Y( I ) » Y ( I + D# HY ( I )# MY(I + 1 )# TENS) 

DUMX ( J l-FUNC (DA*T(I)»T(I+1)#X(I)#X( I+l)# MX ( I)*MX(I41)*TEN$) 

9 DSKJ )-SQRT((DUMX(J)-X(I) )**2+( OUMY ( J )-Y ( I) )**2) 

C 

DO 10 J-1*NT 
K- NT41-J 

10 DS2 ( J l-SQRT ( ( DUM X(K)-X(I+l))**2+( OUMY ( K) -Y ( 1 + 1 ) )**2) 

C FIND X AND Y LOCATION WHERE SYMBOL AND CURVE INTERSECT 

DELTA-H(ISIZE) * 

IF (IVAR.EQ.O) DA-X(I) 

IF (IVAR.EQ.l) DA-T(I) 

C 

DO 14 J-2*NTA 
DA-DA4DF.LTA 

IF (IVAR.EQ.l) GO TO 11 
IF (DA.GE.X (I+l) ) GO TO 15 
Xl-DA-X(I) 

Y1-FUNC(DA*X(I)*X(I+1)*Y(I)*Y(I41)*MX(I)*MX(I+1)*TENS)-Y(I) 

GO TO 12 


CU 81 
CU 82 
CU 83 
CU 84 
CU 85 
CU 86 
CU 87 
CU 88 
CU 89 
CU 90 
CU 91 
CU 92 
CU 93 
CU 94 
CU 95 
CU 96 
CU 97 
CU 98 
CU 99 
CU 100 
CU 101 
CU 102 
CU 103 
CU 104 
CU 105 
CU 106 
CU 107 
CU 108 
CU 109 
CU 110 
CU 111 
CU 112 
CU 113 
CU 114 
CU 115 
CU 116 
CU 117 
CU 118 
CU 119 
CU 120 


i 
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121 

11 

IF ( DA.GE.T ( 1+1 ) ) GO TO 35 

CU 

121 



Y 1-FU NC(DA, T(I),T(I+1),Y( I), Y(I+1)>MY(I ) > MY< I +1 ) , TENS > -Y (I) 

cu 

122 



Xl«FUNC(DA»T(I)»T(I*l)«Xm;X(I+l)»MXm»HX(I*l)»TENS)-Xm 

CU 

123 


12 

DS*SORT(Xl*Xl+Yl*Yl) 

CU 

124 

125 


T 1-ATANF ( Yl> XI ) 

cu 

125 



CALL SYMBS (NSYM,ISI ZF>XS1# YS1# Tl) 

cu 

126 



DSS1-S0RT(XS1+XS1*-YS1*YS1) 

cu 

127 



IF (ABS(DS-DSSl) .LE.EPS) GO TQ 15 

cu 

128 



IF (DS.GT.0SS1) GO TO 13 

cu 

129 

130 


GO TO 14 

cu 

130 


13 

DA-DA-DELTA 

cu 

131 



DELTA-DELTA/2. 

cu 

132 


14 

CONTINUE 

cu 

133 


C 


cu 

134 

135’ 

15 

xsi-xsi+xm 

cu 

135 



YS1-YS1+Y ( I ) 

cu 

136 



DELTA-H(ISIZE* 

cu 

137 



IF (IVAR.EQ.O) 0 A-X ( I +1 ) 

cu 

138 



IF (IVAR.EQ.il OA-TCI+1) 

cu 

139 

140 ! 

C 


cu 

140 



DO 19 J -2* NTA 

cu 

141 



DA-DA-DELTA 

cu 

142 



IF ( I VAR. EQ. 1 ) GO TO 16 

cu 

143 



IF (DA.LE.X (II) GO TO 20 

cu 

144 

145, 


X 1-DA— X{ 1+1) 

cu 

145 



Y1-FUNC(DA>X(I)#X(I+1)>Y(I).Y(I*1).MX(I),MX(I+1),TENS)-Y(T+1) 

cu 

146 



GO TO 17 

cu 

147 


16 

IF (DA.LE.T(I)i GO TO 20 

cu 

148 



Y1«FUNC(DA*T(I)#T(I + 1) »YCn»Y(I+l)»MY(X)#HY(I + l)#TEN$)-Y(I«l) 

cu 

149 

150 


Xl-FUNC (DA»T(I)»T(I+l)jX(I)#X(I+l)#MX(I)#MX(I+l)#TENS)-X(T+l) 

cu 

150 


17 

DS-SQRT(X1*X1+Y1*Y1) 

cu 

151 



T1-ATANF(Y1/X1) 

cu 

152 



CALL SYMBS <NSYM,ISI ZE. XS2> Y$2,T1) 

cu 

153 



DSS2-SQRT (XS2*XS2+YS2*YS2) 

cu 

154 

155 


IF { A BS ( DS-D SS2 ) .LF.EPS) GO TO 20 

cu 

155 



IF (DS.GT.DSS2) GO TO 18 

cu 

156 



GO TO 19 

cu 

157 


18 

DA-DA +DELTA 

cu 

158 



DELTA-DELTA/2. 

cu 

159 

160’ 

19 

CONTINUE 

cu 

160 
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LISTING 
CARD NO 
161 

165 
170 
175 
180 
185 
190 . 
195. 
200 . 


OF DECK t CURPLT 


PAGE 


C 

20 


21 

C 


22 

23 


C 

24 
C 

25 
C 

C 

C 

C 

26 

27 

28 


XS2-XS2 + XU + 1) 

YS2-YS2+YCI+1) 

NP«1 

NDS1-0 

DO 21 4«2*NT 

IF UDSl(J).LT.DSSl) .AND . (NDS1 .EO.O ) ) 
IF UDSl(J).GE.OSSl) .AND. (NDS1.E0.0) ) 
NP-NP41 

DUHX(NP)-DUMX(J) 

DUMY(NP)-DUMYCJ) 

CONTINUE 

PLOT CURVE BETWEEN POINTS 
DUMX ( 1)»XS1 
DUMY(1)»YS1 
DO 22 J-1>NT 

IF (DS2(J) •LE.0SS2) GO TO 22 

GO TO 23 

NP-NP-1 

NP-NP-H 

DUMX (NP )*X$2 

DUMY ( NP ) ■ YS2 

DUMX ( NP + 1 ) *DUMY (N P+1 ) *0*0 
DUMX{NP*2)-DUMY(NP+2)-1.0 
CALL LINE (DUMX# DUMY »NP#1*0*0*0) 

PLOT SYMBOL AT POINT 1+1 
CALL PNTPLT (X ( I ♦! )» Y( 1*1 ) >NSYM# IS IZF ) 
ADVANCE TO NEXT POINT 
CONTINUE 

RETURN TO CALLING PROGRAM 


GO TO 21 
NDS1-1 


RETURN 

FORMAT (1H1//5X,36HX 
l,13X#lHY/> 

FORMAT <5X*2F15.4> 
FORMAT ( 1H1//5X* 36HT 
1,13X#1HX>13X*1HY/) 


IS NOT STRICTLY INCREASING FOR 


IS NOT STRICTLY INCREASING FOR 


CU 161 
CU 162 
CU 163 
CU 164 
CU 165 
CU 166 
CU 167 
CU 168 
CU 169 
CU 170 
CU 171 
CU 172 
CU 173 
CU 174 
CU 175 
CU 176 
CU 177 
CU 178 
CU 179 
CU 180 
CU 181 
CU 182 
CU 183 
CU 184 
CU 185 
CU 186 
CU 187 
CU 188 
CU 189 
CU 190 
CU 191 
CU 192 
CU 193 
CU 194 
CU 195 

RUN* I 7/16X* 1HX CU 196 
CU 197 
CU 198 

RUN,I7/16X»1HT CU 199 
CU 200 
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LISTING CF DECK: CURPLT 
CARO NO. 


PAGE 6 


201 


29 FORMAT (5X*3F15.A) 
END 


CU 201 
CU 202- 
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LISTING OF DECK* ATANF PAGE 1 

CARO NO, 


1 


5 


10 . 


15 



FUNCTION ATANF CDY,DX) 



AF 

1 

c 




AF 

2 

c 

ROUTINE TO COMPUTE ARCTANGENT FOR 

ANGLES FROM 0 TO 360 

DEGREES 

AF 

3 

c 




AF 

A 

c 

COOED BY — HARRY L. MORGAN 

nasa/larc'/tad/aab 

1983 

AF 

5 

c 




AF 

6 


IF (DX,EQ.0,0) GO TO 1 



AF 

7 


ATANF-ABS (DY/DX) 



AF 

8 


ATANF* AT AN (ATANF) *57,2957795131 



AF 

9 


IF < (DX,GT,0.0> ,ANO, (DY.LT.O.O) ) 

ATANF-360, -ATANF 


AF 

10 


I F ( (DX.LT.O.O).AND.(DY.GT.O.O) ) 

ATANF-180. -ATANF 


AF 

11 


IF { (DX,LT ,0,0), AND, (DY,LT ,0,0) ) 

AT ANF*180 ATAVlF 


AF 

12 


RETURN 



AF 

13 

1 

IF (DY.LT.O.O) ATANF-270. 



AF 

1A 


IF (DY.GE,0.0) ATANF-90, 



AF 

15 


RETURN 



AF 

16 


END 



AF 

17 
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LISTING OF DECK! SYMBS 


PAGE 


1 


CARO NO. 


1 

5 


10 


15 


20 


25 


30 


35 


40 


SUBROUTINE SYMBS (NO, I S, X, Y, T) SY 1 

c SY 2 

C ROUTINE TO COMPUTE THE X AND Y INTERCEPTS OF A RADIAL LINE DRAWN SY 3 

C FROM THE CENTER OF A GIVEN SYMBOL WITH THE SIDE OP THAT SYMBOL SY 4 

C SY 5 

C. CODED BY -- HARRY L. MORGAN NASA/L ARC/TAD /AAB 1083 SY 6 

c SY 7 

C PARAMETER DEFINITION SY 8 

C NO - SYMBOL NUMBER SY 0 

C IS - SYMBOL SIZE 1-SMALL 2-MEDIUM 3-LARGE SY 10 

C X AND Y - INTERCEPT OF SYMBOL AND RADIAL DRAWN FROM SYMBOL SY 11 

C CENTER AT ANGLE T. SY 12 

C , T - ANGLE OF RADIAL DRAWN FROM SYMBOL CENTER - SY 13 

c SY 14 

DIMENSION SC ALE( 3) SY 15 

DATA RAD/57. 2957795131/, PI/3. 141502654/ SY 16 

DATA DA/1. 4142135624/, R/l. 7320508076/, SCALE/. 13, .16, .19/ SY 17 


OATA Tl/213. 6900675260/, T2/326. 3099324740/, T3/U6. 5650511771/, T4/3 SY 18 
133.4349488220/, T7/11 . 3099324740/, T8/168. 6900 675260 /, T9/ 218. 6598082 SY 19 
2541/, SI /321. 3401917450/, S2/185. 7105931375/, S3/354. 2 894068625/, S4/1 SY 20 
31 .3099324740/, S5/ 168 .6900675260/, S6/218 . 6598 082541 /, S 7/321. 3401917 SY 21 


4459/ SY 22 

C SY 23 

IF < (NO.EQ.l) .OR. (NO.EQ.ll) ) GO TO 1 SY 24 

IF { (NO.EQ.2J.0R. (N0.E0.12) ) GO TO 2 SY 25 

IF ( (N0.EQ.3).0R.(N0.E0.13)) GO TO 8 SY 26 

IF ( (N0.E0.4).0R. (N0.EQ.14U GO TO 13 SY 27 

IF { ( N0.EQ.5) .OR. (N0.EQ.15) ) GO TO 17 SY 28 

IF ( (N0.EQ.6).0R .(N0.EQ.16) ) GO TO 21 SY 29 

IF ( (N0.E0.7) .OR.fNO .E0.17)) GOTO 26 SY 30 

IF UNO. EQ. 8). OR, (NO. SO. 18)) GO TO 31 SY 31 

IF ((NO.EQ.9).OR,(NO.EQ.19)) GO TO 35 SY 32 

IF ( (N0.EQ.10) .OR. (NO, EQ. 201 ) GO TO 40 SY 33 

IF (<N0.EQ.21).0R.(ND.E0.22)) GO TO 45 SY 34 

IF (N0.GT.22) RETURN SY 35 

C SY 36 

C SYMBOL NUMBER 1 OP 11 SY 37 

C SY 38 

1 X-SCALF(IS)*.5525*C0S(T/RAD) SY 39 

Y*SCALE(IS)*.5525*SIN(T/RAD) SY 40 
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LISTING CF DECK: SYMBS 


PAGE 


2 


41 RETURN SY *1 

C S Y 42 

C SYMBOL NUMRER 2 OR 12 SY 43 

C SY 44 

45 2 IF ((T.GE.0.).AND.(T.LT.45.>) GO TO 3 SY 45 

IF ( (T.GE.45,0) .AND. (T.LT.135.) ) GO TO 4 SY 46 

IF ( (T.GE.135.).AND. (T.LT.225.)) GO TO 5 SY 47 

IF ((T.GE.225.).AND.(T.LT.315.)) GO TO 6 SY 48 

IF ((T.GE.315.).AND.(T,LE.360,)) GO TO 7 SY 49 

50 3 X-SCALE(IS) /2. SY 50 

Y-X*T AN ( T/RAD) SY 51 

RETURN SY 52 

4 Y-SCALEIIS )/2. SY 53 

IF (T.E0.90.) X-0.0 SY 54 

55 IF (T.NE.90.) X- Y /TA N (T/R AD ) SY 55 

RETURN SY 56 

5 X— SCALEt IS>/2. SY 57 

Y-X*TAN(T/RAD) SY 58 

RETURN SY 59 

60 6 Y — SCALEUSW2. SY 60 

IF (T.EQ.270.) X-0.0 SY 61 

IF (T.NE.270.) X-Y/TAN (T/RAD) SY 62 

RETURN SY 63 

7 X-SCALE(IS)/2. SY 64 

65 Y-X*T AN(TZRAD) SY 65 

RETURN SY 66 

C SY 67 

C SYMBOL NUMBER 3 OR 13 SY 68 

C SY 69 

70 8 IF ( (T.GE.O.).AND. (T • LT. 90. ) ) GO TO 9 SY 70 

IF ((T.GE.90.).4ND.<T.LT.1B0.)) GO TO 10 SY 71 

IF ( (T.GE.180. ). AND. (T.LT.270.)) GO TO 11 SY 72 

IF ((T.GE.270.).AND.(T,LE.360.) ) GO TO 12 SY 73 

9 X-SCALE(I$)*DA/?./(TAN(T/RAD)«-l. ) SY 74 

75 Y--X*SCALE(IS)*DA/2. SY 75 

RETURN SY 76 

10 IF (T.E0.90.) X-0.0 SY 77 

IF (T.NE.90.) X-SCALE(I$)*DA/2./ (TAN( T/RAD)— 1 . ) SY 78 

Y-X+SCALE(IS)*DA/2. SY 79 

80 RETURN SY 80 
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LISTING CF DECK » SYMBS 


PAGE 3 


CARD NO* 


81 

11 

X— SC ALE(IS)*DA/2./(TAN(T/RAD)+l.) 

SY 

81 



Y— X-SCALE(IS)*DA/2. 

SY 

82 



RETURN 

SY 

83 


12 

IF ( T *E0* 270* ) X-0.0 

SY 

84 

85 


IF (T.NE.270.) X— SC AL E ( IS ) *DA /2 ./( TAN{ T/RAD )-l . ) 

SY 

85 



Y-X-SCALF(IS)*DA/2* 

SY 

86 



RETURN 

SY 

87 


C 


SY 

88 


C 

SYMBOL NUMBER 4 OR 14 

SY 

89 

90 

C 


SY 

90 


13 

IF ( (T*GE.0.)*AND.(T.LT.90* > ) GO TO 14 

SY 

91 



IF ((T.GE.90.).AND.(T.LT.T1)) GO TO 15 

SY 

92 



IF ((T*GE*T1).AND.(T.LT.T2>) GO TO 16 

SY 

93 


14 

X-(2./3. )*SCALE(IS)*1.105/(TANCT/RAD)+2. ) 

SY 

94 

95' 


Y--2. *X*2. *SC ALE (IS) *1.105/ 3 f 

SY 

95 



RETURN 

SY 

96 


15 

IF (T *E0*90* ) V-0.0 

SY 

97 



IF (T *NE • 90* ) X-(2./3. >*SCALE (I S ) *1. 105/ (TAN (T/ RAO )-2 . ) 

SY 

98 



Y-2.*X+2.*SC ALE (IS)* 1.105/3. 

SY 

99 

XOO’ 


RETURN 

SY 

100 


16 

Y— SCALE( IS )*1. 105/3. 

SY 

101 



IF (T.NE.270.) X-0.0 

SY 

102 



IF (T.NE.270.) X-Y/T AN (T/RAD ) 

SY 

103 



RETURN 

SY 

104 

105 

C 


SY 

105 


C 

SYMBOL NUMBER 5 OR 15 

SY 

106 


C 


SY 

107 


17 

IF ((T.GE.0.).AND.(T.IT.T3)) GO TO 18 

SY 

108 



IF ( (T.GE*T3).AND.(T.LT*225*) ) GO TO 19 

SY 

109 

110 


IF ( (T.GE.225. ). AND* (T.LT.T4) ) GO TO 20 

SY 

110 


18 

IF (T.EQ.90.) X-0.0 

SY 

111 



IF ( T *NE *90 • ) X-SCALE(IS)*1. 22222/3. /(TAN(T/RAD)+1.) 

SY 

112 



Y—X* SC ALE (IS) *1.222 22/3. 

SY 

113 



RETURN 

SY 

114 

115 

19 

X--SC ALE (IS) *1.22222/3. 

SY 

115 



Y»X*TAN (T/RAD) 

SY 

116 



RETURN 

SY 

117 


20 

Y— SC ALE (IS) *1.22222/3. 

SY 

118 



IF (T.EQ.270.) X-0.0 

SY 

119 

120 


IF (T.NE.270.) X-Y/TAN (T/RAD) 

SY 

120 
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LISTING OF DECK* SYMRS 


PAGE 


4 


CARD NO. 


121 


RETURN 

SY 

121 


C 


SY 

12? 


C 

SYMROL NUMBER 6 OR 16 

SY 

123 


C 


SY 

124 

125 

21 

A-4.*SCALE(IS)*1.22222/(3.*PI) 

SY 

125 



B-SCALE(IS)*1.2???2-A 

SY 

126 



T5«ATAN(A/B)*RAD 

SY 

127 



T6*360*-T 5 

SY 

128 



T5-90.+T5 

SY 

129 

130 


IF UT.GE.0.).AND.(T.LT.T5) ) GO TO 22 

SY 

130 



IF ( (T.GE. T5). AND. (T.LT. 225. ) ) GO TO 24 

SY 

131 



IF ( (T.GE. 225. ). AND. (T.LT. T6)> GO TO 25 

SY 

132 


22 

IF (T.EQ.90.) GO TO 23 

SY 

133 



BB«-?.*A*(1.+TAN(T/RAD) ) 

SY 

134 

135' 


AA»TAN(T/RAD)**2+1. 

SY 

135 



CC«2.+A*A-(SCALE(I$)*1.22222)**2 

SY 

136 



X-SQRT(BB*BB-4.*AA*CC)/(2.*AA) 

SY 

137 



IF ( (T.GE.O.) . AND. (T .LT.90. )*) X -X + BB / { 2 .* AA) 

SY 

138 



IF ( (T.GE.90.).AND.(T.LT.180. )) X— X + BB/ {2 .* AA ) 

SY 

139 

140 ^ 


IF ( (T.GE. 270.) .AND. (T.LE.360.)) X-X+BB/ < 2.+AA J 

SY 

140 



Y-X*TAN(T/RAD) 

SY 

141 



RETURN 

SY 

142 


23 

X-0.0 

SY 

143 



Y--A+SORT (( SCALE ( IS) *1.2222?) **?-A*A> 

SY 

144 

145' 


RETURN 

SY 

145 


24 

X — A 

SY 

146 



Y-X*TAN(T/RAD) 

SY 

147 



RETURN 

SY 

148 


25 

Y — A 

SY 

149 

150 


IF (T.EQ.270.) X.0.0 

SY 

150 



IF CT.NE.270.) X-Y/TAN (T/RAO ) 

SY 

151 



RETURN 

SY 

152 


C 


SY 

153 


C 

SYMBOL NUMB FR 7 OR 17 

SY 

154 

155 

C 


SY 

155 


26 

IF ( (T.GE.O. 0>. AND. (T.LT. T7)> GO TO 27 

SY 

156 



IF (( T.GE .T7>. AND. (T.LT. T8) ) GO TO 26 

SY 

157 



IF ( (T.GE. T8). AND. (T.LT. T9)) GO TO 29 

SY 

156 



IF ( (T.GE. T9). AND. (T.LT.S1) ) GO TQ 30 

SY 

159 

160 i 

27 

X«SCALE(I5)/2. 

SY 

160 
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LISTING OF DECK: SYMBS 
CARD NO, 


161 

165. 

170 

175 

180 i 

185 

190 

195 

200 


28 


29 

30 


C 

C 

C 

31 


32 


33 


34 


C 

C 

C 

35 


36 • 


Y-X*TAN(T/RAD) 

RETURN 

X*SCALE(IS)*.1*SIN(T/RAD)*C0S(T/RAD)+SCALE( IS) *COS( T/RAD ) *SORT( • 25 
l-(#l*CnS(T/RAD) ) **2 ) 

Y- X*T AN ( T /RAD ) 

RETURN 

X— SCALE(IS)/2. 

Y*X*TAN( T/RAD) 

RETURN 

Y— SCALE( I$)*.4 

IF (T.EQ.270.) X-0.0 

IF (T.NE.270.) X»Y/TAN (T/RAD) 

RETURN 

SYMBOL NUMBER 8 OR IB 

IF ( { T.GE'O'O) .ANO.C T,LT'S2) ) GO TO 32 
IF ( (T.GE.S2).AND. (T.LE.270.) ) GO TP 33 
IF ( (T.GT.270. ).AND. (T.LE.S3)) GO TO 34 

X— .05*SCALE( IS) *SIN(T/RAD)*COS (T/RAD )+SCALE( IS) *COS (T/RAD )*SORT(. 
125-(.05*CQS(T/RAD))**2) 

Y" X*TAN ( T /RAD) 

RETURN 

IF (T.EQ'270.) X-0,0 

IF (T aEQ'270' ) Y— SC ALE(IS)*,55 

IF ( T*NE • 270» ) X— SC AL E ( I S )* .55/ (T AN ( T/RAD )+l . ) 

IF (T.NE.270.) Y-X^T AN (T/RAD) 

RETURN 

X— SC ALE (IS) *.55/ (TAN (T/RAD )-l.) 

Y*X*TAN (T/RAD ) 

RETURN 

SYMBOL NUMBER 9 nR 19 

IF ( (T.GE.O.O) .AND. (T.LT.90. ) ) GO TO 36 
IF ( (T.GE.90.).AND.(T.LT.180«)) GO TO 37 
IF ((T.GE.180.J.AND. (T.LT.270.)) GO TO 38 
IF ( (T.GE.270.). AND. (T.LE.360.)) GO TO 39 
X«SCALE(IS)*R/2./(TAN(T/PAD)+R) 

Y« X*T AN ( T/RAD) 


SY 161 
SY 162 
SY 163 
SY 164 
SY 165 
SY 166 
SY 167 
SY 168 
SY 169 
SY 170 
SY 171 
SY 172 
SY 173 
SY 174 
SY 175 
SY 176 
SY 177 
SY 178 
SY 179 
SY 180 
SY 181 
SY 182 
SY 183 
SY 184 
SY 185 
SY 186 
SY 187 
SY 188 
SY 189 
SY 190 
SY 191 
SY 192 
SY 193 
SY 194 
SY 195 
SY 196 
SY 197 
SY 198 
SY 199 
SY 200 
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LISTING OF DECKt SYMBS 


PAGE 6 


CARO NO« 


201 

205 

210 

215 

220 

225 

230 

235 

2A0 


RETURN 

37 IF (T.EQ.90.) **0.0 

IF (T.NE.90.) X«SCALE(IS)*R/2./(TAN(T/RAD)-R) 

IF (T.EQ.90.) Y-SCALE (IS )*R/2. 

IF (T.NE.90.) Y«X*TAN (T/RAD ) 

RETURN 

3 8 X— SC ALE (IS )*R /2 . / ( TAN ( T/RAD ) +R ) 

y«X*TAN(T/RAD) 

RETURN 

39* IF (T.EQ.270.) X-0.0 

IF (T.EG.270.) Y— SCALF(lS)*R/2. 

IF (T. HE. 270.) X— SC ALE(IS)«<R/2./ (TAN (T/RAD)-R ) 
IF (T.NE.270.) Y-X*TAN (T/RAO) 

RETURN 

C 

C SYMBOL NUMBER 10 OP 20 

C 

AO IF { (T.GE.O.O) .AND. (T.LT.S4) ) GO TO *1 
IF ( (T.GE.SA).AND. (T.LT.S5) ) GO TO A2 
IF ( (T.GE.S5) .AND.(T .LT.S6) f GO TO A3 
IF ( (T.GE.S6). ANO.CT.LT. S7) ) GO TO AA 
A1 X-SCALE(IS)/2. 

Y»X*TAN(T/RAD) 

RETURN 

A 2 IF (T.LT.90.) X-SC AL E ( IS ) *.6/ (TAN ( T/R AD ) +1. ) 

IF (T.EQ.90.) X-0.0 

IF (T.GT.90.) X»SCALF( IS ) *. 6/ (TAN ( T/R AD ) -1 . > 

IF (T.EQ.90.) Y*SCALE( I5)*.6 
IF (T.NE.90.) Y*X*TAN(T/RAD) 

RETURN 

A3 X— SCALE(IS)/2. 

Y«X + T AN(T /RAD) 

RETURN 

AA Y*-SC ALE ( IS )*• A 

IF (T.EQ.270.) X-0.0 

IF (T. ME. 270.) X-Y/TAN (T/RAD) 

RETURN 

C . 

C SYMBOL NUMBER 21 OR 22 

C 


SY 201 
SY 202 
SY 203 
SY 20A 
SY 205 
SY 206 
SY 207 
SY 206 
SY 209 
SY 210 
SY 211 
SY 212 
SY 213 
SY 21A 
SY 215 
SY 216 
SY 217 
SY 218 
SY 219 
SY 220 
SY 221 
SY 222 
SY 223 
SY 22A 
SY 225 
SY 226 
SY 227 
SY 228 
SY 229 
SY 230 
SY 231 
SY 232 
SY 233 
SY 23A 
SY 235 
SY 236 
SY 237 
SY 238 
SY 239 
SY 2A0 
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LISTING OF DECKt SYMBS 
CARD NO. 

241 45 X-.1+SCALE(IS)*COS(T/RAO> 

Y-.1*SCALE(IS)*SIN(T/RA0) 

RETURN 

END 


PAGE 7 


SY 241 
SY 242 
SY 243 
SY 244 
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LISTING OF DECK* CUBSPL 


PAGE 


1 


CARD NO. 


1 


SUBROUTINE CUBSPL (X# Y, N# Y PP# TENS I ON# TENS# A) 

CB 

1 


C 


CB 

2 


C 

THIS SUBROUTINE FITS A CUBIC SPLINE TO A SET OF Y VS X INPUT 

CB 

3 


c 

POINTS 

C 8 

4 

5 

c 


CB 

5 


c 

CODED BY — HARRY L. MORGAN NASA/L ARC/TAD/AAB 1983 

CB 

6 


c 


CB 

7 


c 

IN CALLING PROGRAM DIMENSION X# Y# YPP# AND A BY N 

CB 

8 


c 


CB 

9 

10 


DIMENSION X (N )# Y(N), YPP(N)# A (N) 

CB 

10 


c 


CB 

11 


c 

COMPUTE TENSION PARAMETER 

CB 

12 


c 


CB 

13 



TENS*TENSION*FLOAT(N-l)/(X(N)-X(l)) 

CB 

14 

15 


IF (TENS. LT. 0.0) TENS-0.0 

CB 

15 


c 


CB 

16 


c 

COMPUTE SECOND DERIVATIVE AT END POINTS BY FITTING 

CB 

17 


c 

Y-A*X**2+B*X+C TO THE LAST THREE POINTS AND SOLVE FOR A. 

CB 

18 


c 

SECOND DERIVATIVE AT END POINT IS THEN EQUAL TO 2.*A 

CB 

19 

20 

c . 


CB 

20 



H1«X (21 —X 

CB 

21 



H2-X(3)-X(l> 

CB 

22 



H3«X( 1)-X(2> 

CB 

23 



YPP (1 )-2.*(Y(l)*Hl+Y(2)*H2+Y(3)*H3)/(Hl*X(l)**2+H2*X(2) **2+H3*X(3 ) 

CB 

24 

25 


1**2) 

CB 

25 



H1-X(N-1)-X(N) 

CB 

26 



H2«X(N)-X(N-2) 

CB 

27 



H3-X( N-2 ) -X (N-l ) 

CB 

28 



YPP(N)-2.*( Y(N-2)*H1 +Y (N-l )*H2+Y ( N ) *H3 ) / ( H1*X ( N-2 ) **2+H2*X (N-l )**2 

CB 

29 

30 . 


1*H3*X (N )**2 ) 

CB 

30 


c 


CB 

31 


c 

PERFORM FORWARD ELIMINATION 

CB 

32 


c 


CB 

33 



Nl-N-1 

CB 

34 

35 


Hl-X(2)-X(l) 

CB 

35 



IF (TENS. NE. 0.0) GO TO 1 

CB 

36 



C 0N«6 . 

CB 

37 



E«H1 

CB 

38 



GO TO 2 

CB 

39 

40 

1 

CON«TENS*TENS 

CB 

40 
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LISTING OF DECK! CUBSPL 


PAGE 


2 


CARD NO. 


41 


45 


50 


55 


60 


65 


70 


H3-TENS/SINH(TENS*H1) 

CTH-H3*CaSHtTENS*Hl) 

E-1./H1-H3 

0P-CTH-1./H1 

2 A ( 1) "0*0 

DO 5 I-2»N1 
H2-X(I+1)-X(I) 

C-CON*( <Y(H-1)-Y(I ))/H2-(Y(I)-Y(I-l) ) / HI > 

IF (TENS. NE. 0.0) GO TO 3 

F-H2 

D-H1*(2.-A(I-1)>+2.*H2 
GO TO 4 

3 H3-TENS/SINH(TfcNS*H2) 

CTH-H3*CQSH(TENS*H?) 

F-1./H2-H3 

DN-CTH-1./H2 

D*DN+DP“E*A(I-1) 

DP-ON 

4 A(I)-F/D 
YPP(I)-(C-E*YPP(I-1) )/0 
E-F 

5 H1-H2 
C 

C PERFORM BACK SUBSTITUTION 

C 

J-N 

DO 6 I-2#N1 
J-J-l 

6 YPP(J )-YPP( J)-A( J)*YPP (J+l) 

C 

C RETURN TO CALLING PROGRAM 

C 

RETURN 

END 


CB 

41 

CB 

42 

CB 

43 

CB 

44 

CB 

45 

CB 

46 

CB 

47 

CB 

4P 

CB 

49 

C8 

50 

CB 

51 

CB 

52 

CB 

53 

CB 

54 

CB 

55 

CB 

56 

CB 

57 

CB 

58 

CB 

59 

CB 

60 

CB 

61 

CB 

62 

CB 

63 

CB 

64 

CB 

65 

CB 

66 

CB 

67 

CB 

68 

CB 

69 

CB 

70 

CB 

71 

CB 

72 

CB 

73 

CB 

74 
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PAGE 


1 


CARD NO. 
1 


5 


1Q- 


15 


FUNCTION FUNC (X,X1, X2, Yl# Y2, YPP1> YPP2,TENS) FC 1 

C FC 2 

C ROUTINE TO COMPUTE Y-VALUF AT A GIVEN X-VALUE ALONG SPLINE CURVE FC 3 

C FC 4 

C CODED BY — HARRY L. MORGAN NASA/LARC/TAD/ AAB 1983 FC 5 

C FC 6 

DX 1*X 2-X FC 7 

DX 2*X-X1 PC 8 

DH* X2-X1 PC 9 

IF (TENS. EQ. O.O) GO TO 1 PC 10 

F1-YPP1/(TENS*TENS) PC 11 

F2*YPP2^( TENS* TENS) PC 12 

FUNC-(Fl*$INH(TENS*DXl)+F2*SINH(TENS*DX2))/SINHCTENS*DH)+( CY1-F1)* FC 13 

ID XI* ( Y2-F2 )*DX 2 ) /DM . FC 14 

RETURN FC 15 

1 FUNC«YPP1*DX1**3/(6.*DH)*YPP2*DX2**3/(6.*DH)+(Y1-YPP1*DH*DH/6.)*DX FC 16 

ll/DH*(Y2-YPP2*DH*0H/6. )+DX2fDH FC 17 

RETURN PC 18 

END PC 1 Q 
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PAGE 1 


CARO NO, 


1 


SUBROUTINE DASHLN (T, X,Y, N# NSYM# IOP# ISIZE# IRUN, TENSION) 

DH 

1 


C 


DH 

2 


c 

ROUTINE TO FAIR DATA WITH DASH LINES 

DH 

3 


c . 


DH 

4 

5 

c 

CODED BY — HARRY L, MORGAN NASA/L ARC/TAD/ AAB 1983 

DH 

5 


c 


DH 

6 


c 

PARAMETER DEFINITION 

DH 

7 


c 

T - ARRAY CONTAINING INDEPENDENT VARIABLE 

DH 

8 


c 

X - ARRAY CONTAINING X VARIABLE 

DH 

9 

10 

c 

Y - ARRAY CONTAINING Y VARIABLE 

DH 

10 


c 

N - NUMBER OF POINTS TO BE PLOTTED 

DH 

11 


c 

NS YM - LINE OR SYMBOL NUMBER 

DH 

12 


c 

IOP - PLOTTING OPTION I0P-0 PLOT SYMBOLS ONLY 

OH 

13 


c 

IOP-1 PLOT DASH LINE 

DH 

14 

15 

c 

ISIZE - SYMBOL SIZE 

DH 

15 


c 

IRUN - RUN NUMBER 

DH 

16 


c 

TENSION - SPLINE TENSION 

DH 

17 


c 


DH 

18 


c 

T# X# AND Y MUST BE PROPERLY DIMENSIONED IN CALLING ROUTINE 

DH 

19 

20 


DIMENSION T ( 1) # X(l), Y ( 1 )# NL(22)# NS<22) 

DH 

20 


c 

WORK ARRAYS 

DH 

21 



COMMON /WORK/ $ (166) #X I ( 166) # YI ( 166) 

DH 

22 



COMMON / PLT/ MX (50 )* MY ( 50)# A ( 50) 

DH 

23 


c 


DH 

24 

25 


REAL MX# MY 

DH 

25 


c 

ROUTINE CONSTANTS 

DH 

26 


c 

NL ■ NUMBER OF LONG DASHES AND NS • NUMBER OF SHORT DASHES 

DH 

27 



DATA NL/1#1#1»1#1#2#2#2#2#3#3»3#3#4#4»4#4#5#5#5#5#5/ 

DH 

28 



DATA NS /O# 1# 1#2# 3# 1# 2#3#4#1#2#3#4#1#2#3#4#1# 2# 3#4#5/ 

DH 

29 

30 ' 

c 


DH 

30 


c 

ROUTINE CONSTANTS 

DH 

31 


c 


DH 

32 



NP-166 

DH 

33 



EPS-,02 

OH 

34 

35’ 


SL-.3 

DH 

35 



IF ( NS YM • EO ,2 ) SL-.2 

DH 

36 



SS-.l 

DH 

37 



SP-,1 

DH 

38 


c . 


DH 

39 

40 

c 

SELECT PLOTTING OPTION 

DH 

40 
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CARD NO. 


41 

C 


DM 

41 



IF (I0P.EQ.0.0R.N.LT.3 ) GO TO 1 

DM 

42 



GO TO 3 

DM 

43 


C 


DH 

44 

45 

c 

PLOT SYNBOLS ONLY 

DH 

45 


c 


DH 

46 


1 

DO 2 1*1, N 

DH 

47 


2 

CALL PNTPLT ( XC I } # YC I > ,NSYM, I SIZE ) 

DH 

48 



RETURN 

DH 

49 

50. 

C 


DH 

50 


C 

PLOT DASH LINE 

DH 

51 


C 


DH 

52 


C 

CHECK TO SEE IF T IS STRICTLY INCREASING 

DH 

53 


3 

DO 4 1*2# N 

DH 

54 

55 


IF (T ( I ) • LT • T( 1-1 ) ) GO TO 5 

DH 

55 


4 

CONTINUE 

DH 

56 



GO TO 6 

DH 

57 


C 

PRINT ERROR MESSAGF IF T IS NOT STRICTLY INCREASING 

DH 

58 


5 

PRINT 15# IRUN 

DH 

59 

60 


PRINT 16, (T(I),X(I),Y(I),I*1,N) ? 

DH 

60 


C . 

IF T IS NOT STRICTLY INCREASING, PLOT SYMBOLS ONLY 

DH 

61 



GO TO 1 

DH 

62 


c 

FIT SPLINE CURVE THROUGH DATA POINTS 

DH 

63 


6 

CALL CUBS PL (T, X, N,M X, TENSION, TENS, A ) 

DH 

64 

65 


CALL CUBSPL (T,Y, N, MY, TENSION, TENS, A > 

DH 

65 


C 

COMPUTE NP EQUALLY SPACED POINTS ALONG CURVE 

DH 

66 



IF (NSYM.EQ.l) NP-NP-2 

DH 

67 



DT*<T(N)-T(1))/FL0AT(NP-1) 

DH 

68 



TI-T(l) 

DH 

69 

70 . 


S(l)-0. 

DH 

70 



M*2 

DH 

71 



DO <5 1*1, NP 

OH 

72 



IF (TI.LT.T(l)) TI-T(l) 

DH 

73 



IF (TI.GT.T(N)) TI*T (N ) 

DH 

74 

75] 


DO 7 J*M, N 

DH 

75 



K-J-l 

DH 

76 



IF (TI.GE.TtKJ.AND.TI.LE.TtJ)) GO TO 8 

DH 

77 


7 

CONTINUE 

DH 

78 


8 

M-K 

DH 

79 

80. 


IF (M.LE.l) M*2 

DH 

80 


84 



3 


LISTING OF DECK* DASHLN 


PAGE 


CARO NO. 


81 

85 

90 

95 

100 

105 

110 

115 

120 


XI (I )«FUNC <TI#T(K ).T(K+1)#X(K)#X(K+1)#MX(K)#MX(K+1)#TFNS) 
YI(I)*FUNC(TI#T(K)#T(K+1)#Y(K)#Y(K+1)#MY(K)#MY(K+1)#TENS) 

IF (I .EQ.l.OR.NSYM.FQ.l) GO TO 9 

sm-sti-n+soRTuxi m-xiu-i) j**2+(yk imhm))**2) 

9 TI-TI+DT 
ST-S(NP) 

C IF NS YM-1# PLOT CONTINUOUS LINE 

IF (NSYM.EQ.l) GO TO 13 

C COMPUTE LENGTH OF LONG AND SHORT DASH AND SPACE BETWEEN EACH 

NLG-NL(NSYM) 

NSG-NS(NSYM) 

DT-F LOAT ( NLG) 4SL+FL0 AT (NSG )*SS+FLO AT( NLG+NSG) *SP 
NLT-IFIX< (ST+SP) /DT) 

IF (NLT.LE.O) NLT-1 
QTN-CST+SP) /FLOAT(NLT) 

TI-DTN/DT 

SL«SL*TI 

SS-S$*TI 

SP-SP*TI 

C POSITION PFN AT START OF LINE . 

SI-0. 

CALL CALPLT ( XT ( 1 ) # YI (1 ) #3 ) 

JSTART-1 

C PLOT NL LONG DASHES 

10 DO 11 I-lfNLG 
SI-SI+SL 

IF (SI+EPS.GE.ST) SI -ST 

CALL LINEAR (SI# XX# YY»NP# S#XI#YI# JSTART) 

CALL CALPLT (XX,YY#2) 

SI-SI4SP 

IF (SI.GE.ST) GO TO 14 

CALL LINEAR ( SI, XX, YY, NP, S# XI# YI, JSTART) 

CALL CALPLT (XX#YY#3) 

11 CONTINUE 

C PLOT NS SHORT DASHES 

IF (NSG.EQ.O) GO TO 10 

DO 12 I-l# NSG 

SI-SI4SS 

IF (SI+EPS.GE.ST) SI-ST 

CALL LINEAR ( ST# XX# YY, NP# S# XI » YI # JSTART) 


DH 

81 

DH 

82 

DH 

83 

DH 

84 

DH 

85 

DH 

86 

DH 

87 

DH 

88 

DH 

89 

DH 

90 

DH 

91 

DH 

92 

DH 

93 

DH 

9 A 

DH 

95 

DH 

96 

DH 

97 

DH 

98 

DH 

99 

DH 

100 

DH 

101 

DH 

102 

DH 

103 

DH 

104 

DH 

105 

DH 

106 

DH 

107 

DH 

108 

DH 

109 

DH 

110 

DH 

111 

DH 

112 

DH 

113 

DH 

114 

DH 

115 

DH 

116 

DH 

117 

DH 

118 

DH 

119 

DH 

120 
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121 


CALL CALPLT <XX,YY,2) 

DH 

121 



SI-SI+SP 

DH 

122 



IF (SI.GE.ST) GO TO 14 

DH 

123 



CALL LINEAR (S I# XX, YY, NP# S, XI. YI* JSTA RT) 

DH 

124 

125 


CALL CALPLT (XX*YY#3) 

DH 

125 


12 

CONTINUE 

DH 

126 


C 

CONTINUE PLOTTING DASH LINE 

DH 

127 



GO TO 10 

DH 

128 


C 

PLOT CONTINUOUS LINE 

DH 

129 

130 

13 

XI(NP4l)-0.0 

OH 

130 



YI (NP + 1)«0.0 

DH 

131 



XI (NP42)«YI(N»4-2)*1.0 

OH 

132 



CALL LINE (XI#YI*NP# 1#0*0*0.) 

DH 

133 


C 

* 

DH 

134 

135 

C 

RETURN TO START OF CURVE 

DH 

135 


c 


DH 

136 


14 

CALL CALPLT ( XI ( 1 )» YI ( 1 ), 3 ) 

DH 

137 


c 


DH 

138 


c 

RETURN TO CALLING PROGRAM 

DH 

139 

140 

c 


DH 

140 



RETURN 

DH 

141 


c 


DH 

142 


15 

FORMAT (1H1//5X# 36HT IS NOT STRICTLY INCREASING FOR 

RUN# 17/1 6X# 1HT DH 

143 



1*13X*1HX#13X#1HY/J 

DH 

144 

145 

16 

FORMAT (5X# 3F15. 4) 

DH 

145 



END 

DH 

146- 
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CARD NO. 





L 


SUBROUTINE LINEAR (T TN T , X INT, Y INT , N, T, x» Y, J STA RT ) 

LR 

1 


C 


LR 

2 


C 

LINEAR INTERPOLATION MPTHHD 

LR 

3 


C 


LR 

4 

5 

c 

CODED BY — HARRY L* MORGAN NAS A/ L ARC / TAD / AAB 

1983 LR 

5 


c 


LR 

6 



DIMENSION T ( 1 ) # X(l), Y ( 1 ) ' 

LR 

7 



J"1 

LR 

8 



IE (TINT.LE.T(l)) GO TO 3 

LR 

9 

10 


J ■ N 

LR 

10 



IF (TINT.GE.T(N) ) GO TO 3 

LR 

11 



JEND* JSTART 

LR 

12 



DO 1 I«JSTART»N 

LR 

13 



J"I 

LR 

14 

15. 


IF (T(J)-TINT) 1,3,2 

LR 

15 


1 

JEND-J 

LR 

16 


2 

DT»(T INT-T ( J-1) ) / ( Tf J ) -T{ J -1 ) ) 

LR 

17 



XINT-Y( J-1) + (X(J 1 -X f J-1) )*DT 

LR 

18 



YINT-Yt J-1 > ♦ ( V ( J )-Y( J-l) )*0T 

LR 

19 

zo. 


JSTART-JEND 

LR 

20 



RETURN 

LR 

21 


3 

YINT-Y(J) 

LR 

22 



XINT«X(J ) 

LR 

23 



J START* J 

LR 

24 

25 


RETURN 

LR 

25 



END 

LR 

26' 
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PAGE 

1 

D NO, 





1 


SUBPOUT IN E LSOPLT (T ,X , Y,N>N SYM, IS I Z E, IOP> I R UN ,TE*JS I ON, U SO#NDOL 

>D LS 

I 



IF) 

LS 

2 


C 


LS 

3 


C 

ROUTINE TO FAIR DATA WITH LF AST SQUARFS CURVE FIT 

LS 

4 

5 

C 


LS 

5 


C 

CODED BY — HARRY L. MORGAN NASA/LARC/TAD/ AAB 1983 

LS 

6 


c 


LS 

7 


c 

PARAMETER DEFINITION 

LS 

8 


c 

T - ARRAY CONTAINING INDEPENDENT VARIABLE 

LS 

9 

10 

c 

X - ARRAY CONTAINING X VARIABLE 

LS 

10 


c 

Y - ARRAY CONTAINING Y VARIABLE 

LS 

11 


c 

N - NUMBER 0* POINTS TO BE PLOTTFD 

LS 

12 


c 

NS YM - LINF AND SYMBOL NUMBER 

LS 

13 


c 

ISIZE - SYMBOL SIZE 

LS 

14 

15 

c 

IOP - PLOTTING OPTION I0P-0 PLOT SYMBOLS ONLY 

LS 

15 


c 

IOP-1 PLOT SYMBOLS AND CURVr 

LS 

16 


c 

IRUN - RUN NUMBER 

LS 

17 


c 

TENSION - SPLINE TENSION 

LS 

18 


c 

ILSO - LEAST SQUARES CURVE PIT OPTION 

LS 

19 

20 

c 

ILSO-O LEAST SQUARES POLYNOMIAL CURVE FIT 

LS 

20 


c 

ILSQ-1 LEAST SQUARES CUBIC SPLINE FIT 

LS 

21 


c 

NPOL - ORDER OF LEAST SQUARES POLYNOMIAL CURVE FIT 

LS 

22 


c 

DF - STANDARD DEVIATION FOR LEAST-SQUARES CUBIC-SPLINE FIT 

LS 

23 


c 

DIMENSION T, X# AND Y IN CALLING PROGRAM 

LS 

24 

25 


DIMENSION TCI), X(l) , Y(l> 

LS 

25 


c 

WORK ARRAYS 

LS 

26 



DIMENSION XNEW( 50 ), YNEW(50) 

LS 

27 



DIMENSION COEF ( 50# 4) 

LS 

28 



COMMON /WORK/ WK(360),SDV(50)#CX(11),CYC11) 

LS 

29 

30 

c 


LS 

30 


c 

ROUTINE CONSTANTS 

LS 

31 


c 


LS 

32 



NMAX-50 

LS 

33 



WT-1.0 

LS 

34 

35 


I PRINT*© 

LS 

35 


c 


LS 

36 


c 

SELECT PLOTTING OPTION 

LS 

37 


c 


LS 

3 6 



IF (IOP. EQ.0.0R.N.LT.3) GO TO 1 

LS 

39 

40 


GO TO 3 

LS 

40 
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PAGE 

2 

C 


LS 

41 

C 

PLOT SYMBOLS (INLY 

LS 

42 

C 


LS 

43 

1 

DO 2 I "1 * N 

LS 

44 

2 

CALL PNTPLT ( X (I) , Y ( I ) #N S YM # I S IZ E ) 

LS 

45 


RETURN 

LS 

46 

C 


LS 

47 

C 

CHECK TO SEE IF T IS STRICTLY INCREASING 

LS 

48 

c 


LS 

49 

3 

DO 4 1*2# N 

LS 

50 


IF (T(I).LT.T(T-l) ) GO TO 5 

LS 

51 

4 

CONTINUE 

LS 

52 


GO TO 6 

LS 

53 

5 

PRINT 20; IRUN 

LS 

54 


PRINT 21# (T(I)#X(I)#Y(I),I-1,N) 

LS 

55 


GO TO 1 

LS 

56 

C 


LS 

57 

C 

PLOT DATA POINTS WITH SYMBOLS 

LS 

50 

c 


LS 

59 

6 

DO 7 I* 1# N 

LS 

60 


SDV( I ) ■ WT 

LS 

61 


IF (ILSQ.EQ.l) SDV(I)»DF 

LS 

62 


CALL PNTPLT C X ( I ) # Y ( I ) #N S YM# I S I ZE ) 

LS 

63 

7 

CONTINUE 

LS 

64 

C 


LS 

65 

C 

COMPUTE LEAST SQUARES CURVE 

LS 

66 

c 


LS 

67 


IF (ILSQ.GT.0) GO TO 10 

LS 

68 

c 

COMPUTE COEFFICIENTS OF LEAST SQUARES POLYNOMIAL 

LS 

69 


MPOL- NPOl 

LS 

70 


IF (MPOL.GT.N-1) MPOL-N-1 

LS 

71 


IF (MPOL.LE.O) MPOL* 1 

LS 

72 


IF (MPOL •GT*IO) MPQL *10 

LS 

73 


Ml-KPOL+1 

LS 

74 


CALL LSQ (T#X#SDV»N#MP0L#CX) 

LS 

75 


CALL LSO (T#Y#SDV#N#MPOL#CY) 

LS 

76 

c 

COMPUTE NEW X AND Y VALUES AND ERRCR SUMMATIONS 

LS 

77 


EPRX-ERRY*0.0 

LS 

78 


DO 9 I* 1 # N 

LS 

79 


XNEWU J-CX(l) 

LS 

80 
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8 1 
85 ! 

90 

95 

100 

105 ^ 

110 

115 

120 


.8 


9 

C 

C 

10 


C 


11 


12 

13 

C 


C 


14 


YNEW( I ) -C Y ( 1 } 

DELTA-1.0 
DO 8 J-2#M1 
DELTA -DELTA* T( I) 

XNEW( I)-XNEW(I)+DELTA*CX( J) 

YNEW ( I )» YNE W ( I )*D C LTA*C Y ( J ) 

ERRX-ERRX*(X(I>-XNEV(I ))**2 
ERRY-ERRY-MYd )-YNEW(I ))**2 ' 

CONTINUE 
GO TO 17 

COMPUTE COEFFICIENTS OF LEAST SQUARES CUBIC SPLINE 
FIT X 

K- — 1 

SM-FLOAT (N ) 

CALL CSDS (NMAX>N»T*X*SDVfSM>K»COEF»WK>IERR) 

IF { I ERR ,NE .0 ) GO TO 19 

COMPUTE NEW X AND ERROR SUMMATION ' 

ERRX-0 • 

DO 13 I-1,N 

IF ( I .EQ.N) GO TO 11 

XNEW ( D-COEF C I# 1 ) 

GO TO 12 
K-I-l 

DELTA -T ( I )-T(K ) 

XN EW ( I )- ( (C0EF(K> 4)*DELTA+C0EF(K>3 ) )*OELTA+COEF(K>2 ) )*DELTA+C0EF (K 
1#1> 

ERRX-£RRX+(X(I)-XNEW(I ) )**£ 

CONTINUE 

FIT Y 

K* -1 

SH-FLOAT(N) 

CALL CSDS <NMAX»N,T,Y,SDV>SM#K»COEF#WK>IERR) 

IF ( I ERR ,NE .0) GO TO 19 

COMPUTE NEW Y AND ERROR SUMMATION ' 

ERRY-O. 

DO 16 I-1>N 

IF {I. EQ.N) GO TO 14 

YNEW( I)-C0EF(T#1) 

GO TO 15 
K-I-l 


LS 

81 

LS 

82 

LS 

83 

LS 

84 

LS 

85 

LS 

86 

LS 

87 

LS 

88 

LS 

89 

LS 

90 

LS 

91 

LS 

92 

LS 

93 

LS 

94 

LS 

95 

LS 

96 

LS 

97 

LS 

98 

LS 

99 

LS 

100 

LS 

101 

LS 

102 

LS 

103 

LS 

104 

LS 

105 

LS 

106 

LS 

107 

LS 

108 

LS 

109 

LS 

110 

LS 

111 

LS 

112 

LS 

113 

LS 

114 

LS 

115 

LS 

116 

LS 

117 

LS 

118 

LS 

119 

LS 

120 
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121 

125 

130. 

135 

140 ^ 

145 

150 

155 

16 o! 


DELTA-T (I )-T (K ) LS 121 

YNEW(I)-( (C0EF(K,4)*DELTA*C0EF(K,3) ) *DELTA+COEF <K,2 ) )*DELTAKOFF (K LS 122 

1,1) LS 123 

15 ERRY-ERRY+(Y(I)-YNEW(I ))**2 LS 12A 

16 CONTINUE LS 125 

C LS 126 

C PLOT DASHED LINE THRU NEW X AND Y VALUES LS 127 

C LS 128 

17 TEN-TENSION LS 129 

IF (ILSQ.EQ.l) TEN-0.0 LS 130 

CALL DASHLN {T,XNEW> YNEV,N,NSYM,IOP, ISI7E, IRUN,TEN ) LS 131 

C LS 132 

C PRINT RESULTS OF LEAST SQUARES CURVE FIT LS 133 

C LS 13A 

IF (I PRINT .EQ.O) RETURN LS 135 

IF (ILSQ.EQ.O) PRINT 22, IRUN, NSYM, HPOL LS 136 

IF (ILSQ.EQ.l) PRINT 23, IRUN,NSYM»DF LS 137 

DO 18 1-1, N LS 138 

DELTA-X<I)-XNFV< I) LS 139 

SM-Yt I)-YNEW(I) LS 140 

PRINT 24, I,T(I)»X(I),XNEW(I)»DELTA,Y(I),YNEW(I),SM LS 141 

18 CONTINUE LS 142 

PRINT 25, ERR Y, ERR Y LS 143 

RETURN LS 144 

C LS 145 

C PRINT ERROR MESSAGE LS 146 

C LS 147 

19 PRINT 26, IERR, IRUN, NSYM LS 148 

GO TO 1 LS 149 

C LS 150 

20 FORMAT (1H1//5X, 36HT IS NOT STRICTLY INCREASING FOR RUN, I 7/16X, 1HT LS 151 

1,13X,1HX,13X,1HY/J LS 152 

21 FORMAT (5X,3F15.4) LS 153 

22 FORMAT { //5X, 55HR ESULTS OF LEAST SQUARES POLYNOMIAL CURVE PIT FOR LS 154 

IRUN -,I4,3X,6HNSYM -, 1 2, 3X, 6HNP0L -,I3//4X,1HI,6X,1HT,9X, 1HX,9X,4H LS 155 

2XNEW,4X, 6HX-XNEW,6X, 1HY,9X, 4HYNEW, 4X, 6HY-YNEW ) LS 156 

23 FORMAT ( //5X,5 7HRESULTS OF LEAST SQUARES CUBIC SPLINE CURVE FIT FO LS 157 

1R RUN -, I4,3X,6HNSYM -,I2,3X,4HDF -, F10.5//4X, 1HI, 6X, 1HT, 9X, 1HX, 9X LS 158 

2,4HXNEW, 4X,6HX-XNEW,6X»1HY,9X, 4HYNEW,4X, 6HY-YNEW ) LS 159 

24 FORMAT (I5,7F10«4) LS 160 


9 ! 
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161 


25 

26 


FORMAT ( /10X#12HFRRDR IN X 
FORMAT ( /5X#13HERR0R NUMBER 
IT FOR RUN ■# 14# 11H AND NSYM 
END 


#E15«6#10X#12HERR0R IN Y *#F15.6) LS 161 
# 13* 42 H OCCURRED IN CSDS CALL IN LSQPL LS 162 
* ,I3 > LS 163 

LS 164- 
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LISTING OF DECK t LSQ 


PAGE 


1 


CARO NO. 


1 


SUBROUTINE LSQ (X#Y#W#NP#N#C) 


LQ 

1 


C . 



LQ 

2 


c 

LEAST SQUARES WEIGHTED POLYNOMIAL CURVE FIT ROUTINE 


LQ 

3 


c 



LQ 

4 

5 

c 

CODED BY — HARRY L. MORGAN NAS A/L ARC /TAD/ AAB 

1983 

LQ 

5 


,c 



LQ 

6 


c 

PARAMETER DEFINITION 


LQ 

7 


c 

X - ARRAY CONTAINING INDEPENDENT VARIABLE 


LQ 

8 


c 

Y - ARRAY CONTAINING DEPENDENT VARIABLE 


LQ 

9 

10 

c 

W - ARRAY CONTAINING WEIGHTING VALUES 


LQ 

10 


c 

NP - NUMBER OF POINTS 


LQ 

11 


c 

N - ORDER OF POLYNOMIAL 


LQ 

12 


c 

C - COEFFICIENTS OF POLYNOMIAL 


LQ 

13 


c . 



LQ 

14 

15 

c 

NEW Y VALUE - C (1)+C (2>*X*C C3)*X**2+. . .♦C(N-H )*X**N 


LQ 

15 


c 



LQ 

16 


c 

DIMENSION X(NP), Y(NP)# W(NP), AND C(N*1) IN CALLING 

PROGRAM 

LQ 

17 



DIMENSION XU># YU)# W(l)» C(l) 


LQ 

18 


c 

DIMENSION A(N + 1#N*2) N-10 FOR THIS VERSION 


LQ 

19 

20 


COMMON /WORK/ A(U#12) 


LQ 

20 


c 



LQ 

21 


c 

COMPUTE LEAST SOUARES MATRIX 


LQ 

22 


c 



LQ 

23 



Nl-N+1 


LQ 

24 

25 


N2-N+2 


LQ 

25 



DO 1 I-1#N1 


LQ 

26 



DO 1 J«1*N2 


LQ 

27 


i ! 

A ( I> J )-0. 


LQ 

28 



DO 3 K«1#NP 


LQ 

29 

30 


Tl*l • 


LQ 

30 



DO 3 J«1#N1 


LQ 

31 



T2-T1 


LQ 

32 



DO 2 I-1#N1 


LQ 

33 



A( J#I ) -A ( Jil )+T2*W(K ) 


LQ 

34 

35 

2 

T2-T2*X(K) 


LQ 

35 



A ( J# N 2 ) "A ( J# N2)-Y(K)*T1*W(K> 


LQ 

36 


3 

T1 *T1*X (K) 


LQ 

37 


C 



LQ 

38 


c 

SOLVE FOR COEFFICIENTS OF LEAST SQUARES POLYNOMIAL 


LQ 

39 

40. 

c 



LQ 

40 
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LISTING OF DECKi L$Q 


CAPO NO, 

41 DO 4 K-1,N 

00 4 J*K>N 

Tl-A ( J4l»K) /A(K#K) 

DO 4 I *K> N2 

45 4 A(J+l*I)»A(J+i;i)-A(K# I)*T1 

C(N1)»-A(N1#N?)/A(N1#N1) 

DO 5 1-2, N1 
K-N2-I 

C(K)— A(K*N2)/A <K,K) 

50 L«K+1 

DO 5 J*L> N1 

5 C(K)»CCK)~C{J)*A(K#J)/A(KfK) 

RETURN 
END 
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LQ 41 
LO 42 
LO 43 
LO 44 
LQ 45 
LQ 46 
LQ 47 
LQ 48 
LQ 49 
LQ 50 
LQ 51 
LQ 52 
LQ 53 
LQ 54- 



LISTING OF DECK* CSDS 


PAGE 


1 


CARD NO. 
1 


5 


10 


15 


20 


25 


30 


35 


40 



SUBROUTINE 

CSOS ( M AX#IX#X*F#DF#S#IPT#CQEF#WK#IERR) 


CS 

1 

C **************************************** *********** ******************** 

cs 

2 

c* 



+ 

CS 

3 

c* 

PURPOSE* 


* 

cs 

4 

c* 


SUBROUTINE CSDS FITS A SMOOTH CUBIC SPLINE TO A 

* 

cs 

5 

c* 


UNIVARIATE FUNCTION. DATA MAY BE UNEQUALLY SPACED. 

* 

cs 

6 

c* 



♦ 

cs 

7 

c* 

USE* 


♦ 

cs 

8 

c* 


CALL CSDS (MAX, IX#X#F,DF,S*IPT#COEF#WK, IERR ) 

* 

CS 

9 

c* 



* 

cs 

10 

c* 

MAX 

INPUT INTEGER SPECIFYING THE MAXIMUM NUMBER OF DATA 

* 

cs 

11 

c* 


POINTS FOR THE INDEPENDENT VARIABLE. 

* 

cs 

12 

c* 



* 

cs 

13 

c* 

IX 

INPUT INTEGER SPECIFYING THE ACTUAL NUMBER OF DATA 

* 

cs 

14 

c* 


POINTS FOR THE INDEPENDENT VARIABLE. IX3MAX. 

* 

cs 

15 

c* 



* 

cs 

16 

c* 

X 

ONE-DIMENSIONAL INPUT ARRAY DIMENSIONED AT LEAST 

* 

cs 

17 

c* 


IX IN THE CALLING PROGRAM. UPON ENTRY TO CSDS* 

♦ 

cs 

18 

c* 


X(I) MUST CONTAIN THE VALUE OF THE. INDEPENDENT 

* 

cs 

19 

c* 


VARIABLE AT POINT I. 

* 

cs 

20 

c* 



* 

cs 

21 

c* 

. F 

ONE-DIMENSIONAL INPUT ARRAY DIMENSIONED AT LEAST 

* 

cs 

22 

c* 


IX IN THE CALLING PROGRAM. UPON ENTRY TO CSDS# 

* 

cs 

23 

c* 


F(I) MUST CONTAIN THE VALUE OF THE FUNCTION AT 

* 

cs 

24 

C* . 


POINT X(I>. 

* 

cs 

25 

c* 



* 

cs 

26 

c* 

DF 

ONE-DIMENSIONAL INPUT ARRAY DIMENSIONED AT LEAST 

* 

cs 

27 

c* 


IX IN THE CALLING PROGRAM. UPON ENTRY TO CSDS# 

* 

cs 

28 

c* 


DF ( I ) MUST CONTAIN AN ESTIMATE OF THF STANDARD 

* 

cs 

29 

c* 


DEVIATION OF F(I>. 

* 

cs 

30 

c* 



* 

cs 

31 

c* 

S 

A NON-NEGATIVE INPUT PARAMETER WHICH CONTROLS THE 

* 

cs 

32 

c* 


EXTENT OF SMOOTHING. S SHOULD BE IN THE RANGE 

* 

cs 

33 

c* 


( IX-(2*TX)**,5)<S<(IX*( 2* IX) *+.5) . 

* 

cs 

34 

c* 



* 

cs 

35 

c* 

IPT 

INPUT INITIALIZATION PARAMETER. THE USER MUST 

* 

cs 

36 

c* 


SPECIFY I P T» — 1 WHENEVER A NEW X ARRAY IS 

* 

cs 

37 

c* 


INPUT. THE ROUTINE WILL ALSO CHECK TO INSURE THAT 

* 

cs 

38 

c* 


THE X ARRAY IS IN STRICTLY INCREASING ORDER. 

* 

cs 

39 

c* 



* 

cs 

40 
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LISTING OF DECK* CSD$ 


PAGE 2 


CARD NO. 


41 

c* 

C OFF 

A TWO-DIMENSIONAL OUTPUT ARRAY DIMENSIONFD (MAX, 4) 

* 

cs 

41 


c* 


IN THE CALLING PROGRAM. UPON RETURN, COEF(I,J) 

* 

cs 

42 


c* 


CONTAINS THE J-TH COEFFICIENT OF THE SPLINE FOR 

* 

cs 

43 


c* 


THE INTERVAL BEGINNING AT POINT X(I). THE 

* 

cs 

44 

45 

c* 


FUNCTIONAL VALUE OF THE SPLINE AT ABSCISSA XI, 

* 

cs 

45 


c* 


WHERE X < I )<X1<X(H*1>, IS GIVEN BY* 

* 

cs 

46 


c* 


F(X1 )*((C0EF(I,4)*H+C0EF(I#3))*H+C0EF(I,2))*H 

* 

cs 

47 


c* 


♦CQEF ( 1,1) 

* 

cs 

48 


c* 


WHERE H-Xl-X(I) 

* 

cs 

49 

50 

c* 



* 

cs 

50 


c* 

WK 

A ONE-DIMENSIONAL WORK AREA ARRAY DIMENSIONED AT 

* 

cs 

51 


c* 


LEAST (7*IX+9) IN the CALLING PROGRAM. 

* 

cs 

52 


c* 



* 

cs 

53 


c* 

I ERR 

OUTPUT ERROR PARAMETER* 

* 

cs 

54 

55 

c* 


*0 NORMAL RETURN. NO ERROR D ETECTFD. 

* 

cs 

55 


c* 


•J THE J-TH ELEMENT OF THE X ARRAY IS NOT IN 

* 

cs 

56 


c* 


STRICTLY INCREASING ORDER. 

* 

cs 

57 


c* 


—1 THERE ARE LESS THAN FOUR VALUES IN THE X ARRAY 


cs 

58 


c* 



♦ 

cs 

59 

60 

c* 


UPON RETURN FROM CSDS, THIS PARAMETER SHOULD BE 

* 

cs 

60 


c* 


TESTED IN THE CALLING PROGRAM. 

* 

cs 

61 


c* 



* 

cs 

62 


c* 



♦ 

cs 

63 


c* 



* 

cs 

64 

65 

c* 

REQUIRED 

ROUTINES -NONE 

♦ 

cs 

65 


c* 



* 

cs 

66 


c* 

SOURCE 

IMSL ROUTINE ICSSMU MODIFIED BY 

* 

cs 

67 


c* , 


COMPUTER SCIENCES CORPORATION 

* 

cs 

68 


c* 



* 

cs 

69 

70 

c* 

LANGUAGE 

-FORTRAN 

* 

cs 

70 


c* 



* 

cs 

71 


c* 

DATE RELEASED SEPTEMBER 5, 1973 

* 

cs 

72 


c* . 



* 

cs 

73 


c* 

LATEST REVISION MARCH 1975 

* 

cs 

74 

75 

C* ********************* ************************************************* 

cs 

75 


c 




cs 

76 


c 




cs 

77 



DIMENSION 

X<1>, F(1 ) , DF ( 1 ), COEF (MAX, 4) , WKtl) 


cs 

78 


c 




cs 

79 

so! 

c 


SET UP WORKING AREAS 


cs 

80 
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LISTING OP DECK* CSDS 


PAGE 


3 




CARO NO* 


81 

85 


90 


95 


100 , 


105 , 


110 


115 


120 


C 

I ERR "0 

IP (IPT*NE*-1) GO TO 4 
IPT-0 

IF (IX.LT.4) GO TO 1 
GO TO 2 

1 IERR—1 
RETURN 

2 IX1«IX-1 

DO 3 I-1#IX1 

IF (X<I+l)-X(T>*GT.O) GO TO 3 

IERR-I+1 

RETURN 

3 CONTINUE 
NP1-IX+1 
IB1-NP1 

IB 2*1 Bl+NPl 
IB 3*1 B2 + NP1 + 1 
IB 4*1 B3+NP1 
I B5* I B4+N P 1 
IB 6* I B5+N Pl+1 
WK (1 ) *0* 

WK ( 2 ) *0 • 

UK( I B2) *0* 

WK (IB 3 ) *0 • 

IJK2-I82+NP1 

WK(IJK2)-0. 

I JK5-IB5+1 
WK(IJK5)*0. 

IJK5-IB5+2 

WK(IJK5>-0. 

WK(IB6)*0. 

IJK5-IB5+NP1 
WK (I JK5) *0* 

4 CONTINUE 
P-0* 

H-X<2)-X(ll 

F2--S 

FF-(F (2>-F (1))/H 
IF UX.LT.3) GO TO 10 


cs 

31 

cs 

32 

cs 

83 

cs 

84 

cs 

85 

cs 

86 

cs 

87 

cs 

88 

cs 

89 

cs 

90 

cs 

91 

cs 

92 

cs 

93 

cs 

94 

cs 

95 

cs 

96 

cs 

97 

cs 

98 

cs 

99 

cs 

100 

cs 

101 

cs 

102 

cs 

103 

cs 

104 

cs 

105 

cs 

106 

cs 

107 

cs 

108 

cs 

109 

cs 

110 

cs 

111 

cs 

112 

cs 

113 

cs 

114 

cs 

115 

cs 

116 

cs 

117 

cs 

118 

cs 

119 

cs 

120 
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CARD NO. 


121 

125. 

13o! 

135 i 

140' 

145 

150 

1 55* 

16<V 


DO 5 1-3# IX 
G-H 

H-X( I )-X ( I-l ) 

E-FF 

FF-(F ( I ) -F { I -1 ) ) / H 
COFF (I-1#1)-FF-E 
IJK3-IB3+I 

WK (I JK3) * ( G+H)* .66666666666667 

IJK4-IB4+I 

WK(I J K4) -H/3. 

IJK2-IB2+I 
WK(IJK2)-DF(I-2)/G 
WK (I)-DFCI)/H 
IJK1-IB1+I 

WK(IJKl)— DF(I-1)/G-DF(I-1)/H 

5 CONTINUE 

DO 6 1-3# IX 

IJK1-IB1+I 

IJK2-IB2+I 

C0EF(I-1#2)-WK(I)*WK{I)+WK(IJK1)*WK(IJK1)+WKCIJK2)*WK(IJK2) 
COEF ( 1-1# 3 ) -WK (T )*WK (IJK1+1)+WK(IJK1)*WK (IJK2 + 1) 

C OEF ( I—l# 4) - WK ( I ) *WK {I JK2+2 ) 

6 CONTINUE !' 

C ’ 

C NEXT ITERATION 

C 

7 IF CIX.LT.3) GO TO 10 

DO 8 1-3# IX 
IJK1-IB1+I-1 
IJKO-I-1 

WK(IJK1)-FF*WK(IJK0) 

IJK2-IB2+I-2 

IJKO-I-2 

WK(IJK2)-G*WK(IJK0) 

IJKO-I 

IJK3-IB3-H 

WK (I JK0)-1./{P*C0FP (I-l#2 )+WK(IJK3)-FF*WK( IJK1 )-G*WK(IJK2) ) 
IJK5-IB5+I 
IJKN-IJK5-1 
I JKO- I4KN-1 


CS 121 
CS 122 
CS 123 
CS 124 
CS 125 
CS 126 
CS 127 
CS 128 
CS 129 
CS 130 
CS 131 
CS 132 
CS 133 
CS 134 
CS 135 
CS 136 
CS 137 
CS 138 
CS 139 
CS 140 
CS 141 
CS 142 
CS 143 
CS 144 
CS 145 
CS 146 
CS 147 
CS 148 
CS 149 
CS 150 
CS 151 
CS 152 
CS 153 
CS 154 
CS 155 
CS 156 
CS 157 
CS 158 
CS 159 
CS 160 
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PAGE 


5 


CARD NO. 


161 

165 

170 

175 

180. 

185 

190 

195 

200 


8 


9 

10 

C 

C 

C 


11 


WK(IJK5)-C0EF (1-1,1) -WK( I JK1 ) *WK (I JKN)-WK( I JK 2 ) * WK ( IJ KO ) 
I J K4 * IB 4+1 

FF-P+COEF (I-1,3)+WK( I J K4 )-H*WK { I JKl ) 

G» H 

H-CQE F ( I— 1» 4 ) *P 
CONTINUE 
DO 9 1-3, IX 
J-IX-I+3 


IJK5-IB5+J 

IJK6-IJK5+1 

IJK7-IJK6U 

IJK1-IB1+J 

IJK2-IB2+J 


WK (I J K5 ) ' 

CONTINUE 

E-0 


WK ( J } *WK ( I J K5 )— WK (IJK1)*WK(IJK6)— 


WK(IJK2)+WK(IJK7) 


H-0 


COMPUTF U AND ACCUMULATE E 


DO 11 I- 2, 1 X 
G-H 

IJK5-IB5+I 

H-(WK(IJK5+1)-WK(IJK5) ) / ( X ( I )-X ( 1-1 ) ) 
IJK6-IB6+I 

WK(IJK6)«(H-G)*DF(I-1)*DF(I-1> 
E-E+WK(I JK6)*(H-G) 

CONTINUE 

G--H*DF(IX)*DF(TX) 

I J K6- IB6+NP1 
WK (I JK6) -G 
E-E-G4H 
G-F2 

F2-E + P*P 

IF (F2.GE.S.0R.F?.LF.G) GO TO 14 
FF-O. 

EJK6-IB6+2 

H- (WK ( I JK6+1 )— WK (IJK6) )/(X(2)— X(l)) 

IF (IX.LT.3) GO TO 13 
DO 12 1-3, IX 


CS 161 
CS 162 
CS 163 
CS 164 
CS 165 
CS 166 
CS 167 
CS 168 
CS 169 
CS 170 
CS 171 
CS 172 
CS 173 
CS 174 
CS 175 
CS 176 
CS 177 
CS 17B 
CS 179 
CS 180 
CS 181 
CS 182 
CS 183 
CS 184 
CS 185 
CS 186 
CS 187 
CS 188 
CS 189 
CS 190 
CS 191 
CS 192 
CS 193 
CS 194 
CS 195 
CS 196 
CS 197 
CS 198 
CS 199 
CS 200 
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LISTING OF DECK* CSDS 
CARO NO, 


PAGE 


6 


201 

205' 

210 

215 

220 

225’ 

230 

235 


12 

13 ; 

C 

c 

c 

c 


c 

c 

c 

c 

14 


15 


16 


G-H 

IJK6- IB6+-I 

H-(WK(IJK6+1) — WK { I JK 6 ) ) / ( X ( I )— X ( T— 1 ) > 

IJK1-IB1+I-1 

IJK2-IB2-H-2 

G«H-G-VK(IJK1)*WK (1-1 )-WK ( I JK2 ) *WK (1-2) 

FF«FF*G*WK(I)*G 

WK(I)-G 

CONTINUE 

H-E-P+FF 

IF (H.LE.O) GO TO 14 

UPDATE THE LAGRANGE MULTIPLIER P 
FOR THE NEXT ITERATION 

P-P + < S-F2)/( (SQRT(S/E)*P)+H) 

GO TO 7 


IF E LESS THAN OR EQUAL TO S, 

COMPUTE THE COEFFICIENTS AND RETURN. 


DO 15 1*2, NP1 
IJK6-IB6+I 

COEF ( I-1,1)-F(I-1)-P*WK(IJK6) 

IJK5-IB54I 

C0EF(I-1,3)-WK(TJK5) 

CONTINUE 
DO 16 1-2, IX 
H-Xd)-X(I-l) 

COEF (1-1,4)- (COEF (I* 3) -COEF ( 1-1, 3 ) ) fi 3.*H) 

COEF ( 1-1,2)- (COEF (1, 1J-C0FF ( 1-1,1) ) /H-(H*CQEF (1-1, 4)+C0EF (1-1,3) )* 
1H 

CONTINUE 

RETURN 

END 


CS 201 
CS 202 
CS 203 
CS 204 
CS 205 
CS 206 
CS 207 
CS 208 
CS 209 
CS 210 
CS 211 
CS 212 
CS 213 
CS 214 
CS 215 
CS 216 
CS 217 
CS 218 
CS 219 
CS 220 
CS 221 
CS 222 
CS 223 
CS 224 
CS 225 
CS 226 
CS 227 
CS 228 
CS 229 
CS 230 
CS 231 
CS 232 
CS 233 
CS 234 
CS 235 
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STING OF 

DECK t 

AVI SLR 

PAGE 

1 

,RD NO. 





1 


SUBROUTINE AXISLB ( I AX I S > XOR G* YQRG , PPG* SC ALE/ D 1ST » HT, Nn IG ) 

XL 

1 


C 


XL 

2 


C 

ROUTINE TO DRAW SCALE VALUES ON AXIS 

XL 

3 


C 


XL 

4 

5 

c 

CODFD BY — HARRY L. MORGAN NASA/LARC/TAO/AAB 

1983 XL 

5 


c 


XL 

6 


c 

PARAMETER DEFINITION 

XL 

7 


c 

IAXIS ■ 1 - DRAW Y-AXIS SCALES 

XL 

8 


c 

IAXIS ■ 2 - DRAW Y-AXIS SCALES 

XL 

9 

10 

c 

XORG* YQRG - X,Y COORDINATES Op BEGINNING OF AXIS 

XL 

10 


c 

ORG - SCALE VALUE AT BEGINNING OF AXIS 

XL 

11 


c 

SCALE - SCALE VALUE °ER INCH 

XL 

12 


c 

DIST - LENGTH OF AXIS 

XL 

13 


c 

HT - HEIGHT OF VALUES DRAWN ON AXIS 

XL 

14 

15 

c 

ND IG - NUMPER OF SIGNIFICANT FIGURES TO RIGHT OP DECIMAL 

POINT XL 

15 


c 

NDIG — 1 WILL DROP DECIMAL POINT 

XL 

16 


c 


XL 

17 


c 

COMPUTE ROUTINE CONSTANTS 

XL 

18 



WADJ-0.0 

XL 

19 

20 


IF (HT.LE.0.1) WA0J-0.01 

XL 

20 



C0NS»6.*HT/7.*WADJ 

XL 

21 



C 0NS2 *2 • +HT/7* 

XL 

22 



AD J»0 .4*HT 

XL 

23 



IF (ADJ.LT.0.1) AD J * 0* 1 

XL 

24 

25 


NL-IFIX<DISH-0.0001)+1 

XL 

25 



XO-XORG 

XL 

26 



YO-YORG 

XL 

27 



VALUE-ORG 

XL 

28 


c 

POSITION PEN AT BEGINNING OF AXIS 

XL 

29 

30 


CALL CALPLT (XOPG.YORG, 3) 

XL 

30 


c 

DRAW AXIS VALUE AT FACH INCH MARK ALONG AXIS 

XL 

31 



DO P I* I# NL 

XL 

32 



VALUE -SCALE* FLOAT (I-l)+nRG 

XL 

33 


c 

COMPUTE NUMBER OF SIGNIFICANT FIGURES IN VALUE 

XL 

34 

35 


NDG-NDIG 

XL 

35 



V-ABS (VALUE) 

XL 

36 



IF (V-l.OE-5) It 1> 2 

XL 

37 


1 

NDG— 1 

XL 

38 



NSIG-1 

XL 

39 

40 


V-0.0 

XL 

40 
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CARD NO. 


PAGE 2 


41 


45 


50 


55 


60 


65 


VALUE-0.0 
GO TO 5 

2 NSIG-NPG+1 

3 S-V+S IGN( 1, E-7*V) 

N-IFIX(S) 

IF (N) 5*5,4 

4 NSIG-NSIG+1 
V-0. 1C*V 

GO TO 3 

5 IF ( VALUE .LT. 0.0 ) NSIG-NSIG+1 

S-C0NS*FL0AT(N$IG)-C0NS2 

C POSITION AND DRAW AXIS VALUE 

IF (IAXIS.E0.2) GO TO 6 
XT -X0-0 .5*S 
YT-YQ-HT-ADJ 
X0«X0*1. 

GO TO 7 

6 XT-XO-S-ADJ 
YT-Y0-0.5+HT 
YO-YO+1. 

7 CALL NUMBER ( XT* YT* HT, VALUE, 0.* NDG ) 

C ADVANCE TO NEXT INCH MARK 

8 . CONTINUE 

C . RETURN TO BEGINNING OF AXIS 

CALL CALPLT < XOR G* YORG# 3 ) 

C RETURN TO CALLING PROGRAM 

RETURN 
END 


XL 

41 

XL 

42 

XL 

43 

XL 

44 

XL 

45 

XL 

46 

XL 

47 

XL 

48 

XL 

49 

XL 

50 

XL 

51 

XL 

52 

XL 

53 

XL 

54 

XL 

55 

XL 

56 

XL 

57 

XL 

58 

XL 

59 

XL 

60 

XL 

61 

XL 

62 

XL 

63 

XL 

64 

XL 

65 

XL 

66 

XL 

67 

XL 

68- 
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PAGE 


CARO NO, 
1 


35 


40 


SUBROUTINE COEFSY (XO# YO*HT#ISYM) 

ROUTINE TO DRAW A c RODYNAMIC COEFFICIENTS 


CODED BY — HARRY L. MORGAN 
PARAMETER DEFINITION 


NASA/LARC/TAD/AAB 


1983 


DRAW LETTER C 

CALL CHARACT (X0#Y1* HT#2H)C*0. >2*0.0) 
CALL CHARWH <WD# XI# YU HT> 2H)C# 2*0.0) 
X1»XQ+1. 15+WD 
Yl* YQ-HT 

DRAW SUBSCRIPT OF C 
IF (ISYM.EQ.l) CALL 
IF (ISYM.EQ.2) CALL 
IF t ISYM.E0.3) CALL 
IF (ISYM.EQ.4) CALL 


CHARACT (XI, Y1*HT#2H>L#0. #2*0.0) 
CHARACT <X1#Y1#HT*2H)D#0. #2*0.0) 
CHARACT (X1*Y1#HT#2H(M#0.#2#0.0) 
LAM ( XI# Y1#HT) 


CY 

CY 

CY 

CY 

CY 

CY 

CY 


CY 

CY 



C 

XO - X COORDINATE OF LEFT-HAND 

EDGE OF SYMBOL 

CY 

8 


C 

YO - Y COORDINATE OF CENTERLINE 

OF MAIN LETTER OF SYMBOL 

CY 

9 

10 

c. 

HT - SYMBOL HEIGHT 



CY 

10 


C 

ISYM • 1 - LIFT 

WIDTH ■ 

2.0*HT 

CY 

11 


C 

ISYM • 2 - DRAG 

WIDTH - 

2 . 1*HT 

CY 

12 


C 

ISYM ■ 3 - PITCHING MOMENT 

WIDTH « 

2.2*HT 

CY 

13 


C . 

ISYM - 4 - ROLLING MOMENT 

WIDTH » 

1.7*HT 

CY 

14 

15. 

c ; 

ISYM « 5 - YAWING MOMENT 

WIDTH - 

1.9+HT 

CY 

15 


c 

ISYM ■ 6 - SIDE FORCE 

WIDTH - 

2.1+HT 

CY 

16 


c 

ISYM - 7 - LIFT/DRAG 

WIDTH - 

2.6+HT 

CY 

17 


c 

ISYM - 8 - ALPHA# DEG 

WIDTH - 

3.5*HT 

CY 

18 


c 

ISYM - 9 - BETA*DEG 

WIDTH - 

3.4*HT 

CY 

19 

20 

c 




CY 

20 


c 

INITIALIZE CHARACTER SETS 



CY 

21 



CALL CHARST1 



CY 

22 



CALL CHARST7 



CY 

23 



CALL CHARST2 



CY 

24 

25 

c 

POSITION PEN AT START OF SYMBOL 



CY 

25 



CALL CALPLT (X0*Y0*3) 



CY 

26 



Y1-Y0-0.5*HT 



CY 

27 



IF (ISYM.EQ.7) GO TO 1 



CY 

28 



IF (ISYM. EQ.8.0R. ISYM. EQ.9) GO TO 

2 


CY 

29 

30 


IF ( ISYM. GT .9) GO TO 3 



CY 

30 


CY 31 
CY 32 
CY 33 
CY 34 
CY 35 


36 

37 


CY 38 
CY 39 
CY 40 
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41 


IF (ISYH.EQ.5) CALL CHARACT CX1,Y1#HT#2H(N#0.# 2*0.0) 

CY 

41 



IF (ISYM.EQ.6) CALL CHARACT {X1#Y1#HT#2H)Y#0.#2#0.0) 

CY 

42 



GO TO 4 

CY 

43 


C 

DRAW L/D 

CY 

44 

45 

1 

CALL CHARACT (Xrj# Yl, HT,4H) L/D#0.#4,0.25) 

CY 

45 



GO TO 4 

CY 

46 


C 

DRAW ALPHA# DEG OR BETA #DFG 

CY 

47 


2 

IF (ISYM.EG.8) CALL CH AR ACT ( XO# Yl# HT# 11H*7( A*2 < • DEG# 0.# 1 1# 0. 20) 

CY 

48 



IF ( I SYM.EQ.9) CALL CHARACT ( XO# Yl# HT# I1H$7 ( B$2 ( #D6G# 0. # 11# 0. 20) 

CY 

49 

50. 


GO TO 4 

CY 

50 


C 


CY 

51 


C 

ADD ADDITIONAL SYMBOLS HERE AS DESIRED 

CY 

52 


3 

CONTINUE 

CY 

53 


C 


CY 

54 

55 

c 


CY 

55 


c 

RETURN PEN TO START OF SYMBOL 

CY 

56 


4 

CALL CALPLT (XO#YO,3) 

CY 

57 



CALL CHNGSET (1) 

CY 

58 


c 

RETURN TO CALLING PROGRAM 

CY 

59 

60* 


RETURN 

CY 

60 



END 

CY 

61- 
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LISTING OF DECK: LAM 


PAGE 


1 

CARD NO. 


1 


SUBROUTINE LAM (XOfYOfHT) 

LM 

1 


C 


LM 

2 


c 

ROUTINE TO DRAW ROLLING MOMENT SUBSCRIPT 

LM 

3 


c 


LM 

4 

5 

c 

CODED BY — HARRY L. MORGAN NASA/L ARC/TAO / AAB 1983 

LM 

5 


c 


LM 

6 


c 

PARAMETER DEFINITION 

LM 

7 


c 

XOfYO - X# Y COORDINATES OF LOWER LEFT-HAND CORNER OF SYMBOL 

LM 

8 


c 

HT - HEIGHT OF SYMBOL 

LM 

9 

10 

c 


LM 

10 



DIMENSION XM ( 17 ) » YM(17> 

LM 

11 



DATA XM/3.#13.*0.f0.f4.fl0.fl3«jl4.fl3.fl2.2fl0.f4.#l.#l.fl4.4*3.6 

LM 

12 



If 3. / 

LM 

13 



DATA YM/43.#43.>7»#2*j0»f0.f2*f 5*f5*f2.6* 1>> 1 • > 2. 6* 7 *f44.*44.> 43*/ 

LM 

14 

15 


HTH-HT/44. 

LM 

15 



CALL CAL PL T (X0fYO,3) 

LM 

16 



DO 1 I » 1 f 1 7 

LM 

17 



K-2 

LM 

18 



IF (I.EQ.l) K-3 

LM 

19 

20 


X*X0+(XM(I)+2.)*HTH 

LM 

20 



Y«YO+YM(I)*HTH 

LM 

21 


1 . 

CALL CALPLT <X,Y,K) 

LM 

22 



CALL CALPLT (X0fY0f3) 

LM 

23 



RETURN 

LM 

24 

25 


END 

LM 

25- 
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1 


SUBROUTINE GRIOLN ( XO, YO, XL, YH, NDI V) 

GL 

1 


C 


GL 

2 


C 

ROUTINE TO DRAW AREA OF GRID LINFS 

GL 

3 


C 


GL 

4 

5 

C 

CODED BY — MARRY L. MORGAN NA 5A/LARC/TAD/ AA0 

1983 GL 

5 


C 


GL 

6 


C 

PARAMETER DEFINITION 

GL 

7 


C 

XO,YO - X, Y coordinates of lower left-hand corner of 

GRID AREA GL 

8 


C 

XL - LENGTH OF GPID AREA 

GL 

9 

10 

c 

YH - HEIGHT OF GRID APFA 

GL 

10 


c 

ND IV - NUMBER OF GRID LINES PER INCH 

GL 

11 


c 


GL 

1? 


c 

COMPUTE NUMBER OF GRID LINES 

GL 

13 



XD-FLOAT(NDIV) 

GL 

14 

15 


IF (XD.LE.O.O) XD«1. 

GL 

15 



NGX«IFIX(XD*XL*0.0001)+1 

GL 

16 



NGY»IFIX(XD*YH+0.0001)+1 

GL 

17 


c 

DRAW GRID LINES PARALLEL TO Y-AXIS 

GL 

18 



CALL CALPLT (XO,Y0,3) 

GL 

19 

20 


D-0. 

GL 

20 



IF (NGX.GT.l) D-XL/FL0AT(NGX-1) 

GL 

21 



Xl-XO 

GL 

22 



Yl-YO 

GL 

23 



Y2-Y0+YH 

GL 

24 

25 


L *0 

GL 

25 



DO 3 I-1,NCX 

GL 

26 



L-L + l 

GL 

27 



X2-X14D 

GL 

28 



IF (L-l) 1,1,? 

GL 

29 

30 

1 

CALL CALPLT (X1,Y2,2) 

GL 

30 



IF (I.EQ.NGX) GO TO 3 

GL 

31 



CALL CALPLT (X2,Y2,3) 

GL 

32 



GO TO 3 

GL 

33 


2 

CALL CALPLT (X1,Y1,2) 

GL 

34 

35 


IF ( I .EO.NGX) GO TO 3 

GL 

35 



CALL CALPLT (X2,Y1,3) 

GL 

36 



L-0 

GL 

37 


3 

X 1*X 2 

GL 

38 


c 

DRAW GRID LINES PARALLEL TO X-AXIS 

GL 

39 

40 


CALL CALPLT (X0,Y0,3) 

GL 

40 
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2 


CARD NO. 
41 

45^ 


D*0 • 

IF (NGY.GT.l) D*YN/FL OAT( NGY- 

Yl-YO 

Xl-XO 

X2-X0+XL 



L-0 

DO 6 I-1>NGY 

50 

4 

L-L+l 
Y2 »Y1 *D 

IF (L-l) 4,4,5 
CALL CALPLT <X2,Y1,2> 
IF (I.EO.NGY) GO TO 6 
CALL CALPLT <X2,Y2,3) 
GO TO 6 

55 

5 

CALL CALPLT (Xl,Yl,2) 
IF (I.EO.NGY) GO TO 6 
CALL CALPLT (X1,Y2,3) 
L»0 


6 

Y1*Y2 

60 

C 

FRAME GRID AR FA 

Xl-XO+XL 

Yl-YO+YH 

65 


CALL CALPLT (X0,Y0,3) 
CALL CALPLT (X1,Y0,2) 
CALL CALPLT (X1,Y1,2) 
CALL CALPLT (X0,Y1,2) 



CALL CALPLT (X0,Y0,2) 
CALL CALPLT <XO,YO,3) 

70 

c ! 

RETURN TO CALLING PROGRAM 

RETURN 

END 


GL 41 
GL 42 
GL 43 
GL 44 
GL 45 
GL 46 
GL 47 
GL 48 
GL 49 
GL 50 
GL 51 
GL 52 
GL 53 
GL 54 
GL 55 
GL 56 
GL 57 
GL 58 
GL 59 
GL 60 
GL 61 
GL 62 
GL 63 
GL 64 
GL 65 
GL 66 
GL 67 
GL 68 
GL 69 
GL 70 
GL 71 
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LISTING OF DECK* RUNKEY 


PAGE 


1 


CARD NO. 
1 


5 


10 


15 


20 


25' 


30 


35 ' 


40 


SUBROUTINE RUNKEY (XO# Y0# RUN, LSYM, ISI ZE> HT#NRNMAX> 

C 

C ROUTINE TO DRAW. TITLE BLOCK CONTAINING SYMBOLS AND RUN NUMBERS 

C CODED BY — HARRY L. MORGAN NASA/LARC /TAD/AAB 1983 

C 

C PARAMETER DEFINITION 

C XO# YO - X#Y COORDINATES OF LOWER LEFT-HAND CORNER OF KEY 

C RUN - RUN NUMBERS 

C LSYM - SYMBOL ORDER 

C ISIZE - SYMBOL SIZE 

C HT - SYMBOL HEIGHT 

C NRNMAX - MAXIMUM NUMBER OF ALLOWABLE RUNS PER SHEET 

C 

DIMENSION RUN f 1 ) # LSYMC1) 

C WORK ARRAYS 

COMMON / PLT/ KRUN(10)#NRUN(10)#LRUN(22)#JSYM(22) 

C 

INTEGER RUN 

C INITIALIZE ROUTINE CONSTANTS 

DATA XM/1.0/#XMARG/0.1/#XNSP/2.5/#SPC/0.2/ 

C DEFINE LETTERING SIZE 

SZ *HT +XM 

C INITIALIZE CHARACTER SET 

CALL CHARST2 
C 

C DETERMINE THE NUMBER OF NON-ZERO RUNS AND CORRESPONDING SYMBOLS 

C 

NR-0 

DO 2 1*1# NRNMAX 
IF (RUN(I)) 1 # Zf 1 

1 NR*NR4l 

LRUNt NR)-RUN(I) 

JSYM ( NR )*LSYM (I ) 

2 CONTINUE 

IF (NR.EQ.O) RETURN 

C . 

C COMPUTE WIDTH OF WORDS SYMBOL AND 99999 

C 

CALL CHARWH (SYMV#HA#HBS# SZ# BH) S(YMBOL# 8#$PC) 


KY 1 
KY 2 
KY 3 
KY 4 
KY 5 
KY 6 
KY 7 
KY 8 
KY 9 
KY 10 
KY 11 
KY 12 
KY 13 
KY 14 
KY 15 
KY 16 
KY 17 
KY 18 
KY 19 
KY 20 
KY 21 
KY 22 
KY 23 
KY 24 
KY 25 
KY 26 
KY 27 
KY 28 
KY 29 
KY 30 
KY 31 
KY 32 
KY 33 
KY 34 
KY 35 
KY 36 
KY 37 
KY 38 
KY 39 
KY 40 
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41 


CALL CHARWH (WD,HA, MB, SZ, 5H99999, 5 $ SPC ) 

KY 

41 



XS-X0+0.5*SYMW+XMARG 

KY 

42 



XR-X0*2.*SYMW+0. 5*WD+XMARG 

KY 

43 



YSR*YO+XMARG 

KY 

44 

45 

C 


KY 

45 


C 

ORAW SYMBOL AND RUN NUMBER 

KY 

46 


C 


KY 

47 



DO 6 L-1,NR 

KY 

48 



I-NR+l-L 

KY 

49 

50 ' 


ENCODE ( 10# 7# MRUN) LRUN(I) 

KY 

50 



DECODE (10,8,MRUN) KRUN 

KY 

51 


C 

FIND NUMBER OF DIGITS IN RUN NUMBER 

KY 

52 



DO 3 J-l>10 

KY 

53 


3 

NR UN ( J ) «1H 

KY 

54 

55 


KDIG-0 

KY 

55 



DO 4 J*l>10 

KY 

56 



IF (KRUNtJ ).E0.1H ) GO TO 4 

KY 

57 



KDIG-KDIG+1 

KY 

58 


C 

SHIFT RUN NUMBER FROM RIGHT TO LEFT JUSTIFIED 

KY 

59 

60 


NRUNt KDI G )«KRUN (J ) 

KY 

60 


4 

CONTINUE 

KY 

61 


C 

REFORMAT RUN NUMBER 

KY 

62 



ENCODE <10,8,MRUN> NRUN 

KY 

63 



DECODE { 10# 9* MRUN ) IRUN 

KY 

64 

65 

c 

DRAW SYMBOL 

KY 

65 



Y2*YSR+0. 5*SZ 

KY 

66 



CALL PNTPLT ( XS, Y2# J SYM (T ) , ISIZ E ) 

KY 

67 


c 

DRAW RUN NUMBER 

KY 

68 



IF (KDIG.EQ.O) GO TO 5 

KY 

69 

70, 


CALL CHARWH t WD,HA#HB> SZ, IRUN, KOI G, SPC ) 

KY 

70 



XX *XR-0. 5* WD 

KY 

71 



CALL CHARACT ( XX, YSR , SZ,I RUN,0., KDIG, SPC > 

KY 

72 


5 

YSR-YSR+XNSP*SZ 

KY 

73 


6 

CONTINUE 

KY 

74 

75 

C 


KY 

75 


c 

DRAW AND UNDERLINE WORD SYMBOL AND RUN 

KY 

76 


c 


KY 

77 



CALL CHARWH <WD,HA,H8, SZ, 5H)R (UN,5»SPC ) 

KY 

78 



IF (HB.LT.HBS) HBS-HB 

KY 

79 

80 


XX-XS-0.5* (SYMW4-XMARG) 

KY 

80 
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81 


CALL CALPLT (XX, Y$R,3) 

KY 

81 



XX-XS40.5*<SYMW+XHARG) 

KY 

82 



CALL CALPLT (XX,YSR,2) 

KY 

83 



XX-XR-0.5MV0+XHARG) 

KY 

84 

85 


CALL CALPLT (XX,YSR,3) 

KY 

85 



XX-XR+ 0.5MWD+XMARG) 

KY 

86 



CALL CALPLT (XX,YSR,2) 

KY 

87 



YSR«YSR-HBS+0.5*XMARG 

KY 

88 



XX*XS-0.5*SYMW 

KY 

89 

90 


CALL CALPLT <XX,YSR,3) 

KY 

90 



CALL CHARACT ( XX, YSR,SZ,BH)S CYMBOL, 0.,8>SPC ) 

KY 

91 



XX«XR-0.5*WD 

KY 

92 



CALL CHARACT (XX, YSR,SZ,5H)R (UN,0.,5,SPC) 

KY 

93 


C 


KY 

94 

95 

C 

RETURN TO TITLE BLOCK ORIGIN 

KY 

95 


C 


KY 

96 



CALL CALPLT (XO,Y0,3) 

KY 

97 



CALL CHARST1 

KY 

98 


C 


KY 

99 

100 


RETURN 

KY 

100 


C 


KY 

101 


7 

FORMAT (110) 

KY 

102 


8 

FORMAT (10A1) 

KY 

103 


9 

FORMAT (A10) 

KY 

104 

105 


END 

KY 

105 
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APPENDIX B 


DESCRIPTION OF INPUT DATA FOR WIND-TUNNEL DATA 
PLOTTING PROGRAM PLOTWD 

This appendix contains a description of the input requirements for the 
wind-tunnel data plotting program PLOTWD. The input data is divided into two 
basic parts: (1) the plotting setup information and (2) the plotting sheet 

information. The setup information defines the data variables, axis scales 
and positions, and the type of data fairing desired. All setup variables are 
input as floating point quantities with a format of F10.0, except the YLABEL 
and XLABEL variables which are axis labels and are input with a format of A10. 


Plotting Setup 


CARD 

VARIABLE VALUE 

DESCRIPTION 

1 

TEST 

Test number 


2 

SHEETW 

- 

Total width of plotting sheet, in. 


SHEETH 

- 

Total height of plotting sheet, in. 


SPACE 

— 

Space between plotting sheets, in. 

3 

ISYM 

- 

Starting symbol number (See table I.) 


ISIZE 

1 

Small symbol size 



2 

Medium symbol size 



3 

Large symbol size 

4 

IOP 

0 

Plot data points only 



1 

Plot data points and fair with tension 
spline 


TENSION 

- 

Tension factor, a 

5 

IEDIT 

0 

Do not sort or edit data 



1 

Sort and edit data 


TOLR 


Edit tolerance of independent 
test variable (T0LR=0 for sort only) 

6 

IW 

- 

Data array location of independent 
test variable 
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CARD VARIABLE VALUE 


DESCRIPTION 


7 NPLOT 


Number of plots per sheet 


8 YTAPE 

YOFFSET 

YSCALE 

YLABEL 

XTAPE 

XOFFSET 

XSCALE 

XLABEL 


Data array location of y-axis variable 
Sheet height location of y-axis 
origin, in. 

Scale value per inch for y-axis 
Label for y-axis 

Data array location of x-axis variable 
Sheet width location of x-axis 
origin, in. 

Scale value per inch for x-axis 
Label for x-axis 


Card 8 is input NPLOT times and NPLOT is limited to a value of 10 for each 


setup • 


NAMELIST VARIABLE 
$ SHEET 
NO 
RUN 

NEWCASE 

$END 


Plotting Sheet Namelist 

VALUE DESCRIPTION 

- Namelist label 

Sheet number 


0 or default 

1 


Run numbers to be plotted on 
sheet 

A new SHEET namelist follows 
A new setup deck follows 

End of namelist 


Note that no part of the namelist may be in column 1 and that a maximum of 10 
run numbers may be input per sheet. 
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DESCRIPTION OF OUTPUT FOR WIND-TUNNEL DATA PLOTTING PROGRAM PLOTWD 
This appendix contains a description of the typical output for the wind- 
tunnel data plotting program PLOTWD. A sample four-page output is presented 
in table IV and was generated from the upper set of sample input data pre- 
sented in table III. Page 1 of the output is a printout of the input setup 
data as described in appendix B. This page should be carefully checked after 
each program execution to insure that the desired plotting variables and 
scales have been properly input. Page 2 of the output is a list of the run 
numbers contained in the data copied from the input file TAPE1 to the random- 
access file TAPE2. The total number of data points copied to the random- 
access file is printed following the list of run numbers. The quantity NPMAX 
and the dimensioned size of the array NPT in the main program PLOT should be 
equal to or greater than the total number of data points copied to file TAPE2 • 
.Page 3 of the output is a list of the sheet and corresponding run numbers 
plotted. A run number with the value zero indicates that no data were plotted 
with the symbol corresponding to the order of the run number in the list. If 
the user— specif ied run number is not available on the random— access file* a 
message will be printed stating that the specified run could not be found. In 
addition, if during the plotting of the scaled data any data points fall out- 
side the sheet boundaries, a message will be printed stating the number of 
data points outside the boundaries. These outside data points will not be 
plotted and, therefore, may require that the user redefine either the axis 
scale factors or the sheet height and width. 

If the user selects the sort and edit option, page 4 will be output which 
contains a summary of the sort and edit information for each run called during 
the program execution. The first set of values listed in this summary are the 
input values of the independent variable t for the particular run. The sorted 
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t values listed next represent the results of the sorting procedure which 
simply reorders the input values in a monotonically increasing order. The 
sorted t values for a particular run will be the same regardless of the number 
of times the run is plotted. The edited and deleted t values listed last 
represent the results of the editing procedure which deletes all but one of 
the multiple data points within the specified tolerance TOLR. The data points 
deleted are selected based on the input y values for the particular part of 
the figure and, therefore, may be different for each part. The edited and 
deleted t values listed are only for the first part of the figure. The 
deleted points are not used during the spline curve fairing process, but 
nevertheless, they are plotted with the appropriate symbol. 

If the user replaces the call to subroutine CURPLT in the main program 
PLOT with a call to subroutine LSQPLT, the sample output summary information 
presented in table V will be listed for the least-squares polynomial curve fit 
option and in table VI for the least-squares cubic-spline curve fit option. 

The parameter IPRINT in subroutine LSQPLT must also be set equal to 1 for the 
summary information to be listed. This information will appear in the output 
immediately following the print of the particular sheet and corresponding run 
numbers on page 3 of the basic output. The summary information will be listed 
each time the subroutine is called; therefore, the user is cautioned that a 
large amount of printed output can easily be generated even if a moderate 
number of plots and runs are specified. 

The summary information presented in tables V and VI consists of a tabu- 
listing of the input t, x, and y values, the new x and y values 
generated by the least-squares curve fit, and the differences between input 
and new x and y values. Immediately following the tabulated listing, the 
sum-of-the- squares of the differences between the input and new x and y values 
are listed as an error-in-x and -y, respectively. If the independent 
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variable t and the x-axis variable are the same as it is for the sample output 
case, the error-in-x will be very small. If the least-square polynomial curve 
fit option is chosen, the error-in-x and -y will be a function of the degree 
of the polynomial. If the least-square cubic-spline option is chosen, the 
error-in-x and -y will always be equal to the product of the number of input 
points times the square of the standard deviation DF. 
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TABLE 1. - LIST OF STANDARD SYMBOL AND DASHED LINE COMBINATIONS 
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TABLE II.- LIST OF STANDARD AERODYNAMIC SYMBOLS DRAWN BY SUBROUTINE COEFSY 


ISYM I COEFFICIENT SYMBOL I WIDTH/HEIGHT 


Lift 

Drag 

Pitching Moment 
Rolling Moment 
Yawing Moment 
Side Force 
Lift- Drag Ratio 
Angle of Attack 
Angle of Sideslip 


c L 

C ra 

Cz 

c„ 


2.0 

2.1 

2.2 


1.9 

2.1 

2.6 

3.5 

3.4 









»] 


TABLE III.- LISTING OF SAMPLE INPUT CASES 




CASE 1 - Single Plot with Three Variables 

496. 

20. 23.. 10. 

1. 3. 

1. 0. 

1. .5 

1. 

3. 

6 . 2 . .4 CL. 

7. 11. .2 CD 

8. 20. .4 CM 

SSHEET NO=206*PUN=39,59S 

1. 4. 2. ALPHA 

l. A. 2. ALPHA 

1. 4. 2. ALPHA 

CASE 2 - Four Plots with One Variable 

496. 

12. 12. 10. 

1. 2. 

1. 0. 

1. .5 
1. 

4 . 

8. 2. .4 CM 

9. 1. 4. L/D 

2. 8. .8 CL 

3. 7. .2 CD 

SSHEET N0 = 206 ♦ RUN = 3 9 * 59$ 

2 • 1 . • 8 CL 

U 8. 10. ALPHA 

1. 1. 10. ALPHA 

2 • 8 . .8 CL 




TABLE IV.- SAMPLE OUTPOT FOR CASE 1 INPUT 


PAGE 1 OUTPUT 

1 PLOTTING SETUP DECK FOR 

TEST - 

A96 







SHEET WIDTH ■ 

20.00 

SHEET 

HEIGHT ■ 

23.00 SPACE 

BETWEFN SHEETS 

■ 10.00 



STARTING SYMBOL 

NUMBER ■ 

1 

SYMBOL SIZE IS LARGE 






DATA WILL BE FAIRED AND 

PLOTTED 

USING SPLINE WITH TENSION 

« 0. 

00 




DATA WILL BE EDITEO WITH 

AN EDIT 

TOLERANCE 

• .500 






TAPE LOCATION OF 

INDEPENDENT VARIABLE • 

1 






NUMBER OF PLOTS 

PER SHEET « 3 








PLOT NO. 

YT APE 

YOFFSET YSCALE YLABEL 


XTAPE 

XOFFSET 

XSCALE 

XL ABEL 

1 

6 

2.00 

• 

AOOO CL 


1 

A. 00 

2.0000 

ALPHA 

2 

7 

11.00 

• 

2000 CD 


1 

A. 00 

2.0000 

ALPHA 

3 

8 

20.00 

« 

AOOO CM 


1 

A. 00 

2.0000 

ALPHA 

PAGE 2 OUTPUT 

THE FOLLOWING 

LIST OF 

RUNS ARE 

AVAILABLE 

ON THE DATA DISK FOR 

TEST NO. 

A96 



2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

1? 

13 

1A 

15 

16 

17 

18 

19 

20 

21 

22 

23 

2A 

25 

27 

28 

29 

30 

31 

32 

33 

3A 

35 

36 

37 

38 

39 

AO 

A1 

A2 

A3 

AA 

A5 

A6 

A7 

A8 

A9 

50 

51 

52 

53 

5A 

55 . 

56 

57 

58 

59 

60 

61 

62 

63 

6A 

65 

66 

67 

68 

69 

70 

71 

72 

73 

7A 

75 

76 

77 

78 

79 

80 

81 

8 2 

83 

BA 

85 

86 

87 

88 

89 

90 . 

91 

92 

93 

9A 

95 

96 

97 

98 

99 

100 

101 

102 

103 

10A 

105 

106 

107 

108 

109 

110 

111 

112 

113 

11A 

115 

116 

117 

118 

119 

120 

121 

122 

123 

12A 

125 

126 

127 

128 

129 

130 

131 

132 

133 

13 A 

135 

136 

137 

138 

139 

1 AO 

1A1 

1 A 2 

1 A 3 

1 AA 

1A5 

1A6 

1A7 

1A8 

1A9 

150 

151 

152 

153 

15 A 

155 

156 

157 

158 

159 

160 

161 

162 

163 

16 A 

165 

166 

167 

168 

169 

170 

171 

172 

173 

17A 

175 

176 

177 

178 

179 

180 

181 

182 

183 

18A 

185 

186 

187 

188 

189 

190 

191 

192 

193 

19A 

195 

196 

197 

198 

199 

200 

201 

202 

203 

20A 

205 

206 

207 

208 

209 

210 

211 

212 

213 

21A 

215 

216 

217 

218 

219 

220 

221 

222 

223 

22A 

225 

226 

227 

228 

229 

230 

231 

23 2 

233 

23A 

235 

236 

237 

238 

239 

2A0 

2A1 

2A2 

2A3 

2AA 

2A5 

2A6 

2A7 

2A8 

2A9 

250 

251 

252 

253 

25A 









TOTAL NUMBER OF 

DATA POINTS ■ 

A901 







PAGE 3 OUTPUT 

DATA PLOTTED 










SHEET NO = 206 RUNS = 

39. 59. 0. 

0* 

0. 

0. 0. 

0* 

0. 0 









TABLE IV.- CONCLUDED 


PAGE 4 OUTPUT 



THE FOLLOWING 

RUNS HAVE BEEN EDITED 




PUN NO. 

- 39 

TOLERANCE 

■ .500 





INPUT 

SORTED 

EDITED 

deleted 

POINT NO 

. T-VALUF 

POINT NO. 

T-VALUE 

POINT NO. 

T-VALUE 

POINT NO. T-VALUE 

1 

-6.1700 

1 

-6.1700 

1 

-6.1700 

4 .3500 

2 

-4.1800 

2 

-4.1800 

2 

-4.1800 


3 

-1.8800 

3 

-1.8800 

3 

-1.8800 


4 

.3500 

20 

.1700 

20 

.1700 


C 

2.0100 

4 

.3500 

5 

2.0100 


6 

4.2800 

5 

2.0100 

6 

4.2800 


7 

6.3700 

6 

4.2800 

7 

6.3700 


8 

8.4500 

7 

6.3700 

8 

8.4500 


9 

10.3200 

8 

8.4500 

9 

10.3200 


10 

12.4100 

9 

10.3200 

10 

12.4100 


11 

14.4000 

10 

12.4100 

11 

14.4000 


12 

16.8000 

11 

14.4000 

12 

16.8000 


13 

18.6700 

12 

16.8000 

13 

18.6700 


14 

20.8600 

13 

18.6700 

14 

20.8600 


15 

22.7500 

14 

20.8600 

15 

22.7500 


16 

24.7800 

15 

22.7500 

16 

24.7800 


17 

26.6700 

16 

24.7800 

17 

26.67C0 


18 

28.7500 

17 

26.6700 

18 

28.7500 


19 

30.7000 

18 

28.7500 

19 

30.7000 


20 

.1700 

19 

30.7000 




PUN NO. 

a 50 

TOLERANCE 

■ .500 





INPUT 

SORTED 

EDITED 

DELETED 

POINT NO 

l. T-VALUE 

POINT NO. 

T-VALUE 

POINT NO. 

T-VALUE 

POINT NO. T-VALUE 

1 

-6.0200 

1 

-6.0200 

1 

-6.0200 

4 .2800 

2 

-4.4000 

? 

-4.4000 

2 

-4.4000 


3 

-2.2200 

3 

-2.2200 

3 

-2.2200 


4 

.2800 

4 

.2800 

20 

.2900 


5 

2.0000 

20 

.2900 

5 

2.0000 


6 

4.8700 

5 

2.0000 

6 

4.8700 


7 

6.3000 

6 

4.8700 

7 

6.3000 


a 

8.2000 

7 

6.3000 

8 

8.2000 


9 

10.5000 

8 

8.2000 

9 

10.5800 


10 

12.6400 

9 

10.5800 

10 

12,6400 


11 

14.9700 

10 

12.6400 

11 

14.9700 


12 

16.6100 

11 

14.9700 

12 

16.6100 


13 

18.4600 

12 

16.6100 

13 

18.4600 


14 

20.8900 

13 

18.4600 

14 

20.8900 


15 

22.8300 

14 

20.8900 

15 

22.8300 


16 

25.3700 

15 

22.8300 

16 

25.3700 


17 

26.9800 

16 

25.3700 

17 

26.9800 


18 

28.9200 

17 

26.9800 

18 

28.9200 


19 

31.3400 

18 

28 .9200 

19 

31.3400 


20 

.2900 

19 

31.3400 





AIL PLOTTING COMPLETED 


120 




TABLE V.- SAMPLE OF LEAST -SQUARES POLYNOMIAL CURVE FIT 


DATA PLOTTED 

SHEET NO * 206 RUNS * 39, 59 , 0 , 0 , 0 , 0 , 0 , 0 , 


RESULTS OP 

LEAST SQUARES POLYNOMIAL CURVE 

FIT FOR RUN 

- 39 

NSYM • 1 

I 

T 

X 

XNEW 

X-XNEW 

Y 

YNEW 

Y-YNFW 

1 

-6.1700 

.9150 

.9150 

-.0000 

1.1329 

1.0588 

.0741 

2 

-4.1800 

1.9100 

1.9100 

.0000 

1.4455 

1.4068 

.0387 

3 

-1.8800 

3.0600 

3.0600 

.0000 

1.7651 

1.8479 

-.0828 

4 

.1700 

4.0850 

4.0850 

.0000 

2.2144 

2.2678 

-.0534 

5 

2.0100 

5.0050 

5.0050 

.0000 

2.6249 

2.6597 

-.0348 

6 

4.2800 

6.1400 

6.1400 

.0000 

3.1164 

3.1541 

-.0377 

7 

6.3700 

7.1850 

7.1850 

.0000 

3.5735 

3.6117 

-.0382 

8 

8.4500 

8.2250 

8.2250 

.0000 

4.0115 

4.0610 

-.0495 

9 

10.3200 

9.1600 

9.1600 

.0000 

4.4649 

4.4533 

.0116 

10 

12.4100 

10.2050 

10.2050 

0.0000 

4.9459 

4.8711 

.0747 

11 

14.4000 

11.2000 

11.2000 

-.0000 

5.4040 

5.2417 

.1623 

12 

16.8000 

12.4000 

12.4000 

-.0000 

5.8305 

5.6430 

.1874 

13 

18.6700 

13.3350 

13.3350 

-.0000 

5.9419 

5.9138 

.0281 

14 

20.8600 

14.4300 

14.4300 

-.0000 

6.0995 

6.1757 

-.0763 

15 

22.7500 

15.3750 

15.3750 

-.0000 

6.2170 

6.3472 

-.1301 

16 

24.7800 

16.3900 

16.3900 

-.0000 

6.3325 

6.4677 

-.1352 

17 

26.6700 

17.3350 

17.3350 

.0000 

6.4900 

6.5143 

-.0243 

IB 

28.7500 

18.3750 

18.3750 

.0000 

6.3650 

6.4852 

-.1202 

19 

30.7000 

19.3500 

19.3500 

.0000 

6.5800 

6.3744 

.2056 


ERROR 

IN X « 

•154875E— 22 


ERROR IN Y 

■ • 189547E+00 





ro 

ro 


TABLE VI.- SAMPLE OF LEAST-SQUARES CUBIC-SPLINE CURVE FIT 


DATA PLOTTED 








SHEET NO * 

206 RUNS * 39* 

59, 

0, 

0, 

0, 

% 

o 

% 

o 

«k 

o 

RESULTS OF 

LEAST SQUARES CUBIC SPLINE CURVE FIT FOR 

RUN « 39 

NS YM * 

1 DF - .10000 

I 

T 

X 

XNEW 

X-XNEV 

Y 

YNEW 

Y-YNEW 


1 

-6.1700 

.9150 

.9150 

0.0000 

1.1329 

1.0207 

.1122 


2 

-A. 1800 

1 .9100 

1.9100 

.0000 

1.AA55 

1.A16A 

.0291 


3 

-1.8800 

3.0600 

3.0600 

0.0000 

1.7651 

1.8793 

-.11A2 


A 

.1700 

A. 0850 

A. 0850 

0.0000 

2.21AA 

2.3021 

-.0877 



2.0100 

5.0050 

5.0050 

.0000 

2.62A9 

2.6911 

-.0662 


6 

A. 2800 

6.1 AOO 

6.1A00 

.0000 

3.116A 

3.1807 

— . 06 A3 


7 

6.3700 

7.1850 

7.1850 

.0000 

3.5735 

3 • 63 5 A 

-.0619 


e 

8.A500 

P.2250 

8.2250 

0.0000 

A. 0115 

A. 0837 

-.0722 


9 

10.3200 

9.1600 

9.1600 

0.0000 

A.A6A9 

A. A7AA 

-.0095 


10 

1 2 • A 1 0 0 

10.2050 

10.2050 

0.0000 

A.9A59 

A. 8855 

• 060 A 


n 

1A.A000 

11.2000 

11.2000 

.0000 

5.AOAO 

5.2396 

. 16AA 


12 

16.8000 

12.A000 

12.A000 

.0000 

5.8305 

5.6055 

.2250 


13 

IP. 6700 

13.3350 

13.3350 

.0000 

5.9A19 

5. 8 AOO 

.1016 


1A 

20.8600 

1A.A300 

1A.A300 

.0000 

6.0995 

6 • 062 A 

.0371 


15 

22.7500 

15.3750 

15.3750 

.0000 

6.2170 

6.2162 

.0009 


16 

2 A. 7800 

16.3900 

16.3900 

.0000 

6.3325 

6.3511 

-.0186 


17 

26.6700 

17.3350 

17.3350 

0.0000 

6.A900 

6.A569 

.0331 


18 

28.7500 

IP. 3750 

18.3750 

.0000 

6.3650 

6.5606 

-.1956 


19 

30.7000 

19.3500 

19.3500 

.0000 

6.5800 

6.6538 

-.0738 



ERROR 

IN X ■ 

• 575553E-25 


ERROR IN 

Y ■ • 190000E+00 





TABLE VII. - COMPUTER CODE FOR CASE 1 


C DEFINE HEIGHT OF COEFFICIENTS 

HT«0 • A 

C DEFINE HEIGHT OF AXIS LABELS 

SZ*0.75*HT 

C PROVIDE A BLANK SPACE AT BOTTOM OF PLOT FOR FIGURE TITLE 

CALL CALPLT <0.0*4.0#-3> 

C DRAW GRID LINES 

CALL GRIDLN < 0. 0# 0. 0* 20 . 0, 23 .0, 1 ) 

C DRAW ZERO LINES 

CALL ZEROLN ( 4 . 0# 0. 0* 23 . 0# 2 > 

CALL ZEROLN ( 0. 0 , 2 . 0* 20 . 0 * I ) 

CALL ZEROLN ( 0. 0* 1 1. 0* 20* 0# 1 ) 

CALL ZEROLN ( 0. 0, 20.0* 20.0# 1 ) 

C DRAW AND LABEL X-AXIS ( ALPHA > 

CALL AXISLB ( 1# 0 .Of 0. 0* -0 . 0# 2. 0# 20.0 # SZ*-1 ) 
X0«(20.0-3.5*HT) /2*0 
YO— 1.0 

CALL COEFSY (X0#Y0*HT,8) 

C DRAW AND LABEL LOWER PORTION OF Y-AXIS ( CL ) 

CALL AXISLB (2#0.0#0.0*-0.8#0.4#10.0#SZ#1) 

XO— 1.8 
Y0-5.0 

CALL COEFSY (X0,Y0#HT,1) 

C DRAW AND LABEL MIDDLE PORTION OF Y-AXIS ( CD ) 

CALL AXISLB (2# 0.0# 11 . 0* 0. 0# 0. 2# 5 .0# SZ# 1) 

XO— 1.8 
YO-13.5 

CALL COEFSY (X0,Y0*HT*2> 

C DRAW AND LABEL UPPER PORTION OF Y-AXIS ( CM ) 

CALL AVISLB ( 2. 0 . 0 # 1 7 . 0, — 1 . 2 * 0 . A. 6 . 0 » S 7 , 1 ) 

XO— 1.8 
Y0-20.0 

CALL COEFSY (X0#Y0#HT#3) 

C DRAW AND LABEL RUN KEY 

XD-2.0 
Y0-HT+23.0 

CALL RUNKFY < XO* YO* RUN# LSYM# ISI ZE#HT* 10) 


123 




non o o noon o o noon o o noon o o noon 


TABLE VIII.- COMPUTER CODE FOR CASE 2 


C DEFINE HEIGHT OF COEFFICIENTS 

HT-0.25 

C DEFINE HEIGHT OF AXIS LABELS 

SZ-0.75+HT 

C PROVIDE A BLANK SPACE AT BOTTOM OF PLOT FOR FIGURE TITLE 

CALL CALPLT (0.0, 3. 0,-3) 

LOWER LEFT-HAND PLOT ( CM VS CL ) 

DRAW GRID LINES AND ZERO LINES 
CALL GRIDLN (0.0,0. 0,5. 0, 5.0,1 ) 

CALL ZEROLN (1 . 0, 0 . 0, 5 . 0, 2 > 

CALL ZEPOLN (0.0,2.0,5.0,11 
DRAW ANO LABEL X-AXIS 
CALL AX IS LB ( 1, 0. 0, 0. 0, -0 . 8, 0. 8, 5 .0, SZ, 1 ) 

X0-(5.0— 2.0*HT) / 2 . 0 
Y0--0.8 

CALL CHEESY (X0,Y0,HT,1> 

DRAW AND LABEL Y-AXIS 
CALL AX IS LB (2, 0.0, 0.0, -0.8,0. A, 5.0, SZ,1) 

XO—l.O 
YO-2.5 

CALL COEFSY (X0,Y0,HT,3) 

LOWER RIGHT-HAND PLOT ( L/D VS ALPHA 1 

DRAW GRID LINES AND ZERO LINES 
CALL GRIDLN ( 7 . 0, 0. 0, 5 . 0, 5. 0, 1 ) 

CALL ZEROLN (8.0,0.0,5.0,21 
CALL ZEPOLN ( 7 . 0,1 .0, 5 . 0 , 1 ) 

DRAW ANO LABEL X-AXIS 

CALL AXISLB ( 1, 7.0, 0.0, -10 .0,10 .0,5 . 0, SZ,-1 ) 
XH-7.0M5.C-3.5 + HT 1/2.0 
YO — 0.8 

CALL COFFSY (X0,Y0»HT,8) 

DRAW AND LABEL Y-AXIS 

CALL AXISLB ( 2, 7 . 0, 0 . 0,-4 . 0, <♦ . 0, 5 . 0, SZ,-1 ) 

X0-6.0 
YO* 2. 5 

CALL COEFSY (X0,Y0,HT,7) 

UPPER LEFT-HAND PLOT ( CL VS ALPHA ) 

DRAW GRID LINES AND ZERO LINES 
CALL GRIDLN 10.0,7.0,5.0,5.0,1) 

CALL ZEROLN (1 . 0, 7.0, 5 . 0, 2 ) 

CALL ZEPOLN ( 0. 0, 8 . 0, 5 .0, 1 ) 

DRAW AND LABEL X-AXIS 

CALL AXISLB ( 1, 0. 0, 7. 0,-10. 0, 10. 0, 5. 0, SZ,-1 ) 
X0-(5.0-3.5*HT) /2.0 
YO-6.2 

CALL COFFSY (X0,Y0,HT,8) 

DRAW AND LABEL Y-AXIS 
CALL AXISLB ( 2, 0 .0 ,7. 0,-0 . 8, 0. 8 , 5 .0, S Z, 1 ) 

XO— 1.0 
YO-9.5 

CALL COEFSY (X0,Y0,HT,1> 

UPPER RIGHT-HAND PLOT ( CD VS CL ) 

DRAW GRID LINES AND ZERO LINES 
CALL GRIDLN (7.0,7.0,5.0,5.0,11 
CALL ZEROLN (8.0,7.0,5.0,21 
CALL ZEROLN (7.0,7.0,5.0,11 
DRAW AND LABEL X-AXIS 
CALL AXISLB (1, 7 . 0, 7. 0,-0 .8, 0. 8,5 .0, SZ, 1 1 
X0-7.0+(5.0-2.0*HT) /2.0 
YO-6.2 

CALL COEFSY (X0»Y0,HT»1) 

DRAW AND LABEL Y-AXIS 
CALL AXISLB ( 2, 7 . 0, 7. 0, 0 . 0, 0 . 2, 5 . 0, S Z, 1 1 
X0-6.0 
YO-9.5 

CALL COEFSY (X0,Y0,HT,2) 

DRAW AND LABEL RUN KEY 

XO-1.0 
Y0-12.0+HT 

CALL RUNKEY ( XD, YO, RUN, LS YM, ISI ZE, HT, 10 1 






RLPHfl 

2.0000 

TEST '496 .SHEET 206 .RUNS 39 59 0 0 0 


Figure 1. - Plot of Case 1 input with standard symbols and cubic-spline curve fairing. 




2.0000 


TEST 496 .SHEET 206 .RUNS 39 


59 0 0 0 0 


Figure 2.- Plot of Case 1 input with standard dashed line combination and cublc-spllne curve fairing 



ALPHA 

2.0000 

TEST 496 , SHEET 206 .RUNS 39 59 0 0 0. 0. 0 


a) Least-Squares Polynomial 


ng. 


Figure 3.- Plot of Case 1 input with least-squares curve fairi 


127 



2.0000 


TEST 436 -SHEET 206 -RUNS 33 53 

b) Least-Squares Cubic Spline 


Figure 3. - Concluded. 
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39 
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Figure 4.- Plot of Case 1 input with grid lines, axis scale values and labels, and legend drawn. 
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bels, and legend drawn. 
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